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This review discusses instabihties of the Fermi-Uquid state of conduction electrons in metals with 
particular emphasis on magnetic quantum critical points. Both the existing theoretical concepts 
and experimental data on selected materials are presented; with the aim of assessing the validity 
of presently available theory. After briefly recalling the fundamentals of Fermi-liquid theory, the 
local Fermi-liquid state in quantum impurity models and their lattice versions is described. Next, 
the scaling concepts applicable to quantum phase transitions are presented. The Hertz-Millis- 
Moriya theory of quantum phase transitions is described in detail. The breakdown of the latter 
is analyzed in several examples. In the final part experimental data on heavy-fermion materials 
and transition-metal alloys are reviewed and confronted with existing theory. 
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I. INTRODUCTION 

The Fermi-liquid description of metals is one of the 

most successful theories in condensed matter physics. It 
can be applied to describe vastly different systems, rang- 
ing from liquid ^He to simple metals like copper or gold to 
complicated compounds like CeCue, where the Coulomb 
interaction in strongly localized /-electron shells leads to 
gigantic; interaction effects and a hundred-fold increase of 
the effective masses. Deviations from Fermi-liquid behav- 
ior are a central topic in the experimental and theoreti- 
cal studies of correlated electronic systems, triggered by 
the discovery of high-temperature superconductivity, the 
success in synthesizing effectively low-dimensional mate- 
rials, and the study of compounds which can be tuned 
through zero-temperature phase transitions. 



A. Outline and scope of the review 

In this review we want to give a combined theoretical 
and experimental overview of the breakdown of Fermi- 
liquid (FL) behavior in the vicinity of magnetic quantum 
phase transitions. After a summary of Fermi-liquid the- 
ory (Sec. II), including the Kondo effect in local-moment 
and Kondo-lattice systems, we will describe in Sec. Ill 
the established theoretical approach to continuous quan- 
tum phase transitions (QPT) in metallic systems. This 
approach was pioneered by Hertz (1976). Recent theories 
have addressed its inadequacy in a number of important 
situations, and we will summarize the current status. In 
Sec. IV, we will then turn to a variety of experimental 
systems where quantum criticality and non-Fermi-liquid 
(NFL) behavior have been observed. By carefully exam- 
ining available results we shall attempt to state where 
the standard approach of Hertz applies, and for which 
systems other theories have to be considered. 

Due to space restrictions, a number of interesting top- 
ics in the field of metallic quantum criticality will be 
omitted. We will almost exclusively concentrate on three- 
dimensional (3d) metals, i.e., we will not touch upon 
high-temperature superconductors and other quasi-two- 
dimensional (2d) and quasi-one-dimensional (Id) mate- 
rials. We will only focus on QPT involving magnetic or- 
der, this removes genuine metal-insulator transitions as 
well as charge-density wave transitions from our agenda. 
Most of our discussion will be restricted to the paramag- 
netic and quantum critical regimes of the transitions; we 
will say little about the long-range ordered phases which 
pose additional complications (like e.g. non-trivial Gold- 
stone modes). Further, our primary interest is in clean 
materials where the effect of quenched disorder is weak. 
We will therefore leave out metallic spin glasses, and will 
only briefly mention Kondo disorder and quantum Grif- 
fiths scenarios as sources of NFL behavior - we refer the 
reader to recent reviews (Miranda and Dobrosavljevic, 
2005; Vojta, 2006b). A comprehensive compilation of ex- 
perimental NFL data was given by Stewart (2001, 2006). 



B. Non-Fermi-liquid behavior vs. breakdown of the 
Fermi-liquid concept 

Before we focus on some theoretical models and ex- 
perimental systems where the Fcrmi-liquid phenomenol- 
ogy appears to fail, it is important to discuss the some- 
times confusing terminology in this field. The meaning 
of phrases like "non- Fermi liquid" , "Fermi-liquid instabil- 
ity" or "breakdown of Fermi-liquid theory" varies quite 
substantially depending on the context, the community, 
or on theoretical prejudices. This review will not be able 
to avoid this problem completely, especially as we try to 
emphasize the open questions in this field. 

Conceptually, one should carefully distinguish between 
two quite different statements. The first is the observa- 
tion of "non-Fermi liquid behavior", i.e., apparent devi- 
ations from the Fermi- liquid phenomenology, e.g., from 
a constant specific-heat coefficient or a dependence 
of the resistivity at low temperatures. This experimen- 
tal definition should not be confused with the theoreti- 
cal statement of a "breakdown of Fermi-liquid theory" , 
which implies that the concept of a FL and its underlying 
assumptions (see Sec. II) have become invalid. 

From this point of view it is not surprising that some- 
times "non-Fermi-liquid behavior" can be explained us- 
ing FL concepts. For example, in a disordered FL the 
low-temperature resistivity displays a \fT cusp (Alt- 
shuler and Aronov, 1985) instead of the quadratic T de- 
pendence of the weakly disordered case. While this non- 
analytic behavior is related to the existence of diffusive 
modes in the disordered system, which are absent in the 
usual phenomenology of a Fermi liquid, this does not im- 
ply a "complete breakdown" of the FL concept. Another 
example is the theory of Hertz (1976) for magnetic quan- 
tum phase transitions in three dimensions (see Sec. III.C) 
where the relevant low-energy excitations are the usual 
fermionic quasiparticles and their collective excitations. 
The magnetic collective excitations become soft at the 
quantum critical point (QCP), mediating a singular in- 
teraction between the quasiparticles and therefore induc- 
ing NFL behavior. Nevertheless, methods and concepts 
of FL theory can still be applied - this underlying as- 
sumption is actually the basis of this theory. 

What are the theoretical concepts that can replace the 
Fermi-liquid paradigm in cases where Fermi-liquid the- 
ory breaks down and the low-energy excitations do not 
carry the quantum numbers of fermionic quasiparticles? 
In many cases, the answer to this question is not known. 
One can envision at least two possible scenarios. One is 
that new weakly interacting quasiparticles with different 
quantum numbers and interactions can be found. Fa- 
mous examples are Luttinger liquids with purely bosonic 
excitations, or the fractional quantum Hall effect with 
quasiparticles with fractional charge. The other possi- 
bility is more difficult to treat theoretically, and much 
less is known in this case: it is conceivable that no well- 
defined quasiparticles exist at all and all excitations are 
incoherent. 
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C. Exponent puzzles 

Non-Fermi-liquid behavior manifests itself in the 

power-law behavior of physical quantities, with expo- 
nents different from those of a Fermi liquid. In this sec- 
tion we discuss under what conditions power-law behav- 
ior is expected and we comment on the problem of how 
exponents can be extracted from experimental data. 

Power-law behavior of physical quantities is generally 
expected in the absence of any close-by scale. For exam- 
ple in an ordinary metal one type of power laws - e.g. 
a linear temperature dependence of the resistivity p(T) 

- is usually observed in the phonon- dominated regime 
u>D <^T <^ ep, where ojd stands for the Debye frequency, 
and ep is the Fermi energy (ks = h = 1). A different set 
of Fermi-liquid exponents governs physical quantities in 
the regime T ^ ujd,€f, with, e.g., Ap = p{T) — po ^ 

or Ap ~ in regimes dominated by electron-electron 
or electron-phonon scattering, respectively. 

While it is sometimes possible to fit a broad crossover 
regime (e.g., aroimd T k uju) with some effective power 
law, one should not confuse this non-universal and non- 
generic crossover effects with a true power-law scaling be- 
havior. From a theoretical point of view, an exponent is 
well defined, if at least formally a scaling limit exists (e.g., 
T /ljd oo and T /ep — > 0) where the power law can be 
observed in a broad temperature range. This purely for- 
mal condition translates to the experimental requirement 
that an algebraic behavior can only be established if it 
extends over a considerable range of, e.g., temperature 
in a regime where no other relevant scale is expected to 
exist. Note that it is not required that the regime of 
power-law behavior extends down to zero temperature. 
A Fermi liquid, for example, is essentially never the true 
ground state of a metallic system, the Fermi-liquid fixed- 
point can nevertheless govern the physics over several 
decades in temperature down to an exponentially small 
temperature where, e.g., superconductivity sets in. 

From an experimental point of view, there are various 
methods to extract exponents from a measurement. As 
an illustration for the difficulties which arise and for some 
of the methods used in this context, we discuss briefly 
how, e.g., an exponent characterizing the T dependence 
of the resistivity p{T) can be determined. Most com- 
monly, the exponent is obtained from a fit to the data, 
e.g., the T dependence of the resistivity is fitted to a 
model form p(T) = po + AT'^ for a certain temperature 
range Tmin <T < Tmax- A general criterion under which 
conditions such a procedure is reliable cannot be given 

- this depends crucially on both prefactor and T depen- 
dence of the leading corrections to the algebraic behavior 
which in many cases are not known. Usually, such fits 
are believed to be reliable if the exponent depends in a 
certain regime only weakly on Tmax and Tmin, if the fit 
extends over more than one decade in temperature and 
if a plot of p(T) as a function of T" "looks" linear. 

A slightly less biased method which can also be used 
to investigate crossover phenomena, is the calculation of 



an effective T-dependent exponent, defined by the log- 
arithmic derivative of the measured quantity, a(T) « 
din[p{T) ~ po]/dinT. One should keep in mind that, 
in particular in crossover regimes, this effective exponent 
may have little physical significance. However, if the data 
really can be described by a power law in a large tem- 
perature regime, then a{T) will be independent of T and 
coincides with the "true" exponent. A serious problem is 
the dependence of this procedure on the residual resistiv- 
ity po- It is usually chosen in such a way as to get the least 
temperature dependence of a{T) - this is a dangerous 
bias in the interpretation of the data, especially if some 
deviations from power-law behavior can be expected, e.g., 
close to but not directly at a quantum phase transition. 
This bias can be avoided by defining a temperature- 
dependent exponent a{T) = l + Mn[dp{T)ldT]/MnT, or 
equivalently, by fitting power laws to the data in small 
temperature intervals. Obviously, the latter definition of 
a{T) does not depend on a constant background po but 
it is numerically very unstable and requires rather precise 
data. Hence this procedure is rarely employed. 

All of the above methods fail if a small scale A exists 
where the behavior of the measured quantity crosses over 
from one to another power law. This is a very common 
situation close to some quantum critical point, where A 
can be tuned as a function of some control parameter like 
pressure p, magnetic field B or amount of disorder. In 
this situation a critical behavior of some quantity, e.g. 
X{T) = p{T) — Po, is expected, which has the scaling 
form X{T,p) - r"/(T/A(p)) with A cx (p-pc)"' in the 
limit T, A ^ 0, T/A = const (see Sec. III.B for more 
details and the discussion of other quantities). In this 
situation, a scaling analysis is the ideal tool to extract 
the exponents a and a' and the asymptotics in the limit 
T/A > 1 and T/A < 1. If scaling holds, the data 
for various values of the control parameter p can be col- 
lapsed on a single curve by plotting X/T" as a function 
oiT / [p—pc)"' . This scaling collapse is used to determine 
the exponents a and a' . However, meaningful scaling re- 
quires that the data from different curves to be collapsed 
do overlap significantly - a condition not always fulfilled. 



D. Non-Fermi-liquid scenarios 

In a well-defined Fermi liquid, the usual FL expo- 
nents show up below a characteristic scale T*. (T* is 
non-universal and depends on many parameters, e.g., 
the strength of electron-phonon interactions.) What are 
the requirements to observe non-Fermi-liquid exponents 
down to the lowest temperatures? A trivial answer to 
this question is that the scale T* has to disappear. This 
can happen in at least two ways: (i) T* may be tuned to 
zero, e.g., by approaching a QCP; (ii) T* may be elimi- 
nated by strong disorder: if the distribution of T* in the 
system is sufficiently broad, no characteristic energy can 
be defined below which macroscopic Fermi-liquid theory 
is valid and NFL behavior is expected. 
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Possibility (i), namely the suppression of T* in the 
vicinity of a magnetic bulk quantum critical point, is the 
main topic of this review. In Sec. II. F we will also men- 
tion single-impurity critical points which can induce local 
NFL behavior. Possibility (ii) is covered by the recent re- 
view of Miranda and Dobrosavljevic (2005). 

We note a further route to NFL behavior here: in prin- 
ciple, a stable NFL fixed point (corresponding to a NFL 
phase) may exist, where the low-energy excitations do 
not carry the quantum numbers of fermionic quasiparti- 
cles. However, we are not aware of any promising can- 
didate for such a fixed point of a metal in d = 3 (with 
the exception of gauge field theories, see Holstein, Nor- 
ton, and Pincus 1973; Varma, Nussinov, and van Saarloos 
2002, or metallic states with liquid-crystal-like order, see 
Oganesyan, Kivelson, and Fradkin 2001). An extensive 
discussion about NFL fixed points, especially in d = 2, 
can be found in the context of theories of NFL behavior 
in high-temperature superconductors, see e.g. Anderson 
(1997); Varma, Nussinov, and van Saarloos (2002). Ex- 
perimentally, MnSi under pressure shows signatures of a 
NFL phase md = 3, see Sec. IV. B. 2, however, a theoret- 
ical description is not available to date. 



II. LANDAU FERMI-LIQUID THEORY 

Systems of interacting fermions at low temperature 
have been of interest early on in the development of con- 
densed matter theory. Landau put forth a phenomeno- 
logical theory of interacting Fermi systems, the Fermi- 
liquid theory or Landau theory, which is based on the 
concept of quasiparticles (Landau, 1957a, b, 1959). It 
proposed to map the properties of Fermi systems at low 
temperature T onto a dilute gas of strongly interacting 
fermionic excitations. To some extent, a microscopic jus- 
tification of this picture was given by Landau and oth- 
ers, although a rigorous general mathematical proof is 
not available. Recent studies of this problem have used 
the renormalization group (RG) method (Feldman et al., 
1993; Shankar, 1994) which can be used to establish rig- 
orous mathematical bounds on the stability of the FL 
state (Feldman et al., 1993). 

In the following we review the salient assumptions and 
results of FL theory for later reference. 



A. The quasiparticle concept 

Let us start by considering the non-interacting system, 
where the occupation of single-particle states |kcr) with 
momentum k is given by 



T=0 
kcr 



e{kF - k) 



(1) 



where 6{x) is the step function. The Fermi momen- 
tum kp is determined by the density of particles n = 



Let us now imagine that the interac- 



the low-energy excitation spectrum of the interacting sys- 
tem is in one-to-one correspondence with the Fermi-gas 
spectrum, and if the ground state retains the full symme- 
try of the Hamiltonian, the system is termed a "normal 
Fermi liquid" . Note that the interaction will lead to the 
appearance of collective modes. However, these bosonic 
excitations occupy a negligible fraction of phase space in 
the limit of low temperatures and therefore do not spoil 
the principal one-to-one correspondence of single-particle 
states. In an ordered state this one-to-one correspon- 
dence is lost. 

The low-energy single-particle excitations of the Fermi 
liquid, with quantum numbers k and a, are called "quasi- 
particles". In the ground state, their distribution func- 
tion is again riko- (1). The energy of a quasiparticle, eto-, 
is defined as the amount of energy by which the total en- 
ergy E of the system increases, if a quasiparticle is added 
to the unoccupied state Ikcr), 



Ekcr = 



dE 



(2) 



where dnua is the corresponding change of the distri- 
bution function. As a consequence of the interaction, 
the single-particle energies depend on the state of the 
system, Cka = ekalnkv}- The energy of a single low- 
energy quasiparticle added to the ground state may be 
parametrized as 



ekal^k'a' } = 1>F{k - kp) 



(3) 



for an isotropic system at small energies, with vp = ^ 

being the Fermi velocity. The effective mass m* deter- 
mines the density of states per spin at the Fermi level 



No 



m*kF 



(4) 



(here and in the following we use units where h = ks = ^, 
imless explicitly stated). The effect of interactions with 
other excited quasiparticles on the energy of a specific 
quasiparticle may be expressed in terms of an effective 
two-particle interaction function or "Fermi-liquid" inter- 
action /ko-k'o-' 



kV 



(5) 



kcr- 



where Jnto- = Hko- — n\ 

For isotropic systems with short-range interaction the 
FL interaction function only depends on the angle be- 
tween k and k' and on the relative spin orientation of a 
and a', and hence may be parametrized as 



° (.=0 



(6) 



tion between the particles is turned on adiabatically. If 



Here k = k/|k|; cr = ±1, Pi{x) arc the Legen- 
dre polynomials, and and F^ are the dimensionless 
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spin-symmetric and spin-antisymmctric "Landau param- 
eters" , which characterize the effect of the interaction on 
the quasiparticle energy spectrum. For Galileian invari- 
ant systems the Landau parameter and the effective 
mass m* are related through m* /m = 1-1- /3. 

In a crystal, the symmetry of the system is reduced to 
discrete rotations/reflections (the elements of the point 
group of the lattice), and (if spin-orbit interactions can 
be neglected) rotations in spin space. As a conseqiience 
the band structure ek and the FL interaction /ko-kv may 
be strongly anisotropic. The parametrization of ek and of 
/ktrkv then requires additional parameters, which weak- 
ens the predictive power of FL theory. In applications of 
FL theory to metals, it is frequently assumed that an 
isotropic approximation in 3d or quasi-2d systems can 
give a reasonable account of the FL properties. 



B. Thermodynamic properties 

The equilibrium distribution function n^^ at finite 
temperature T follows from the assumed one-to-one cor- 
respondence: 



In the spin-symmetric isotropic case the compressibility 
dn/dfi diverges when —1, which is an indication of 

phase separation into a dense and a dilute phase. More 
common is the case of ferromagnetism, which appears 
when Fq ^ — 1. In the case of an instability at ^ > 
the corresponding susceptibility of an anisotropic density 
excitation in k-space diverges, which is termed Pomer- 
anchuk instability. It may lead to an anisotropic defor- 
mation of the Fermi surface (for the spin-antisymmetric 
sector this has been considered by Akhiezer and Chud- 
novskii, 1976). While spatially uniform Fermi surface de- 
formations can be captured by Fermi-liquid theory, this 
is more difficult in the case of instabilities (charge and 
spin density waves) at finite momentum, as the full mo- 
mentum dependence of fkk'{q) becomes important. The 
critical behavior of the FL properties on approaching a 
Pomeranchuk instability is discussed in Sec. III.G. 

A different class of instabilities is signaled by a singu- 
larity in the quasiparticle scattering amplitude at zero 
total momentum. It usually leads to the formation of 
pair-correlated ordered states, i.e., unconventional super- 
conductors, which are not a subject of this review. 
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(7) 



This is a complicated implicit equation for n^^^ due to 
the dependence of eka- on {n'^,^i}- 

The derivative of the internal energy with respect to 
temperature yields the specific heat at constant volume. 
The leading term at T ^ Tp {Tp = is the Fermi 
temperature) is linear in T, as for the free Fermi gas, 
and given by the (renormalized) density of states 



-NoT = jT. 



(8) 



The spin susceptibility x at !> and the electronic 
compressibility follow as 



X 



1 



dn 
d/j, 



2No 



(9) 



where /x^ is the magnetic moment of electrons, x and 
are affected both by the mass renormalization and by 
Fermi liquid parameters describing an effective "screen- 
ing" of the external fields. 



C. Instabilities within a Fermi-liquid description 

Thermodynamic stability requires that the susceptibil- 
ities X and dn/dji be positive, which leads to the require- 
ments F^'* > —1. A general analysis of the stability of 
the system with respect to any variation of rika results 
in the stability conditions (Pomeranchuk, 1958) 



>-(2£+l) , £ = 0,1,. 



(10) 



D. Finite-temperature corrections to Fermi-liquid theory 

The leading corrections to Fermi-liquid theory at low 

temperatures T <§; Tp have been the subject of extensive 
theoretical and experimental study. While for a Fermi 
gas these corrections are of relative magnitude (T/Tp)^, 
collective effects in an interacting Fermi system generally 
lead to much larger corrections. Thus one finds that the 
specific-heat coefficient varies with temperature as 



7(T) - 7(0) 



-g^THniT) for d = 3 
-g2T for d = 2 



(11) 



The coefficients 92 and 93 have been calculated exactly 
in lowest order perturbation theory (for a review and 
references see Chubukov, Maslov, and Millis 2006). Sub- 
leading corrections for a Fermi gas with weak repulsion 
have been derived within an RG approach by Aleiner 
and Efetov (2006). The leading corrections to the spin 
susceptibility have been found as 



X{T) - X(0) 



-C3T2 for d = 3 
-C2T for d = 2 



(12) 



In contrast to the specific heat, in d = 3 a non-analytic 
contribution is absent (Belitz, Kirkpatrick, and Vojta, 
1997; Chubukov, Maslov, and Millis, 2006). There is, 
however, a non-analytic dependence on the wavevector 
(see Sec. III.H.l) and on magnetic field. In two dimen- 
sions the leading T-power of the correction is again re- 
duced from to T by singular interaction processes 
(Chubukov et al, 2005). 

Experimentally, the best evidence for the above finite- 
temperature corrections has been reported in ^He, where 
a T'^£n{T) contribution in 7(T) could be identified (see 
Grey wall 1983 and references therein). 



6 



E. Transport properties 

1. Quasiparticle relaxation rate 

At low temperature, T <C Tp, there exists a small num- 
ber of thermally excited quasiparticles, which interact 
strongly. The decay rate of a quasiparticle on top of 
the filled Fermi sea is dominated by binary collision pro- 
cesses: the considered quasiparticle in state |1) = |kic7i) 
scatters off a partner in state |2), the two particles ending 
up in final states |3) and |4). The decay rate is given by 
the golden rule expression 



1 



27r ^' |a(l, 2; 3, 4)pn°(l - n°)(l - nl) (13) 



234 



where a(l,2;3,4) is the transition amplitude. The sum- 
mation over momenta and spins is restricted by conser- 
vation of momentum, energy, and spin. A full evaluation 
of yields (Baym and Pethick, 1991) 



1 

Tk 



{W) = dcos 



277 



(14) 



(|Ao(0,0)P+3|Ai(0, 



The quantities Aq and Ai arc the dimensionless scatter- 
ing amplitudes in the singlet and triplet channel [>lo,i = 
2Aroa(l, 2; 3, 4)], 6 and </> parametrize the angle between 
ki, k2 and the planes (ki,k2), (k3,k4), respectively. In 
2d systems the prefactor of in t^^ is logarithmically 
enhanced, ^ ~ T'^eniTp/T) (Chubukov et al, 2005). 

The forward scattering limit of the quasiparticle scat- 
tering amplitude can be expressed as (Landau, 1959) 



A-{e,ct> = o) = J2Y 



+ F//(2£+l) 



Pe{cos0) (15) 



where a = s,a labels the spin symmetric and antisym- 
metric particle hole channels, respectively. As the sys- 
tem approaches a phase transition to a state governed 
by spatially uniform order, such as a ferromagnet, the 
corresponding component of A" {9,0) tends to diverge. 
In the case of the ferromagnet, —1, and A^ = 

Fq/{1 + Fq) 00. The quasiparticle scattering is then 
dominated by ferromagnetic fluctuations, and the relax- 
ation rate r^^ is expected to scale as Nq/{1 + FqY <^ X^i 
with X the static spin susceptibility. 



momentum of the quasiparticles, i.e., lu ^ T, q ^ T/vp. 
The distribution function satisfies the kinetic equation 

StnicCT + VfcCkcr • V^nka - VrEkcr • VfeUkcT = I{nk<T} (16) 

The left-hand side describes the dissipationless flow of 
quasiparticles in phase space. It goes beyond the Boltz- 
mann equation in that the quasiparticle energy ek(r, t) it- 
self depends on position and time, due to its dependence 
on the distribution function as given by (5). Among 
other things, this gives rise to the appearance of collec- 
tive modes, as well as interesting nonlinear effects (which 
we will not discuss). 

On the r.h.s. of (16) we have the so-called collision in- 
tegral /, which describes the abrupt change of momentum 
and spin of quasiparticles in a collision process. It is given 
by Ik. = -nkjr™({nk.}) + (1 - r7.k,)/T™({l - 
n-ka}): which is the number of quasiholos minus the num- 
ber of quasiparticles in state \ka) decaying per unit time. 
The non-equilibrium relaxation rate l/'''^°"°'^({'T-k<7}) is 



and ek 



by their 



obtained from (13) by replacing nk^, 
non-equilibrium counterparts. 

If the applied external field is weak and one is allowed 
to linearize in the deviation of the distribution function 
from its equilibrium value. The resulting linearized and 
Fourier-transformed kinetic equation is given by 



(w - Vk • q)6nk{q, w) -h Vk • q ^-^<^ek 

CCk 



: i6I (17) 



where (5ek<j (q, i^) is defined as 



hak'a'Snu'a' (q, w) . 



(18) 



In principle, /ko-kv is also a function of q and w, and 
the limit q, a; — > is understood here. For a charged 
system with long-range Coulomb interaction, a classical 
or Hartree-type interaction part must be separated out, 

fkaka' (q) — + fkak'a' ■ 

Assuming the Landau parameter Fq to be dominant 
and collision effects to be small, the kinetic equation may 
be solved in the presence of an external potential —dfjf^^^ 
to obtain the density response function: 



5n 

—^=xM,oj) 



xo(q,t^) 



where 



Xo(q,w) = ^ 

ko- 



l + ^;fXo(q,a;) 
kq 



w — Vk • q -I- iO 9ek 



(19) 



(20) 



2. Kinetic equation and dynamic response 

In the presence of slowly varying disturbances, the sys- 
tem may be described by a quasiclassical distribution 
function n)s_„{T,t). This is possible as long as the en- 
ergy and momentum of the quanta of the external field, 
uj and q, are much smaller than the typical energy and 



This is the well-known random-phase approximation 
(RPA) form of the density response function in the limit 

q ^ kp. 

In analogy to the density response, the dynamical spin 
susceptibility Xmi^.^^) is defined as the response of the 

magnetization to a magnetic field, Xm = SM/6{iJ,rnB), 
and is given by Eq. (19) with Fq replaced by Fq. The 
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dynamical structure factor S'(q, uj), experimentally acces- 
sible via magnetic neutron scattering, is related to the 
dynamical susceptibility Xmi^.,^) through 

^(q,^) = [l+n(w)]lmxm(q,w + i0) (21) 

where n{u!) = {e'^^'^ — 1)~^ is the Bose function. 



A oc Nq and 7 cx A^o • This is indeed observed for a large 
number of heavy fermion systems (Kadowaki and Woods, 
1986), and Af-y"^ is termed the Kadowaki- Woods ratio. A 
corresponding dimensionless quantity may be defined as 



Rkw 



p{T)-po 



(25) 



3. Electrical resistivity 



F. Kondo effect - concept of a local Fermi liquid 



As an example of a transport coefficient we consider 
the electrical conductivity a = fr'^ (where p is the resis- 
tivity) , defined as the response of the electric current den- 
sity j to the (screened) electric field E, j = ctE (assuming 
cubic symmetry). In terms of the linearized distribution 
fiuiction fck — limoj^o hniq^o '5'T.k(q, w) the conductiv- 
ity is given by Uij = e X^ko- ^ki'^'T-kZ-Ej. Here 5ny^ satisfies 
the Boltzmann equation 



evk-E(g)=57{Jnl} 



(22) 



where is the deviation of Uk from local equilibrium. 
The collision integral 51 describes the effect of inter- 
quasiparticle collisions and any other collision processes. 
For a translation invariant system, the quasiparticle colli- 
sions arc momentum conserving and the resistivity would 
be zero. The resistivity is finite in a real solid, if Umk- 
lapp scattering is possible (for a recent discussion of the 
role of Umklapp scattering see Rosch 2005). The most 
important source of momentum dissipation at low T is, 
however, impurity scattering. 

In lowest approximation the collision integral may be 
modelled as 



51 



e5ni 



CiSni 



(23) 



where C^^e and Ci are the linear integral operators 
describing electron-electron and electron-impurity col- 
lisions, respectively. Taking into account that Cg-e ^ 
^ at low T, one finds for the resistivity 



piT)=po + AT^ + .. 



(24) 



Here po is the residual resistivity from impurity scatter- 
ing. The coefficient A is given by a weighted angular av- 
erage of the squared quasiparticle scattering amplitudes 
Ao,i{0,^) = 2A^oa(l, 2; 3, 4), which sensitively depends 
on the anisotropy of the scattering and the band struc- 
ture. The resulting transport scattering rate is in 
general different from the single-particle relaxation rate 
(14) (except for isotropic scattering), as only finite- 
angle scattering affects transport. 

Provided that the transition amplitudes a(l, 2; 3, 4) are 
weakly momentum dependent, i.e., when they are gov- 
erned by local physics, the ratio of the resistivity co- 
efficient and the square of the specific-heat coefficient, 
may be expected to be material-independent since 



The Kondo problem (Hewson, 1993) goes back to the 
discovery of a resistance minimum at low temperatures 

in metals with dilute magnetic impurities. The min- 
imum and low-T increase of the resistance were suc- 
cessfully explained by Kondo (1964) within a pertur- 

bative calculation. Within the so-called Kondo (or s- 
d) model, a magnetic impurity is described by a lo- 
cal spin S (assumed to be 5 i located at r = 0) 
exchange-coupled to the local conduction-electron spin 



density Sq = ^ Ek,k' E, 



Ck', 



Hsd = ^ ek4^Ck<7 + JS ■ 



So 



(26) 



where Taa' is the vector of Pauli matrices, and J is ex- 
change coupling. Kondo found that the electrical resis- 
tivity p due to scattering of conduction electrons off the 
impurity acquired a logarithmic dependence on temper- 
ature in third order in J 



Pb 



1 + 2iVo J in^ 



(27) 



where oc is the usual Born approximation result. 
No is the local conduction electron density of states per 
spin, and D is the half-width of the conduction band. 
The reason for the T dependence in p lies in the resonant 
scattering from the degenerate ground state of the mag- 
netic impurity. As seen from the result (27), perturbation 
theory breaks down below the Kondo temperature 



exp 



1 



(28) 



when the first-order correction term becomes comparable 
to the Born approximation. The prefactor in (28) is cho- 
sen such that a Xj lr?{^ jTy^ correction to the magnetic 
susceptibility ximp 

(T) is absent for T~>Tk- For J ^ 
the Kondo temperature Tk depends on the coupling J in 
a non-analytic way. 

As pointed out by Anderson and collaborators, the 
ground state of the Kondo model (26) is non-degenerate. 
The S = \ impurity spin is fully compensated by a 
"screening cloud" of conduction electron spins contain- 
ing in total one electron spin, bound to the impurity in a 
singlet state. The impurity complex formed in this way 
acts like a potential scatterer. The low-energy physics of 
this system has been formulated by Nozieres (1974) in 
terms of a local Fermi liquid picture. 
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It is \iscful to introduce a model of a magnetic impu- 
rity in which the electronic structure of the impurity is 
displayed directly. This Anderson impurity model (An- 
derson, 1961) consists of an impurity orbital (a d or / 
orbital of an incompletely filled inner atomic shell) hy- 
bridizing with a conduction band, 

k,(T CT 

+VY.^<iah + h.c.) (29) 

where /j! creates an electron with spin projection a in 
the / orbital, n/o- = flfa-, and V is the hybridization 
matrix element. The decisive feature of the model is 
the strong Coulomb interaction U which leads to local- 
moment formation: if the impurity level lies below the 
Fermi energy, ef < 0, while ef + U > and the bare 
hybridization width r = nNoV'^ is small, T <C \ef\,ef+U, 
the impurity level is mainly occupied by a single electron 
(rather than being empty or doubly occupied) and thus 
represents a local moment of spin i. 

By projecting the Anderson Hamiltonian (29) onto the 
subspace of singly occupied impurity states (Schrieffer 
and Wolff, 1966) one is led to the Kondo model with 



2V^ 



1 



1 



> 0. 



(30) 



We note that Kondo physics is not restricted to the 
screening of magnetic degrees of freedom, it can occur 
in any situation where transitions between a multi-level 
impurity are induced through the interaction with a bath 
of fermionic particles. The particular example of the 
quadrupolar Kondo effect will be discussed in Sec. II. F. 3. 



1. The Anderson or Kondo impurity 

The low-T FL regime of the exactly screened Ander- 
son model can be expected to be continuously connected 
to the non-interacting limit of the Anderson model, and 
consequently perturbation theory in powers of U is ap- 
propriate (Hewson, 1993). The dynamical properties of 
the impurity are described by the self-energy T,f^{Lv) of 
the local / electron, which determines the energy shift 
and broadening of the poles of the local Green's function 



Gf{u> + iO) = 



ef +iT -Hfcriio + iO) 



(31) 



where a particle-hole symmetric conduction band has 
been assumed. In the absence of Coulomb interaction, 
[/ — !■ 0, T,f^(uj) = 0. The imaginary part of S/cr measures 
the decay rate of the / electron into particle hole pairs: 
Standard phase-space restrictions dictate the small-w be- 
havior, 



It follows from the Kramers-Kronig relation that 

ReHf^oj) = Ef^iO) + (1 - Z-^)u} + 0{u}^) . (33) 

Thus, the Green's function Gfa in the limit of small oj 
takes the form 



uj-e*j^ + iT* 



(34) 



where e}^ = Z(e/ -h 5]/„(0)) and F* = ZT. Here Z plays 
the role of the quasiparticle weight factor in FL theory. 

The above analysis was restricted to the low-T regime, 
assuming a FL ground state. To capture the crossover 
from high to low temperature, the concept of scaling and 
of the renormalization group (RG) has turned out to be 
the principal tool. In Anderson's "poor man's scaling" 
approach (1970), applied to the Kondo model, the Hamil- 
tonian is projected onto a smaller Hilbert space where the 
conduction electron bandwidth Do has been reduced by 
an infinitesimal amount dD. The effective Hamiltonian 
still takes the form of the Kondo model, with a modi- 
fied coupling constant J{D). (Technically, "poor man's 
scaling" is a momentum-shell RG without the step of 
rescaling the cutoff and fields.) The flow of the dimen- 
sionless coupling j = NqJ is described by the so-called (3 
function, which can be calculated in a power series in j: 



dinD T 



(35) 



— ImS j„ {ui + iO) oc 



(32) 



The radius of convergence of this series is not known a 
priori. However, from (35) we can infer that for small an- 
tiferromagnetic j > the effective coupling grows under 
the reduction of the bandwidth. In contrast, a ferro- 
magnetic coupling {j < 0), j scales to zero for _D ^ 0, 
leaving a free spin and a free conduction band. This 
solves the problem for ferromagnetic coupling; the finite- 
temperature thermodynamics displays logarithmic cor- 
rections to the free-spin behavior due to the flow of j to 
zero. For the antiferromagnetic case, which is the usual 
situation, there are obviously two possibilities: (i) the 
coupling tends to a finite value as D ^ 0, (ii) the cou- 
pling grows indefinitely, j oo. It was recognized early 
on by Anderson and coworkers (Anderson, 1970; Ander- 
son and Yuval, 1970) that the latter scenario is realized 
in the usual (single-channel) Kondo situation. 

Using the fact that the fixed point of j — > oo is in- 
deed stable (Nozicrcs, 1974), one may now discuss the 
local properties of conduction electrons at the impurity, 
taking into account (i) the energy shift provided by the 
scattering off the static potential of the fully screened 
impurity, (ii) the interaction between conduction elec- 
trons at the impurity site induced by virtual excitations 
of the impurity complex. This is elegantly done in terms 
of a phenomenological description following the Landau 
FL theory (Nozieres, 1974), which provides a correct and 
useful picture of the Kondo behavior at small temper- 
atures and fields. In order to describe the crossover 
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from high-temperature local-moment physics to the low- 
temperatme Kondo-screened state, more elaborate meth- 
ods are required. The thermodynamic properties of the 
Kondo model (26) and the Anderson model (29) may 
be calculated exactly by using the Bethe ansatz method 
(Andrei, 1980; Andrei, Furaya, and Lowenstein, 1983; 
Wiegmann, 1980). The dependence of the free energy 
and its derivatives on temperature and magnetic field 
is found to exhibit single parameter scaling behavior in 
(T/Tk) and (S/Tk) from the low-energy regime below 
Tk up to high energies » Tk, but sufficiently below the 
conduction electron bandwidth. Most importantly, there 
is no phase transition upon variation of T or i? for the 
metallic single-impurity Kondo problem (see Sec. II. F. 4 
for phase transitions in impurity models). 

Dynamical properties can be calculated numerically 
using Wilson's renormalization group method (NRG) 
(Krishna-murthy, Wilkins, and Wilson, 1980; Wilson, 
1975). The local /-electron spectral function Af{uj) 
has been determined by NRG as well, both at T = 
(Costi and Hewson, 1990; Sakai, Shimizu, and Kasuya, 
1989) and at finite T (Costi, Kroha, and Wolfle, 1996), 
and by self-consistent diagrammatic methods (Kroha and 
Wolfle, 2005). The electrical resistivity for the Anderson 
model was found to obey single-parameter scaling behav- 
ior in T/Tk (Costi and Hewson, 1992). 

2. Anderson and Kondo lattice models 

So far wc have considered the properties of a single 
quantum impurity in a host metal. We now turn to a 
discussion of materials where "quantum impurity ions" 
are put on a lattice. Generalizing the single-channel An- 
derson impurity model (29) to a lattice of localized or- 
bitals, fi, one obtains the so-called periodic Anderson 
model (PAM), 

k(T ia- 
i ia 

Both direct hopping of and direct exchange between / 
electrons are neglected here. In a situation of large U 
and negative e/ local moments on the / sites become 
well defined. Employing a Schrieffer- Wolff transforma- 
tion as in the single-impurity case, the PAM maps onto 
the Kondo lattice model 

^ = E ^kcLcka + J E Si • Si (37) 

k(T i 

where the Kondo coupling J is related to the parameters 
of the Anderson model through (30) . 

We start by discussing qualitative features of the An- 
derson and Kondo lattice models. For small interaction 
U, the periodic Anderson model (36) can be expected 
to describe a Fermi liquid with two bands. As detailed 



T 




FIG. 1 Doniach's phase diagram of the Kondo lattice, as func- 
tion of the ratio of Kondo (T^') and inter-moment exchange 

(/) energies - here T^' is the single-impurity Kondo scale, 
measuring the strength of the Kondo effect. The heavy Fermi 
liquid is formed below the coherence temperature Tcoh- The 
behavior of Tcoh across the phase transition will be discussed 
in Sec. III.I. 

below, a Fermi liquid at lowest temperatures can also sur- 
vive in the laige-U limit and thus in the Kondo-lattice 
case - this requires the local moments to be screened by 
a lattice generalization of the Kondo effect. The result- 
ing Fermi liquid, formed below a coherence temperature 
Tcoh, will have a Fermi volume containing both c elec- 
trons and local moments (dubbed "large Fermi volume"), 
consistent with the Luttinger theorem (Oshikawa, 2000). 
The resistivity will follow the usual quadratic T depen- 
dence (24). It is interesting to discuss the evolution of 
this state with temperature: For T » Tcoh (but still as- 
suming well-formed local moments) the system can be 
described as c fermions with a "small" Fermi volume in- 
teracting weakly with a paramagnetic system of localized 
spins. The resistivity is often rather low in this situation. 
In the crossover region T ~ Tcoh the Fermi surface fluc- 
tuates strongly, giving rise to a very high resistivity of 
the order of the unitarity limit - experimentally, this re- 
sistivity maximum is often used to define Tcoh- 

The screening of the local moments, required for FL be- 
havior in the Kondo lattice, competes with interactions 
between the local moments. Such interactions can be due 
to direct hopping or exchange between the / orbitals, but 
are also generated due to the polarization of the conduc- 
tion electrons. This indirect Ruderman-Kittel-Kasuya- 
Yosida (RKKY) interaction is given in lowest (quadratic) 
order in J by 

-f^RKKY = ^^hj Sj • Sj , (38) 

lij = NoJ^FikpRij) (39) 

where F{x) = {xcosx — sina;)/x* and Rij is the dis- 
tance between lattice sites i. j. In the Kondo-scrcenod 
state lij is expected to be renormalized, in particular 
at long distances, but a reliable determination of lij is 
not available at present. The competition between the 
Kondo coupling and the inter-moment interaction will 



10 



govern the phase diagram of the Kondo lattice (Doniach, 
1977). The most natural competitor of the Fermi liquid 
is a magnetically ordered metal, Fig. 1, but in the pres- 
ence of strong quantum effects and geometric frustration 
spin-glass and spin-liquid states may also occur. 

We continue with an analysis of the FL phase of the 
periodic Anderson model - here non-local interaction ef- 
fects are believed to be unimportant. A viable method 
is perturbation theory in U (Hewson, 1993). Since the 
interaction U is assumed to act only between electrons 
in the / level, there is only one self-energy Jl^{k.,co). Ex- 
panding S near the Fermi energy lo = and near the 
Fermi surface k = kp (assuming an isotropic band struc- 
ture) and using the fact that ImS in this region, 
one may define a quasiparticle component of the single- 
particle Green's function G^f{'k,u!) of weight 



1 



di:{kF,uj)-\ 



du) 



u>=0 



located at energy 



e/,k = Z ef + E{kF, 0) + (k - ) • V^E 



(40) 



(41) 



and hybridizing with the conduction band with renormal- 
ized strength V = Z^I'^V . One finds two quasiparticle 
bands with dispersion 



= ^ {e/,k + ek ± [(e-/,k - ek)^ + 4^^] } . 



(42) 



These are the same bands as in the non-interacting case 

(JJ = 0) except that the bare energies ej and V are re- 
placed by e/,k and V, respectively. Neglecting the self- 
energy's k dependence, which is expected to be weak 
compared to the u; dependence, one finds for the quasi- 
particle effective mass at the Fermi energy 



m 

= 1 + 



y2 



(43) 



The specific-heat coefficient and the spin susceptibility 
are found as 



27r2 - 
7=^iV(0), 



l + Fo« 



where A^(0) is the renormalized total density of states at 
the Fermi level 



7V(0) = A^W(0) + A^)°'(0)/Z 



(45) 



with N^'^j (0) the densities of states of conduction elec- 
trons and / electrons, respectively, in the limit U = 0. 
The factor R = (1 + Fq)^^ , often called the generalized 
Wilson ratio, in x (44) expresses the effect of quasiparti- 
cle interactions in terms of the Landau parameter Fq. 

There arc various approximation schemes available, al- 
lowing to estimate the parameters Z, if in the above 



scmi-phenomenological quasiparticle theory. The sim- 
plest one uses slave-boson mean-field theory, in the limit 
U ^ oo (Newns and Read, 1987). One finds an approx- 
imate mapping to a model of two non-interacting, hy- 
bridizing fermion bands, with a quasiparticle dispersion 
given in (42). Assuming a flat conduction-band density 
of states iVi°^ (e) = 1/ (2D) , and for D > e/ , F, the renor- 
malized / level e/ is found as a solution of 



cf - if = -Tln{if/D) 

TT 



(46) 



which defines a characteristic energy scale, T* = if, i.e., 

the distance of the renormalized f level to the Fermi en- 
ergy (r = '!tNc^\o)V'^ is the bare hybridization width). 
The renormalized /-electron density of states can be ex- 
pressed as 



liv(o)(0)= --— . 



2 r 1 

^TfY* 



(47) 



The specific heat and the spin susceptibility are given by 
(44), with the Landau parameter Fq = 0. A finite value 
of can be obtained from the contribution of fluctua- 
tions about the mean field (Houghton, Read, and Won, 
1988). If the / states are sufficiently far below the Fermi 
level, the characteristic temperature assumes within this 
approximation the same functional dependence as the 
Kondo temperature of the single impurity problem: 



Dexp 



2r 



(48) 



where T* <C -D, implying an exponentially small quasi- 
particle weight factor Z T*/r -C 1. This is the regime 
of heavy fermion metals, with large effective mass ratio 
m*/m > 1. 

A similar mean-field approach can be taken to the 
Kondo lattice model (Burdin, Georges, and Grempel, 
2000) , showing that (at least) two energy scales are rele- 
vant for the Kondo lattice problem. The onset of Kondo 
screening upon lowering T happens around the 
(single-impurity) Kondo temperature, whereas the Fermi 
liquid is only established below Tcoh- Typically, Tcoh < 
, leading to a crossover regime which is wider than in 

(44) 

the single-impurity case (dubbed protracted screening). 

In the weak-coupling limit, both and Tcoh can be 
obtained analytically, with the result: 



-^coh 



Dexp 
Dexp 



1 



1 

iVoJ 



FK{nc) , 
Fcohiric) 



(49) 



where Fk and Fcoh arc functions of filling (and shape) of 
the conduction band, and Nq is the conduction-band den- 
sity of states for J = 0. In this approximation the ratio 
Tcoh/T^^ is a function of the conduction-band properties 
only, but is independent of the Kondo coupling J. 
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Beyond slave bosons, the local correlation physics of 
Anderson and Kondo lattice model has been studied us- 
ing the Dynamical Mean-Field Theory (DMFT). The 
DMFT makes use of the limit of infinite spatial dimen- 
sions (Metzner and VoUhardt, 1989), in which the self- 
energy S(k, w) becomes independent of momentum k. 
This corresponds to a mapping of the Anderson lat- 
tice model to an effective Anderson impurity model with 
energy-dependent hybridization function A(a;) (Georges 
et al., 1996). The self-consistency equation relates the 
local / Green's function of the lattice G'"*^ to that of the 
effective impurity model G^iam 



Gi°'^(^) = j 



de- 



(50) 



= [^-e/-A(z)-EKz)]"i=G™(z). 

The effective Anderson model has been solved by Quan- 
tum Monte Carlo (QMC, see Georges et al. 1996) and Nu- 
merical Renormalization Group (NRG, see Bulla 1999) 
techniques. The most accurate study using NRG by Pr- 
uschke. Bulla, and Jarrell (2000) shows that near half- 
filling, ric = 1, the coherence temperature is actually 
larger than the single-impurity Kondo temperature, their 
ratio depending on the Kondo coupling J. More impor- 
tantly, at lower fillings ric < 0.8, the ratio Tcoh/T^'^'' was 
numerically found to be independent of J, and in the 
small-ric limit was proportional to rig. These findings are 
in qualitative agreement with earlier QMC studies (Jar- 
rell, 1995; Tahvildar-Zadeh, Jarrell, and Freericks, 1997). 

The results for the J dependence of the coherence scale 
of Burdin, Georges, and Grempel (2000) and Pruschke, 

Bulla, and Jarrell (2000), namely Tcoh « T^^', indi- 
cate that the so-called exhaustion scenario of Nozieres 
(1985), predicting Tcoh oc {T^^f/D, is incorrect. As 
detailed by Nozieres (2005), the original exhaustion ar- 
gument (Nozieres, 1985) is too simplistic, e.g., it does 
not correctly account for the flow of the Kondo coupling. 
We also note that the exhaustion scenario has been fal- 
sified experimentally: In CcjjLai-jjPbs Kondo behavior 
of single-impurity type has been observed down to lowest 
temperatures for x up to 80% , with Tk = 3.3 K essen- 
tially independent of x (Lin et al, 1987), whereas the 
exhaustion arguments would predict a suppression of Tk 
for X > 0.1%. 

The approximate treatments of the Anderson and 
Kondo lattice models discussed so far cannot capture 
the competing magnetic ordering tendencies arising from 
non-local inter-moment exchange. On the mean-field 
level, this competition can be included by decoupling a 
non-local interaction with suitable auxiliary fields (Igle- 
sias, Lacroix, and Coqblin, 1997; Kiselev, Kikoin, and 
Oppermann, 2002; Senthil, Sachdev, and Vojta, 2003; 
Senthil, Vojta, and Sachdev, 2004). In addition, exten- 
sions of DMFT have been devised (Si et al., 2001, 2003). 
We will return to these aspects in Sec. III. I. We also note 
that non-local interactions may enhance Tcoh compared 
to due to mutual screening of local moments. Ex- 



perimental indications for this have been seen, e.g., in 
CeCoIns by Nakatsuji et al. (2002), see Sec. IV.A.4. 



3. Multi-channel Kondo effect 

Nozieres and Blandin (1980) realized the importance of 
the proper matching of the degrees of freedom of impurity 
and environment for local non-Fermi-liquid behavior in 
Kondo models. This can be discussed using a generalized 
Kondo model, where an impurity spin of size S is coupled 
to several identical conduction bands (labeled a, a = 
1,...., M, and called "channels"): 



H = 



(51) 



where is the local conduction electron spin of band 
a. The Hamiltonian (51) has a high degree of symmetry: 
besides spin rotation invariance there is separately invari- 
ance against unitary transformations in channel space. 

For the low-energy behavior of (51) with J > 0, three 
cases have to be distinguished: (i) S = M/2, here the im- 
purity spin is exactly screened by M electrons of spin 1/2, 
yielding local FL behavior; (ii) S > M /2 yields an effec- 
tive uncompensated spin of size S — M/2 - this under- 
screened Kondo effect leads to non-analytic corrections to 
local FL behavior; (iii) S < M/2 results in overscreen- 
ing, associated with true local NFL behavior - this will 
be discussed in the following. 

The occurrence of a new type of ground state for 
S < M/2 can be understood as follows (Nozieres and 
Blandin, 1980). Suppose 5 = |, and M = 2; if the ef- 
fective exchange coupling j would scale to infinity as in 
the single-channel case considered in the last section, the 
impurity would bind one electron with opposite spin ori- 
entation in each of the two channels. The net spin projec- 
tion of the overscreened impurity site would be S'^ = —Sz 
if Sz is the momentary impurity spin projection. Now, 
virtual hopping processes of neighboring electrons (with 
spin projection Sz due to the Pauli principle) induce an 
exchange interaction between the effective spin S' and the 
neighboring electron spins, which is again antiferromag- 
netic. According to the initial assumption, this coupling 
scales to infinity; this process repeats itself and thus does 
not converge. One concludes that the fixed point at in- 
finite coupling is not stable, i.e., there must exist a fixed 
point at some finite coupling strength. In the limit of 
large channel number M the new fixed point is accessible 
by "poor man's scaling" (Nozieres and Blandin, 1980). 
One finds that the f3 function for the scale-dependent 
coupling j now has a weak-coupling expansion: 



M , 

— J +■■■ 
2-' 



(52) 



replacing (35). Thus there exists a fixed point (/? = 0) at 
j* = NqJ* = 2/M < 1, within the range of vaUdity of 
the expansion. 
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At the new fixed point the impurity spin is not ex- 
actly screened, which destroys the Fermi liquid found in 
the exactly screened case, and leads to local non-Fermi- 
liquid behavior. Exact treatments based upon the Bethe 
ansatz method (Andrei and Destri, 1984; Schlottmann 
and Sacramento, 1993; Tsvclik and Wiegmann, 1984) 
bosonization and conformal field theory (Affleck and 
Ludwig, 1991a, b) have been used to determine the limit- 
ing low-temperature behavior. 

For the M-channel, spin 5 = 1/2, Kondo model (51) 
the most striking result is a finite impurity entropy S{T = 
0) indicating some kind of fractional degeneracy of the 
ground state (Andrei and Destri, 1984), 



^(r = 0)=in 



2 cos 



(m + 2) 



(53) 



The low-temperature specific heat and spin susceptibil- 
ity in zero magnetic field follow the power-law behavior 
(Andrei and Destri, 1984; Tsvelik and Wiegmann, 1984) 



Ximp(r)~r^-\ m>2 



(54) 



and a logarithmic law for M = 2 (Schlottmann and 
Sacramento, 1993). The leading low-temperature power- 
law correction of the resistivity is found as 



p{T) - po ~ T° 



(55) 



where 



ioT S 



\ and M > 2. 



~ 2+M 

As noted by Nozieres and Blandin (1980), perfect 
channel symmetry is rather unlikely due to anisotropics 
in realistic crystalline fields. The presence of chan- 
nel asymmetry will always cause a flow to the single- 
channel Kondo fixed point. Accidentally symmetry- 
breaking fields may be small: in such a situation, the 
multi-channel fixed point would dominate over a size- 
able energy/temperature range, before branching off to 
a single-channel fixed point at some low energy scale. 

An interesting route to two-channel Kondo behavior is 
the quadrupolar Kondo effect. Here, the impurity has a 
doubly degenerate non-magnetic ground state, and the 
role of the two internal states is taken by orbital degrees 
of freedom. Then, the conduction electron spin direc- 
tions provide two independent screening channels, with 
channel symmetry being protected by spin symmetry. 

To date, experimental realizations of multi-channel or 
quadrupolar Kondo physics in Kondo alloys have been 
elusive. For a comprehensive discussion of experimen- 
tal candidate systems and a presentation of theoretical 
results available see Cox and Zawadowski (1998). 



4. Impurity quantum phase transitions 

Quantum impurity models can show phase transi- 
tions at zero temperature - these transitions are special 
cases of boundary QPT, with the impurity being a zero- 
dimensional boundary of the system. At such a transi- 
tion only the internal degrees of freedom of the impurity 



become critical. Impurity QPT are of current interest 
in diverse fields such as unconventional superconductors, 
quantum dot systems, quantum computing, and they will 
also play a central role in the scenario of so-called local 
quantum criticality in heavy-fermion systems, described 
in Sec. III.I.l - for a review see Vojta (2006a). In the 
following, we briefly mention a few quantum impurity 
models which display QPT. 

In general, Kondo-typc models can feature a QCP be- 
tween phases with quenched and unquenched (or par- 
tially quenched) impurity moment. A well-studied ex- 
ample is the pseudogap Kondo model, where the den- 
sity of states of the conduction electrons vanishes as \uj\^ 
near the Fermi level. Here, no screening occurs for small 
Kondo coupling J, whereas the impurity spin is screened 
for large J. The critical behavior depends on the value of 
the exponent r and on the presence or absence of particle- 
hole symmetry. The pseudogap Kondo model is in par- 
ticular relevant to impurity moments in d-wave super- 
conductors where the exponent r — 1 characterizes the 
density of states of the Bogoliubov quasiparticles. Inter- 
estingly, the behavior of the standard two-channel Kondo 
model, discussed in Sec. II. F. 3, can also be understood 
in terms of an impurity QPT: tuning the ratio of the 
coupling strengths of the two channels drives the system 
from one Fermi-liquid phase to another, with the equal- 
strength two-channel situation point being the QCP. 

A further class of models with impurity QPT are those 
where the impurity is coupled to a bosonic bath - this can 
represent spin, charge, or lattice collective modes of the 
environment. Bosonic impurity models have been first 
introduced for the description of dissipative dynamics in 
quantum systems (Leggett et ai, 1987). The simplest 
realization is the so-called spin-boson model, describing 
a spin or two-level system linearly coupled to a bath of 
harmonic oscillators. This system has two phases: a delo- 
calized phase with weak dissipation, where the spin tun- 
nels between its two possible orientations, and a local- 
ized phase where large dissipation suppresses tunneling 
in the low-energy limit, leading to a doubly degenerate 
ground state with a trapped spin. As above, the univer- 
sality class of the transition depends on the low-energy 
behavior of the bath spectral density. SU(2)-symmetric 
generalizations of the spin-boson model include so-called 
Bose Kondo models where an impurity spin is coupled to 
magnetic fluctuations of the bulk material. 

In the particular context of strongly correlated elec- 
tron systems, which often feature Fermi-liquid quasipar- 
ticles and strong spin fluctuations at the same time, the 
question of the interplay between fermionic and bosonic 
"Kondo" physics arises. This naturally leads to so-called 
Fermi-Bose Kondo models where an impurity spin is cou- 
pled to both a fermionic and a bosonic bath: 



H = 



Ekcl^Cka + JS-so + H^ + 7oS • 

ka 



(56) 



Here, the fermionic bath is represented by the spin- 1/2 
electrons Cfeo-, with local spin density Sq, and 4> is the 
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spin-1 ordcr-paramctcr field of the host magnet. Its dy- 
namics, contained in H^, is most naturaUy described by 
a quantum (f)^ theory; under certain conditions this can 
be replaced by a model of free vector bosons. 

Most interesting is the case of a bosonic bath with zero 
or small gap, corresponding to the vicinity to a magnetic 
QCP in the d-dimensional bulk. The RG analysis shows 
that the two baths compete: For large J fermionic Kondo 
screening wins, resulting in a fully screened spin. Large 
7o can completely suppress Kondo screening, driving the 
system into the intermediate-coupling fixed point of the 
Bose Kondo problem, corresponding to a bosonic fluc- 
tuating phase with universal local-moment correlations 
(Vojta, Buragohain, and Sachdev, 2000). The competi- 
tion is captured by the RG equations (Sengupta, 2000; 
Si et ai, 2001; Smith and Si, 1999) 

a/ \ 3 - d 3 
/3(7) = ^-7-7 , 

m = f-h^. (57) 

The phase diagram for the Bose-Fermi Kondo model 
thus shows a Kondo-screened phase, a bosonic fluctu- 
ating phase, and a continuous quantum phase transition 
in between. The suppression of Kondo screening by bulk 
spin fluctuations is of relevance in materials with strong 
magnetism, like cuprate superconductors, and near mag- 
netic bulk QCP. A model of the form (56) also appears 
within extended DMFT, where a Kondo lattice model is 
mapped onto a Bose-Fermi Kondo model with additional 
self-consistency conditions, for details see Sec. III.I.l. 



G. Non-Fermi-liquid beliavior from disorder 

The influence of disorder on metallic systems has a 
wide range of aspects, ranging from diffusive transport 
to Anderson localization, as reviewed by Lee and Ra- 
makrishnan (1985). As discussed above, the (perfectly 
screened) Kondo effect leads to a local Fermi liquid 
around a magnetic impurity. Usually, in a random di- 
lute solution of impurities, the microscopic parameters 
J and A^'o which determine the coupling between mag- 
netic impurity and conduction electrons, acquire well de- 
fined values yielding a unique Kondo temperature for the 
system under consideration. However, in disordered sys- 
tems, a distribution of Kondo temperatures may arise 
from statistically fluctuating J and/or Nq. This might 
occur in metals near the metal-insulator transition as 
exemplified by heavily doped semiconductors (Lakner et 
ai, 1994; Paalanen and Bhatt, 1991; Paalanen et ai, 
1986; Sarachik, 1995), or by alloys with two different 
non-magnetic constituents leading to different local envi- 
ronments. A similar situation arises in disordered hcavy- 
fermion systems, where the competition of Kondo screen- 
ing and interaction of magnetic moments takes place in 
a disordered environment, often leading to non- Fermi- 
liquid behavior (Stewart, 2001, 2006). 



That disorder can give rise to phases with anomalous 
behavior has been known for some time (McCoy and 
Wu, 1968). The behavior is not confined to a critical 
point in the phase diagram but may persist over an en- 
tire region. These so-called Griffiths-McCoy phases were 
discussed originally in disordered classical systems (Grif- 
fiths, 1969). Much stronger singularities are found in 
quantum models as has been established theoretically in 
recent years, see Sec. III.J.2. 

The subject of NFL behavior driven by disorder has 
been reviewed recently by Miranda and Dobrosavljevic 
(2005) . It is therefore not necessary to give full coverage 
of this important aspect of NFL physics here. We will 
come back to this question in Sec. III.J where the effect 
of disorder on magnetic quantum phase transitions in 
metallic systems will be discussed. 

III. FERMI-LIQUID INSTABILITIES AT QUANTUM 
PHASE TRANSITIONS: THEORY 

A quantum system can undergo a continuous phase 
transition at T = upon variation of some non-thermal 
control parameter (Sachdev, 1999; Sondhi et al, 1997; 
Vojta, 2003a). Near the critical point of such a quantum 
phase transition in the itinerant electron systems of inter- 
est here, the finite-temperature behavior is characterized 
by scaling laws with temperature exponents different, in 
general, from those of Fermi-liquid theory. This may be 
considered as a breakdown of the Fermi-liquid state in- 
duced by quantum fluctuations near the critical point. 
As is the case with usual classical continuous phase tran- 
sitions, the different systems fall into universality classes, 
depending on the symmetry properties of the phase tran- 
sition. In contrast to classical phase transitions, also the 
dynamics affects critical thermodynamic properties, and 
therefore a larger number of different universality classes 
can be expected, e.g., depending on the presence or ab- 
sence of an efficient coupling between the order parame- 
ter and fermionic quasiparticles. 

A. Classical vs. quantum phase transitions 

Phase transitions at T = are dominated by quantum 
effects, in contrast to classical phase transitions at T > 0, 
even though both may occur in the same physical sys- 
tem. Any continuous flnitc-temperaturc phase transition 
is classical in the following sense: continuous phase tran- 
sitions have divergent correlation length and time. The 
order parameter </> (magnetization, staggered magnetiza- 
tion, etc.) fluctuates coherently over increasing distances 
and time scales as one approaches the transition. The 
latter implies that there exists a characteristic frequency, 
cl)0, for order-parameter fluctuations, which tends to zero 
at the transition. The system behaves classically (even if 
quantum effects are important at short length scales) if 
the transition temperature Tc satisfies ksTc 3> hcOff,. This 
argument shows at once that quantum phase transitions, 
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for which Tc = 0, arc quahtativcly different: their critical 
fluctuations require a gi/aniMm-statistical description. 

Let us start the discussion with the partition func- 
tion Z = Trexp(— /3iJ), /3 = l/ksT. For a classical 
continuous phase transition, the partition function may 
be represented in terms of the relevant time-independent 
order-parameter field (/>(r) as a functional integral over 
all configurations of (j) 



?o j D</.(r)exp(- 



(58) 



where Fd is the Landau-Ginzburg- Wilson (LGW) free 
energy functional 



Fci{4>} = j 



d^rJ^ci W)] (59) 
</'(r)(-5o+^o'V')</'(r)+uo/(r) + . 



and Z^i is the partition function of the non-critical de- 
grees of freedom, and eo is a microscopic energy scale. 
The parameter 5o depends on temperature and tunes the 
system through the phase transition. 

For the quantum system wc can represent the partition 
function in a somewhat similar form, the difference be- 
ing that the quantum nature of the order-parameter field 
requires that one keeps track of the order in which the 
(j) operators appear. Realizing that exp{—(3H) is a time 
evolution operator on the imaginary time axis, the parti- 
tion function Z{f3) may be represented as a path integral 
over all configurations ^(r, r) with (/>(r, 0) = 4>{r, hf3): 

Z(/3) =Zoj V4>{v, r) exp(-5{.^}) (60) 

where 5" is the Euclidian action, 

S{<l>} = j^^ dT j rf<*rr{<^(r,r)} (61) 

with the Lagrange density C given by the LGW free en- 
ergy expression amended by the kinetic energy of quan- 
tum fluctuations of the order-parameter field, Ckin-, 



(62) 



The actual form of C^m depends on the dynamics of the 
system. Wc note that these general arguments, however, 
do not exclude the possibility that £ is a highly non-local, 
non-linear object, such that an efficient description of 
the critical behavior in terms of an order-parameter field 
(/)(r, r) is not any more useful. We will defer a discussion 
of these complications, which arise e.g. in ferromagnetic 
metals, to Sec. III.H. 

If the dependence of 0(r, r) on t can be neglected, as 
is the case at a finite-T phase transition, the contribution 
from quantum fluctuations Ckin ^ 0, and the r-integral 
in (61) gives back the factor of /3 in the expression (58) 



for the classical partition flmction. This can be seen more 
clearly in the Fourier representation, where S, e.g., in the 
case of an insulating magnet, takes the form 



1 



koJn 



K'^^+Si (63) 



with ujn = 2'KnkBT the (imaginary) frequency and k the 
wavevcctor of the order-parameter fluctuation. The term 
Si stands for the fourth-order term in the LGW func- 
tion. The non-thermal control parameter 5o acquires a 
T-dcpendcnt rcnormalization from the interactions; the 
transition occurs at 5q = 5c where the renormalized 5 
vanishes. In the following we will denote the distance to 
the quantum critical point (QCP) by 



r = 5o- Sc{T = 0) , 



(64) 



which can be tuned by varying pressure [then r = {p — 
Pc)/Pc], magnetic fleld, or chemical composition. 

Upon approaching a finite-temperature phase transi- 
tion, the energy of the characteristic order-parameter 
fluctuations oj^ (proportional to the renormalized value of 
S) becomes eventually smaller than ksT. Then only the 
LUn = term in (63) contributes to the critical behavior, 
which is thus governed by spatial order-parameter fluctu- 
ations only, Eq. (59). In contrast, at T = temporal fluc- 
tuations arc not negligible, but appear to be intimately 
intertwined with spatial fluctuations. The representation 
(60,61) of the partition function suggests that the system 
behaves like a (d-|- l)-dimensional classical system, which 
is anisotropic since the "gradient energy" in the time di- 
rection may be of different order (e.g. flrst order, see 
below) as in the spatial directions (usually second or- 
der), and the time axis is restricted to the interval 
The anisotropy of the fictitious classical system may be 
characterized by the so-called dynamical exponent z, de- 
fined by the scaling of frequency with wavevector ui k^. 
For the LGW function (63), uj scales as and hence the 
bare value of z = 1 in this case. For metallic magnets we 
will have z > 1. (Interactions may change the bare value 
of z below the upper critical dimension of a field theory.) 

The effective dimensionality for a system near a quan- 
tum phase transition is thus des = d + z. In many cases 
des is equal to or larger than the upper critical dimension, 
d+, of the respective field theory; df = 4 for the mag- 
netic transitions to be discussed below. While a phase 
transition with d^s < d'^ is controlled by an interacting 
fixed point and usually obeys strong hyperscaling proper- 
ties, a transition with doff > rf+ can be described within 
mean-field theory, and hyperscaling is violated due to the 
presence of dangerously irrelevant operators. 



B. Scaling properties near quantum phase transitions 

The functional integral formulation allows one to em- 
ploy the well-established picture of scaling near a continu- 
ous phase transition. As such a transition is approached. 
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both the order-parameter correlation length ^ and cor- 
relation time (i-e., the correlation length along the 
imaginary time axis) diverge: 



(65) 



where r (64) measures the distance to the QCP. The 
correlation-length exponent v of the quantum transition 
is different from the one of a possible finite-temperature 
transition at Tc in the same system (which describes the 
divergence ^ ^ |T — Tc !"'''=')• finite temperature, the 
(d + 2)-dimensional quantum system has a finite length 
in the time direction, Lr = (61). Its properties can 
then be deduced from finite-size scaling (Privman, 1990). 

We are now in the position to discuss the properties of 
a system near a quantum critical point, located at T = 0, 
r = (Fig. 2). The QCP is usually the endpoint of a line 
of continuous finitc-T transitions. (Exceptions are low- 
dimensional systems where order at finite T is prohib- 
ited by the Mermin- Wagner theorem, or systems where 
no order parameter can be defined for T > 0, as is the 
case for metal-insulator transitions or transitions in the 
topology of the Fermi surface.) In general, the bound- 
ary of the ordered phase follows Tc oc {—r)^ where '0 is 
the so-called shift exponent. In the immediate vicinity of 
this boundary there is a region of classical non-Gaussian 
criticality. The disordered phase of the system at finite 
T can be divided into distinct regimes: For low T and 
r > thermal effects arc negligible {Lt 3> ^t, equiva- 
lently T ^ r^'^)^ and the critical singularity is cutoff by 
the deviation of the control parameter r from criticality. 
This regime is dubbed "quantum disordered" and char- 
acterized by well-defined quasiparticle excitations; for a 
magnetic transition in a metallic system this will be the 
usual Fermi-liquid regime. For r < and T > T^, but still 
Lt ^ £,T, we are in the "thermally disordered" regime; 
here the order is destroyed by thermal fluctuations of the 
ordered state (yet quasiparticles are still well defined on 
intermediate scales) . A completely different regime is the 
high-temperature regime above the QCP where 3> ir- 
In this "quantum critical" regime, bounded by crossover 
lines T ~ \r\''^ , the critical singularity is cutoff by the fi- 
nite temperature. The properties are determined by the 
unconventional excitation spectrum of the quantum crit- 
ical ground state, where the quasiparticles of the stable 
phases are replaced by a critical continuum of excita- 
tions. In the quantum critical regime, this continuum is 
thermally excited, resulting in unconventional power-law 
temperature dependencies of physical observables. 

Assuming that the critical behavior is governed by ^ 
and the critical contribution to the free energy den- 
sity, /cr = / — /reg, should follow the homogeneity law 



/cr(r,T) = 6-(''+^)/er(r6i/^T&^) 



(66) 



where h is an arbitrary scale factor. Note that this "naive 
scaling" (equivalent to hyperscaling) is valid only below 

the upper critical dimension, rfoff < , and we will 
comment on deviations later on. Choosing 6 = ^, Eq. (66) 
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FIG. 2 Generic phase diagram in the vicinity of a continu- 
ous quantum phase transition. The horizontal axis represents 
the control parameter r used to tune the system through the 
QPT, the vertical axis is the temperature T. The dashed 
lines indicate the boundaries of the quantum critical region. 
The lower crossover lines are given by T cx \'r\"^; the high- 
temperature crossover to non-universal (lattice) physics oc- 
curs when the correlation length is no longer large to mi- 
croscopic length scales. The solid line marks the finite- 
temperature boundary between the ordered and disordered 
phases. Close to this line, the critical behavior is classical. 



can be cast into the scaling form = ^ ^'^"'"^Vi(^r/i'r), 
or, equivalently, the ansatz: 



/c] 



T 



(j'l/Ci.z)) 



(67) 



where T is measured in units of To, To and po being non- 
universal constants, while 01,2,3(2;) are universal scaling 
functions. 

From Eq. (67) we can immediately deduce the critical 
contribution to the specific heat, C = TdS/dT, as 



Ccr(r = 0,r)(xT'^/" 



(68) 



in the quantum critical regime. If the quantum disor- 
dered regime of Fig. 2 is a Fermi liquid, then Eq. (67) 
yields for its specific-heat coefficient Cci/T{T — > 0) cx 

As is clear from Fig. 2, a quantum critical point can 
be gcnerically approached in two different ways: r ^ 
atT = OorT^Oatr = 0. The power-law behavior of 
physical observables in both cases can often be related. 
Let us discuss this idea by looking at the entropy S. It 
goes to zero at the QCP (exceptions are impurity transi- 
tions discussed in Sec. II. F. 4 and by Vojta, 2006a), but 
its derivatives are singular. The specific heat C will show 
power-law behavior, as does the observable B = dS/dr. 
At a pressure-tuned phase transition, r = {p — Pc)/Pc, B 
measures the thermal expansion. 



1 ^ 
V df 



1 ^ 
V 'dp 



(69) 
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B/C defines the Griineisen parameter T, 

1 (dS/dp)T 



a 



Cp V„,T {dS/dT)p 



(70) 



where Vm = V/N the molar vohime. Taking the ratio 
of the singular parts of B and C one observes that the 
scaling dimensions of T and S cancel, and therefore B/C 
scales as the inverse of the tuning parameter r. Thus, one 
obtains a universal divergence in the low-T limit (Zhu et 
ai, 2003) 



.(T = 0,r) = ^=Gr\r\-' 



rer(T,r = 0) = GtT- 



(71) 
(72) 



With the help of the scaling ansatz (67) the full scaling 
form of r can be determined, for details see Zhu et al. 
(2003). Remarkably, in the T — > limit, even the prefac- 
tor Gr is universal and given by a combination of critical 
exponents. Further we note that F does not diverge at 
a finite-T phase transition, thus a divergence of F is a 
unique signature of a continuous QPT. 

If the control parameter of the QPT is not pressure 
but an external magnetic field H, the quantity B is the 
T-derivative of the magnetization M, and the role of the 
Griineisen ratio is played by 

{dM/dT)H _ 1 {dS/dH)T _ 1 dT 
" CH T {dS/dT)H T dH g'^ ' 

It can be determined directly from the magnetocaloric ef- 
fect by measuring the change of temperature in response 
to an adiabatic {S = const) change of H. 

As the scaling arguments can be invalid above the up- 
per critical dimension, we will quote concrete results for 
the critical points of metallic magnets in Sec. III.D. 

We finally turn to dynamical scaling. Any physical 
quantity depending on r and t (or equivalently k and uj) 
in the critical region close to the QPT (but sufficiently 
far from the associated finite-T transition) should depend 
on space and time only through the scaled variables fc^ 
and w^T! since ^ is the only length scale and is the 
only time scale in that regime. (Note that multiple time 
scales may be present in a multi-component system, see 
the discussion in Sec. III.H.l.) The Fourier components 
of a physical quantity X affected by the transition are 
thus expected to exhibit the following scaling behavior 



(74) 
(75) 



where dx is the scaling dimension of the observable X. 
Exactly at the quantum critical point this reduces to 

X{k,w;r = 0,T=0) = k-*^" F*{k' /lj) (76) 

We again note that all scaling relations are only expected 
to be valid if the critical point satisfies hyperscaling prop- 
erties, which is true below the upper critical dimension 
. Scaling above d'^ in the presence of a dangerously 
irrelevant variable will be briefly discussed in Sec. III.D. 



C. Itinerant fermion systems 

Quantum phase transitions in itinerant electron sys- 
tems were first studied in a pioneering paper by Hertz 
(1976). Hertz pointed out that near a phase transition 
at T = static and dynamic properties are inextricably 
mixed and applied a renormalization group (RG) treat- 
ment to model systems of this type. This work was later 
reconsidered and extended by Millis (1993). 



1. Definition of the Hertz model 

In the context of strongly correlated electron systems 
one is mainly interested in magnetic phase transitions 
in metals. As prototypes we will consider ferromagnetic 
(FM) and antiferromagnetic (AFM) phase transitions. 
We will assume the collective behavior near the tran- 
sition to be characterized by a real A''-component order- 
parameter field (j), representing the magnetization (for 
the FM) or the staggered magnetization (for the AFM). 
A number of simplifications occur in the limit ^ oo, 
although the actual number of components is A^ < 3. The 
effective action may be derived from the Hamiltonian ei- 
ther by introducing the collective field in functional inte- 
gral representation and integrating out the electron de- 
grees of freedom (Hertz, 1976), or by more conventional 
techniques (Moriya, 1985). Assuming that the resulting 
action S{<p} can be expanded in powers of </> with spa- 
tially local coefficients, one arrives at the Hertz model, 



S = S2 + S4 + ... . 

Here the second-order term is given by 



7(fc) 



k,aJn 



-k, — Wrj 



(77) 



(78) 



where the prefactor of 0^ is nothing but the inverse of 
the dynamical spin susceptibility X~^(k, w„). In this case 
the microscopic correlation length is ^ A;^^, where kp 
is a Fermi wavevector, and eo is the microscopic energy 
scale, given by the Fermi energy ep. The momentum 
summation extends up to a (bare) cutoff Aq. 

The dynamic contribution |a;„|/7(fc) accounts for 
damping of the spin fluctuations (?ik,aj„ by particle-hole 
pairs excited across the Fermi level ( "Landau damping" ) . 
Their phase space increases linear with lu. For a ferro- 
magnetic transition (or other transitions with a Q = 
order parameter), ^{k) = vpk ask ^ 0, i.e., the damping 
rate diverges due to the abundance of particle-hole pairs 
with small momentum. This results in a theory with 
(bare) dynamical exponent z = 3. For an antiferromag- 
netic transition j{k) ~ 70, independent of k, yielding 
z = 2. These forms of "f{k) hold if the wavevector of 
the spin mode in either case is well inside the particle- 
hole continuum, i.e., if the ordering wavevector Q con- 
nects points on a ((i — 2)-dimcnsional manifold of points 
on the Fermi surface. For an antiferromagnetic system 
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with a small Fermi volume and a large ordering vector 
Q > 2kp, the particle-hole pairs decouple from the spin 
fluctuations and u> enters quadratically as in (63). The 
crossover from linear to quadratic to dependence has been 
discussed in detail by Sachdev, Chubukov, and Sokol 
(1995); Sachdev (1999). The special situation where an 
antiferromagnetic mode is tied to wavevector 2kF at the 
edge of the particle- hole continuum ("nesting") will be 
considered separately in Sec. III.H.5. 

The fourth-order term 6*4 of the action accounts for 
the self-interaction of spin excitations, 

Si = uoj dT j d''r[cl>{r,Tff (79) 

with Mo denoting the strength of the interaction. 

Let us point out here that the damping term in the 
Hertz theory has been derived under the assumption 
of Fermi-liquid behavior of the electronic quasiparticles. 
This needs to be justified a posteriori and is discussed in 
Sec. III.H. We also note that in the ordered phase, i.e., 
r < 0, T < Tc in Fig. 3, the action (78) does not ap- 
ply: the form of the damping term will be modified due 
to the appearance of a gap in the electronic band struc- 
ture. (This is already clear from the Goldstone theorem 
which requires the mode damping to vanish as A: —» in 
the ordered phase.) Technically, the limit of vanishing 
order parameter, (c6) 0, docs not commute with the 
long-distance expansion, k,u! ^ 0. A discussion of the 
field theory in the ordered state can be found, e.g., in 
Sachdev, Chubukov, and Sokol (1995). 



2. Pressure vs. field tuning 

Frequently, antiferromagnetic critical points are ac- 
cessed by tuning pressure or magnetic field. The presence 
of a magnetic field changes the universality class of the 
system, as its presence breaks time-reversal invariancc 
which leads to a different dynamics of the order param- 
eter. Provided that the system has spin rotation invari- 
ancc perpendicular to the field, a finite uniform magneti- 
zation leads to a precession of the AFM order parameter 
described by an additional term in the action. 



-I- J 



(80) 



where b is the effective exchange field parallel to the 
magnetization. As this term changes the dynamics, it 
affects the quantum critical behavior. This effect is 
most drastic in an insulator (or an itinerant AFM with 
Q > 2fci?) where the dynamics arises from a term of the 
form J J {dr4>Y ill the absence of magnetic fields. There- 
fore the dynamical critical exponent is given by z = 1 
for b = 0. In contrast, Eq. (80) implies z = 2, and the 
QCP has the same critical properties as the superfluid 
quantum phase transition of bosons driven by a change 
of the chemical potential. 



In itinerant magnets, the precession term iujn4'x4'y 
competes with the Landau damping |cc;„|0^, which both 
have the same scaling dimension. Technically, the term 
in Eq. (80) is an exactly marginal perturbation, and crit- 
ical properties depend on the exact ratio of precession 
and Landau damping, see Fischer and Rosch (2005) for 
an extensive discussion. (Experimentally, many systems 
have a strong Ising anisotropy due to spin-orbit interac- 
tions, rendering the precession term irrelevant.) 

Besides suppressing antiferromagnetism, a uniform 
field can also induce large non-analytic changes in the 
uniform magnetization. These phenomena, usually oc- 
curring in almost ferromagnetic systems, are referred to 
as metamagnetism. Frequently, metamagnetic transi- 
tions, being not associated with a symmetry change of 
the system, are of first order at low T and feature a finite- 
temperature critical endpoint. However, by utilizing ad- 
ditional tuning parameters, the critical endpoint may be 
suppressed down to T = 0, resulting in a quantum criti- 
cal endpoint. Such a scenario has been proposed for the 
bilayer ruthenate Sr3Ru207 (Grigera et al., 2001). The 
theoretical description starts out from the Hertz model 
(77) for a ferromagnetic order parameter, supplemented 
by a sixth-order term in (j). The resiilting phase transition 
is qualitatively similar to a Hertz-type Ising transition in 
a system with 2; = 3; for a detailed discussion we refer 
the reader to Millis et al. (2002). 



3. Scaling equations 

The model defined by the action (77,78) has been stud- 
ied near its critical point (Hertz, 1976; Millis, 1993) us- 
ing the perturbative RG. To define a RG transformation 
one investigates how a change of the cutoff (either in 
momentum or frequency space or both) and their subse- 
quent rescaling can be absorbed in a redefinition of the 
coupling constants, using for example the perturbative 
expressions for the free energy (Millis, 1993). Following 
Millis (1993), we use a scheme where simultaneously the 
cutoffs in momentum space, Afe — > Afe/6, and frequency 
space, A^^/b^, are reduced. The change of 5, u, the 

dimensionless temperature T, and the dimensionsless free 
energy density = F£^q/{ToV) under infinitesimal RG 
transformations are given by (Millis, 1993): 



dT 

dinb 
d5 

dlnb 
du 

dlnb 
dT 

dlnb 



= zT, 



= 26 + 4{N + 2)uf2{T,d)^ 



(81a) 
(81b) 



{A-d-z)u- 4{N + 8)u^f4{T, S) , (81c) 



{d + z):F-^MT,5) 



(Bid) 



with the initial conditions T = T, S = Sq, u = wq, = 0. 
While these equations differ in prefactors from the ones 

of Millis (1993), the derivation of the fi terms can be 
found there. Briefly, the fi are explicit functions of T, 
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FIG. 3 Phase diagram of the Hertz model. I: Fermi-liquid 

regime: II/III: Quantum critical regime; IV: Non-Gaussian 
classical critical regime; V: Magnetically ordered phase. The 
regimes II and III are distinguished in the behavior of the 
correlation length ^, see (97). The quantum critical regime 
also extends into the ordered phase, with singular behavior for 
T > T* similar to regime II. Note that transport properties 
may show more complicated crossovers, see Sec. III.F. 



S, and the bare cutoff Aq (which we will set to unity in 

the following). The dependence of the number of order- 
parameter components N can be absorbed by defining 
/o = Nfo/2, /2 = 4(iV + 2)/2, /4 = 4(iV + 8)/4. At low 
T close to the critical point, 5 = 0, we have 



fo{T) 
UT) 



/o(0) + -iiKa 



in' 



+ ... 



(82) 



with Kd = {2'^-^w'^/^T{d/2))-'^ = l/27r'^-\ d = 2,3. At 
high T, i.e., in the quantum critical regime, 



foiT) = DT, h{T) = CT 



(83) 



where D, C are constants of order unity, J4, does not have 
any critical dependence on T or 5 and may be replaced 
by a (positive) constant. 

The RG equations (81) have a Gaussian fixed point at 
T = u = 5 = Q, which is unstable with respect to the tun- 
ing parameter 5. Of particular interest is the differential 
equation (81c) for the quartic coupling u: if d-l- z > 4, u 
scales to zero, i.e., the upper critical dimension is given by 
(i+ = 4— z rather than = 4 as for classical critical phe- 
nomena. Thus in many cases of interest one can expect 
the critical behavior to be that of the Gaussian model. 
For example, in the cases considered above, z = 2 (anti- 
ferromagnetic metal) or z = 3 (ferromagnetic metal), and 
in dimensions d>2 the Gaussian model applies (the case 
d = z = 2 is marginal and needs special consideration). 



4. Solution for d -I- 2; > 4 

We will now discuss the results obtained by solving 
Eqs. (81). Interpreting 6 as a flow variable, straight- 



forward integration, employing the Gaussian approxima- 
tion, yields the scale-dependent quantities 



T(6) 
u{b) 

5{h) 



uob' 
b^ 



i-d-z 



i: 



So + uo / dbib^ 



,l-{d+z) 



(84a) 
(84b) 

f2m)\. (84c) 



Under the RG process the scale b increases and the sys- 
tem moves away from criticality, i.e., 5{b) increases, until 
S{b) = 1 is reached at 6 = 6o- The regimes I and II/III in 
Fig. 3 are distinguished according to whether the renor- 
malized temperature is smaller or larger than the cutoff 
in energy, T(6o) ^ 1. 

The condition T(6o) ^ 1 defines the quantum disor- 
dered (or Fermi-liquid) regime (Hertz, 1976). To obtain 
the limits of this regime, one may put T = in the 
equation for S{b), to determine 60 from S{bo) = 1. From 
T (60) < 1 one finds with the aid of (84) 



(85) 



where v = 1/2 is the Gaussian correlation- length expo- 
nent (for d + z > 4). and r (64) is the distance to the 
QCP as above. In lowest order perturbation theory we 
have r = Sr with the renormalized control parameter 



Sr = So + 



z + d-2 



(86) 



and (5,. > is required here. The relation T* ~ \r\^^ 
marks the crossover line from the quantum disordered 
regime to tlic quantum critical regime. (We note that 
a line T* ^ \r\^'^ also exists inside the ordered phase, 
Fig. 3, with singular contributions to thermodynamics 
and transport above T*.) 

In the opposite case, when T(6o) ^ 1 at S{bo) ^ 1, it 
is convenient to perform the scaling successively in the 
regimes T(6) <C 1 and T(6) » 1. Starting at small scales 
the results derived above may be used to determine 5(b) 
and u{b) at the scale 61 where T(6i) ~ 1 and hence 61 = 
T-V^. The result is 



^1 



Ml 



w(5i) = woT(''+^-^>/^ 



(87) 



These values provide the initial conditions for the scaling 
in the regime where T(6) > 1 and /2(T) ~ CT. With 
the new scaling variable v = uT one can decouple the 
flow of T and v 



dS 
dMb 



= 25 + Cv, 



dv 
dinb 



= (4 - d)v. (88) 



Integration of these equations starting from the initial 
conditions vi = uiT{bi) = ui,di yields for dimensions 
d> 2 (assuming d + z > 4), 



v{b) 
Sib) 



viib/br)^-'' 
Si{b/b,)' + 



(89) 



(6/6i)2 - (6/&1) 



4-d 



Cvi 

id-2 
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In > 2 dimensions, one finds therefore 

C 



5{b) = 



Sr 



B 



d-2 



where the so-called shift exponent 

z 



d + z-2 



(90) 



(91) 



describes the position of the finite-temperature phase 
transition (see below). For d = 2 (and z > 2) logarithmic 
terms appear in the solution for S{b): 



S{b) = (6/6i)2 di + en{b/bi)Cvi) 



(92) 



For dimensions 2 < d < 4, both v{b) and S{h) increase 
as the scaling proceeds and b is growing. If (5(6o) ^ 1 
is reached when w(6o) ^ 1, the Gaussian approximation 
used here is sufficient. If, however, v{b) becomes of or- 
der unity while 6{b) is still small, the scaling leaves the 
weak-coupling regime and crosses over to a new regime 
characterized by non-Gaussian behavior. The condition 
for Gaussian behavior, the so-called Ginzburg criterion, 
is thus v{bo) <^ 1, or more explicitly, in = 3 dimensions. 



vibo) 



uT 



[(5,. + (B + C)wTi/V']iA 



« 1 



(93) 



correcting the exponent in the numerator of Millis (1993). 
The condition is violated within a narrow strip around 
the temperature for which the denominator in (93) van- 
ishes. One can therefore locate the transition tempera- 
ture Tc((5, 6) from the condition v{bo) = 0{1) which for 
all z > 1 leads to 



(B + C)u 



(94) 



for r = 6r < 0, see Fig. 3). In (94) we have reinstated the 
microscopic energy scale egj to indicate the real temper- 
ature scale. From Eq. (91) we see that in d = 3 the tran- 
sition temperature varies as Tc oc (— r)^/^ and (— r)^/^ 
for antiferromagncts and ferromagnets, respectively; the 
d — 2 result is in Eq. (95) below. 

The width of the critical region around Tc, where the 
Gaussian approximation is expected to fail, can be ob- 
tained from w(6o) = 1- For d -|- z > 4, as assumed, the 
width of the critical region is seen to shrink to zero as 
Tc 0, thus vahdating (94). 

In two dimensions, the presence or absence of the phase 
transition depends on the symmetry of the order param- 
eter. The Ginzburg criterion can be used to locate the 
boundary of the regime where both the correlation length 
is large and the coupling strong. For Ising (XY) sym- 
metry this will give an estimate for the phase transi- 
tion (Kosterlitz-Thouless) temperature Tc- Defining b* 
by the Ginzburg condition v{b*) = 0{1) we obtain from 
6{b*) = 1 for Tc in d = 2, z > 2 



-5r 



u{l + B + Cfe[l/(Virc'""''^'")]} 



(95) 



Up to a logarithmic correction, the transition tempera- 
ture in c? 2 is linear in (— r). 

The correlation length ^ of the order-parameter fluc- 
tuations is given by S{bo) ~ (6o/C)^- I^i the quantum 
disordered regime, T <C Ir]^" , the correlation length is 
determined by the distance to criticality. 



(96) 



with IV = 1/2 and r = 5r (86). In the quantum critical 
regime, T |r|^'', one finds for d = 3: 



6r + {B + C)uT 



1/^ 



(97) 



where the first (second) term dominates in regime II (III), 
see Fig. 3. Note that all prcfactors in Eqs. (90 97) which 
involve both B and C are not exact, as they neglect cor- 
rections of order unity which arise from the crossover 

regime, where T{b) k, 1. Including this intermediate 
regime one obtains for c? > 2 (Garst, 2003) 

C^ = 5r + A{N + 2 f'''^ ^ T"^- (98) 

zcos(%^7r) 

For d = 2, z = 3, Eq. (97) is replaced by 



^-■^ = S^+[B + C£ny ] uT 



(99) 



with bi = T and bo determined from 5{bo) = 1 and 
(92). 



5. Solution for d + z = 4 

The marginal case d = z = 2 requires special con- 
sideration, as in this case logarithmic corrections to the 
Gaussian behavior appear. The scaling of the interaction 
u is now governed by the second term in Eq. (81c). In- 
tegrating (81c) taking only the leading constant term of 
/4 into account gives 



uib) = 



Uo 



1 + uofiinb 



(100) 



For large b (i.e. large bo, at the endpoint of the scaling 
process), the scaled interaction is seen to be independent 
of Uo, and decreases ~ l/£nb. 

Next, integrating the scaling equation (81b) for S{b) 
using (84a) and the scaled interaction (100), one finds 



6{b) = b'^Sr + f3{b, T)T'^ 



(101) 



where the renormalized control parameter 5r is given by 

„-2fa6i 



5r = 5o+ M0/2 



no 

(0) / 
Jo 



dinbi 



1 + uofitnb\ 



(102) 



and /3(6, T) can be expressed in a power series in T^. 

The quantum disordered regime is again defined by 
the condition T(&o) <C 1 at 5{bo) ~ 1. Observing that 
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5(&o) = b^Sr in that regime, the boundary of the quantum 
critical regime is given by T ^ (5^. In the quantum critical 
regime, T ^ Sr, the scaling up to ~ 6t = T~^/'^, where 
Tipr) ~ 1 yields (with G = const) 

6, ^ 6{bT) = I + A (103) 



T 



Ui 



uibr) 



< 1. 



(104) 



instead of Eqs. (87) found in the case d + z > A. 

The correlation length in the quantum disordered 
regime is still given by Eq. (96), while in the quantum 
critical regime we have 



5r I in— 

Or 



N + 2 
' JV + 8 



Jnen{c/T) 



(105) 



en{c/T) ' 

where the first (second) term dominates in regime II (III). 

It should be emphasized that the expansion in regime 
III is not very revealing as it is only valid if £n£n{c/T) ^ 
1. Obviously, and as pointed out by Sachdev and Dunkel 
(2006), this condition will never be satisfied in practice. 
Thus, the d = z = 2 quantum critical theory is in general 
not in a weak-coupling regime at any T > 0. Instead, a 
strongly coupled effective classical model emerges which 
can be used to determine the critical dynamics, for details 
see Sachdev and Dunkel (2006). However, in the con- 
text of metallic antiferromagnets of interest here, even 
more serious complications arise, which we discuss in 
Sec. III.H.2. 



D. Thermodynamic quantities 

As shown above, in the cases where d + z> A the sys- 
tem scales to the non-interacting (Gaussian) fixed point. 
Although the quartic coupling u is formally irrelevant, 
it still affects a number of physical quantities, like the 
location of the critical line or the order-parameter sus- 
ceptibility at r = 0, T > 0. To describe this within a 
scaling approach, as in Sec. III.B, one has to include ex- 
plicitly the quartic coupling u in the scaling ansatz for 
the free energy 



^rp(d+z-A)/z 



(106) 



which replaces Eq. (67). While u is irrelevant, it has to be 
kept as 04 can become a singular function of u. In such 
a case the "naive" scaling relations derived from Eq. (67) 
are modified, and u is called "dangerously irrelevant" . 

The free energy in the theory of Hertz is obtained by 
integrating the RG equation (81d) for the free-energy 
density up to the scale where 5 (bo) = 1 and adding the 
Gaussian free energy 



Fg(6) = - j d'^kj ^cothj^ 



2T(6) 



X arctan 



7(fc)[(5(&) -hfc2 



(107) 



d = 2 
z = 2 



2 = 2 



d = 2 
z = 3 



d = 3 
z = 3 



C'ci 



r..„ = (rini) ' -(2r)-i (2r)-i (rin\) ' 

TABLE I Results for thermal expansion a (69), specific heat 
C, and Griineisen ratio T (70), at a QCP in the quantum 

disordered (Fermi-liquid) regime T <C r'^/^ The left columns 
{z = 2) are for a metallic AFM, the right ones (z = 3) for 
a metallic FM. Non-universal prefactors of Ocr and Ca are 
not shown. The prefactors of Fcr and Fff^cr are (up to the 
logarithmic correction for d = z) universal; these quantities 
are given by Fcr = {dr/dp)Tr,ci/Vm for prcssuro-tunod QCP 
with r — {p — Pc)Ipc, and Vh,ct: ~ {dr/dH)Tr,ci for field- 
tuned QCP with r ^ [H - He) '/He Note that for d = 3, 2 = 
2 the specific heat is dominated by a non-critical fermionic 
contribution C ~ T (Zhu et al., 2003). 



at this scale. The Gaussian contribution to the free en- 
ergy at the scale bo is then Fq = b^^'^'^^^ Fcibo). The 
scale-dependent correction to the free energy, arising in 
the scaling process, is obtained by integrating the scaling 
equation (81d) for J^{b) up to 60: 



Fsc 



ptnho 

Jo 



dinbib^^'^+'^foiTihj) 



(108) 



Depending on the values of d and z, different terms will 
dominate the critical contribution to the free energy; for 
details see Zhu et al. (2003). 

In Tables I and II we present the results for the critical 
contributions to the specific heat C, the thermal expan- 
sion a, Eq. (69), and the resulting Griineisen ratio F, 
Eq. (70), derived from the free energy F = Fq + Fsc, in 
both the quantum disordered and the quantum critical 
regimes. Up to logarithmic corrections the results agree 
with those derived from scaling in Sec. III.B. Note that 
for = ^; the prefactor in (71) vanishes. The 1/r depen- 
dence of ttcr for d = z arises from a T'^lnl/r correction 
to i^cr not captured by scaling. For the quantum critical 



d=2 d=3 
z=2 z=2 



d = 2 
z = 3 



d = 3 
z = 3 



acr ~ Mn^ T^/^ in 



nl/3 



Ccr~ Tin^ -r3/2 T^/s Tin^ 



TABLE II Results for LOW QCP in the quantum critical 
regime T > \r\''^'^ (cf. Table I). 
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regime uy d = \/v = 2 the thermal expansion is loga- 
rithmic. The argument of the logarithm is a power of T 
for d + z > A and is itself logarithmically dependent on 
T for d + 2; = 4; these features reflect the dangerously 
irrelevant or marginal nature of the quartic coupling uq. 

In addition to the critical contributions, the measured 
quantities also contain non-critical background compo- 
nents. We list here the full results for the purpose of com- 
parisons with experiments in heavy-fermion compounds 
undergoing an antiferromagnetic transition [z = 2). Con- 
sider first d = 3. At the QCP (r = 0) 

a = aiT^/^ + azT, (109) 

where the 02 term comes from the (fermionic) back- 
ground contribution. However, approaching the QCP in 
the Fermi- liquid regime gives: 

a = [ai/r^''^ + a2)T . (110) 

For d = 2 and 2 = 2, we have at the QCP (r = 0) 

a = aiin[Mn^] + a2T , (111) 

and in the Fermi- liquid regime approaching the QCP: 

a = {ai/r + a2)T. (112) 

In two dimensions, the thermal expansion coefficient at 
r = diverges in the T = limit, Eq. (Ill), in sharp 
contrast to the textbook statement that a{T ^ 0) = 0. 

Finally, we turn to other thermodynamic quantities. 
The static susceptibility for a ferromagnet is x ^ 
However, x is not critical in the case of an antiferromag- 
net and therefore the calculation is delicate. loffe and 
Millis (1995) performed a careful calculation of the var- 
ious contributions and concluded that in d = 2 (with 
z = 2), the leading low-T dependence of the susceptibil- 
ity is given by x w xo ^ DT with non-universal constants 
Xo and D (probably positive for typical band structures). 
In analogy, we expect x ~ Xo — D'T^^"^ in = 3. 

For a QPT driven by a magnetic field B (see 
Sec. III.C.2), however, the susceptibihty, x = dM/dB = 
—d^F/dB^, is more singular since the uniform magnetic 
field is a relevant perturbation acting as the control pa- 
rameter, r oc B — Be- Below the upper critical di- 
mension one obtains from scaling x T(d+z-2/iy)/z 
the quantum critical regime. For an itinerant magnet, 
d = 3,z = 2, Fischer and Rosch (2005) found in contrast 
that the critical contribution Xc of x is a singular func- 
tion of the irrelevant coupling u, Xc ~ i{T)T ~ T^l"^ j^Ju. 
It therefore strongly violates scaling. 

E. Self-consistent spin-fluctuation theories 

Our discussion in Sees. III.C and III.D was based on 

a renormalization group analysis following the work of 
Hertz (1976). However, many of the main results, e.g.. 



for the phase diagram and susceptibility of nearly FM 
(Moriya and Kawabata, 1973) or AFM metals (Hasegawa 
and Moriya, 1974) in three dimensions, have been ob- 
tained earlier by Moriya and coworkers. 

Their so-called SCR (self-consistcntly rcnormalizcd) 
formahsm, described in detail by Moriya (1985), is a 
self-consistent one-loop approximation for the scattering 
of spin fluctuations (79). Within the RG approach one 
can understand the enormous success of SCR theory and 
where it fails. Technically, one obtains the SCR result 
from the RG equations (81) by neglecting the renor- 
malizations in Eq. (81c) and by replacing the running (5(&) 
on the right-hand side of Eq. (81c) self-consistently by 

lim(,_>oo 5{b)/b'^. The SCR formalism works above the 
upper critical dimension, especially for d = 3 and z = 2 
or 3; these approximations reproduce correctly the lead- 
ing behavior of the; relevant physical quantities including 
prefactors [e.g. of the correlation length, Eq. (98), con- 
trary to claims of Millis (1993)]. The SCR formalism fails 
upon approaching the classical transition as it is blind for 
the Ginzburg criterion. A recent extension of the SCR 
method (Moriya, 2006) has been used to obtain the cor- 
rect behavior at the upper critical dimension, + z = 4. 

The importance of the SCR theory lies in its impres- 
sive success in quantitatively describing a large num- 
ber of different /- and d-electron systems (Lonzarich 
and Taillefer, 1985; Moriya, 1985; Moriya and Takimoto, 
1995), predicting, e.g., transition temperatures by us- 
ing parameters obtained from neutron scattering. It 
is therefore often the method of choice to fit experi- 
ments, in order to find out whether they are described 
by weakly interacting spin fluctuations. Here it turned 
out to be useful to describe the interaction of the spin 
fluctuations not by a (f)^ term as above, but instead - 
in the spirit of a non-linear sigma model - by a con- 
straint set by the sum rule for the dynamic susceptibility 
{SD = (27r)-iEq/'^^coth(a;/2T)ImXaa(q,c^) where 
S is the spin within the unit cell (Kambe et al., 1996; 
Lonzarich and Taillefer, 1985; Moriya, 1985; Moriya and 
Takimoto, 1995). Together with a parameterization of x 
as in Eq. (78), this constraint can be used to determine 
the T dependence of the renormalized mass 6. 

It would be worthwhile to develop the SCR theory fur- 
ther for the ordered side of the phase diagram, where in 
our opinion the formalism (Lonzarich and Taillefer, 1985; 
Moriya, 1985) is not as well justified as on the paramag- 
netic side. 



F. Transport properties 

For quantum critical systems above the upper critical 
dimension, the low-T thermodynamics can be calculated 
reliably using standard methods like RG (see Sec. III.D). 
Much less is known about transport properties. Earlier 
predictions of the resistivity close to an antiferromagnetic 
QCP (see Moriya 1985 and references therein), for exam- 
ple, are not valid for perfectly clean (Hlubina and Rice, 
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1995) or weakly disordered samples (Rosch, 1999, 2000). 
The difficulty arises because scattering from AFM spin 
waves is extremely anisotropic and affects only a small 
fraction of the Fermi surface. Therefore the transport 
properties depend both qualitatively and quantitatively 
on how other scattering mechanisms re-distribute quasi- 
particles and scatter them into these small regions. 

In d = 3 and for small static electric fields, it is possible 
to treat transport within a simple Boltzmann approach, 
as the QCP is above its upper critical dimension, spin- 



J 



spin interactions are irrelevant in the RG sense, and be- 
cause the concept of quasiparticles is still (marginally) 
valid in three dimensions. In d = 2, a quasiparticle de- 
scription is not possible at the QCP, see Sec. III.H.3, and 
a transport theory is more complicated (Kontani, Kanki, 
and Ueda, 1999). In the linear- response regime in d = 3, 
the quasiparticle distribution = — <l>k(c?/k/t^ek) is 
linearized around the Fermi distribution and the col- 
lision term reads (Hlubina and Rice, 1995) 



dt 



coll 



E 



($k - $k' 



2o^ 

fi'?mp'^(ek - ek') + -=^n°^_^^.Imxk-k'(ek-ek' 



(113) 



r 



Here g^^^ and g| arc transition rates for impurity scat- 
tering and inelastic scattering from spin fluctuations, de- 
scribed by the susceptibility Xq(w), and /° and are the 
Fermi and Bosc functions, respectively. In the derivation 
of (113) one has assumed that the spin fluctuations stay 
in equilibrium, i.e., drag effects are neglected. This ap- 
proximation implicitly assumes the presence of sufficient 
momentum relaxation, e.g., by strong Umklapp scatter- 
ing. While this approximation gives probably at least 
qualitatively correct results in the case of a nearly AFM 
metal with a large Fermi volume, it is less clear whether 
it is valid for a FM where small momentum scattering 
dominates. The linearized Boltzmann equations in the 
presence of electric and magnetic fields, E and B, can be 
written in the following form: 



VkE + (vk X B)ak$k = 

2 r°° 



"kk' 



■($k-$k'),(114) 



i;k'|(2^)3 

^kk' = gLp + l'^^''°o.K + l]lmxk-k'(a;) 



where an integration over energy ek' has been performed 
and all k-vectors, Fermi velocities Vk and integrations 
are restricted to the Fermi surface, using JJ dk/\vu\ = 
J d'^k(5(ek — Currents are calculated from = 

JJvlM/[\vu\{27rn 

If quasiparticles are scattered almost equally strongly 
all over the Fermi surface, the k dependence of $k is 
for B = given by $k oc Evk. Only under these as- 
sumptions, one recovers results for the resistivity at the 
QCP which have been derived many years ago by Mathon 
(1968); Moriya (1985); Ueda (1977). 



dk 

\vu\ 



(vkA)' 



(115) 



where n is a unit vector parallel to electric field and cur- 
rent. An almost equivalent form to (115) (including fac- 



tors of e^, h) can be foimd in Moriya (1985). As (115) 
is valid only for approximately uniform scattering, it can 
be used in the case of a FM QCP or if the scattering 
is dominated by (short-range) impurities. In the quan- 
tum critical region, T > T* , of a QCP in d = 3 one 
obtains (Mathon, 1968; Moriya, 1985; Moriya and Taki- 
moto, 1995; Ueda, 1977) 

A ~ / -^^^^ ^™ disorder-dominated, Ap < /'O/ii^.x 

P ^ 1 j.5/3 y ) 



Eq. (115) fails completely in the case of an AFM QCP 
in a weakly disordered metal where quasiparticles scat- 
ter strongly close to lines on the Fermi surface with 
Ek = Ek+Q, the so-called "hot lines". This has been 
first realized by Hlubina and Rice (1995), who argued 
that in a clean metal close to an AFM QCP, the resis- 
tivity is dominated by quasiparticles from regions of the 
Fermi surface far away from the hot lines, where scatter- 
ing rates arc proportional to T^. Accordingly, p cx is 
expected in an ultra-clean metal close to the AFM QCP. 
This effect can be understood in a simple relaxation- 
time approximation, where the resistivity is calculated 
from a Fermi-surface average of the k-dependent scatter- 
ing time, p oc 1/{t\s)fs- Clearly, the longest scattering 
times dominate the resistivity and short-circuit contri- 
butions from the hot lines where r is small. This effect 
is missed in (115), where l/(Tk)FS is effectively replaced 
by (l/rk)F5- However, the predicted dependence of 
p is unobservable in real systems, as tiny amounts of dis- 
order change the picture qualitatively - this has been 
pointed out by Rosch (1999, 2000). Rosch (2000) solved 
the Boltzmann equation (114) numerically and analyt- 
ically (in a certain scaling limit). For sufficiently small 
magnetic fields, weak disorder, and close to an AFM QCP 
in the paramagnetic phase in > 2, the temperature- 
dependent part of the resistivity, Ap = p{T) — po, obeys 
a scaling form as function of B, residual resistivity po, 
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and distance from the QCP, here denoted as r > 0: 



y(d-l)/2 j,{d-l)/2 



B 



PO 



po ' poVt, 



(117) 



Most interesting is the case B = Q,r = where the scal- 
ing function crosses over from /(a, 0, 0) ~ const for low- 
est temperature, i.e. a ^ 0, to /(a, 0,0) ~ a-(^-<^)/(5-<^) 
for temperatures above a crossover scale proportional to 
p2/(d-i)^ Various crossovers in d = 3 are discussed below 
in Eq. (118) for B = 0. The rather complicated magnetic 
field dependence can be foimd in Rosch (2000). 

Rosch (2000) has argued that the magneto-resistivity 
and the sensitivity to disorder can be used to decide 
whether an observed NFL behavior arises only from some 
"hot lines" on the Fermi surface typical for a spin-density 
wave transition or from a breakdown of the Fermi liquid 
on the full Fermi surface which is expected in some other 
scenarios (see Sec. III.H). Sensitivity to weak disorder 
and large non-linear effects in the magneto-resistivity are 
characteristic for the Hertz scenario. 

For B = and d = 3, one obtains for clean samples, 
characterized by a large residual resistivity ratio [RRR, 
defined as p{T = 300K)/po], 



Ap(T) 
P{To) 



ty/x , max[a;, ^/rx\ <t < \fx 

t^/^ ,r<t<x (118) 

t"^ / y/r , t < min[r, y/xr] 



replacing (115). The dimensionless temperature, disor- 
der and distance from the critical point are defined by 
t ~ T/To, X ^ po/p{To) ^ 1/RRR, r ^ l/ikpO^ « P'Pc 
(note ziy = 1 here). Tq is a characteristic temperature 
scale, e.g., an efii'ective Kondo temperature of a heavy- 
fermion system. In a sufficiently clean sample with a 
RRR of the order of 100, the resistivity at the QCP is 
linear in temperature over a large regime To /RRR 
Tq/VRRR due to the interplay of weak impurity scatter- 
ing all over the Fermi surface and strong inelastic scat- 
tering close to the hot lines. Even at some distance from 
the QCP, one has to go to rather low temperatures to 
recover Fermi- liquid behavior. In a very clean sample 
with a large RRR, the crossover temperature t ^ \frx 
can be considerably lower than the characteristic scale 
(t ~ r) below which Fermi-liquid behavior is recovered in 
thermodynamic measurements. 

The picture developed above, fits well to some exper- 
imental trends. In a large number of systems, the resis- 
tivity at the AFM QCP seems to rise with T^-^, e.g. in 
CeCu2Si2 (Gegenwart tt ai, 1998), CeNi2Ga2 (Hauser et 
at, 1998), CeCug-xAg^ (Heuser et at, 1998a) CcNi2Ge2 
(Grosche et ai, 2000), CePd2Si2, or Celns (Julian et ai, 
1996). All of these systems are "dirty" in the sense that 
A/9 <C Pa in the temperature range where the above ex- 
ponents have been fitted. In a few systems, the Lonzarich 
group (Grosche et ai, 2000; Julian et ai, 1996; Mathur 
et ai, 1998) succeeded in preparing high-quality samples 
with RRR values of the order of 100 and the resistivity 



seems to rise with exponents smaller than 1.5. In the 
cleanest samples the resistivity is almost linear in T over 
a substantial temperature range. Especially, the sensi- 
tivity to weak disorder might be interpreted as a signa- 
ture that the systems can be described within the theory 
sketched above. Also, resistivity measurements experi- 
ments in U2Pt2ln and U3Ni3Sn4 (Estrela et al., 2001a, b) 
have been fitted to the theory (118). For further compar- 
ison of theory and experiment and a detailed calculation 
of magneto-transport see Rosch (2000). A number of 
systems do not fit into the scenario above described, an 
example being CeCug-ajAua; where Ap is linear in T in 
a regime with Ap ^ po- Remarkably, one obtains a lin- 
ear T dependence from (115) when one assumes that 3d 
electrons scatter from 2d spin fluctuations (Rosch et al., 
1997) - see Sec. IV. A. 1 for a more detailed discussion. 



G. Approach to the QCP from the Fermi-liquid regime 

In the previous sections we used a description of QCP 
in terms of a bosonic order-parameter field (coupled to 
weakly interacting fermions) and identified the bosonic 
fluctuations as the origin of the singular behavior of ther- 
modynamic and transport quantities. However, away 
from the QCP the system is a Fermi liquid in the low- 
T limit and therefore the approach to QCP for T ^ 
can alternatively be described by Fermi-liquid theory as 
mentioned in Sec. II. C. Within this complementary (but 
equivalent) language, the singularities in thermodynam- 
ics, for example, are not associated with bosonic fluctu- 
ations, but arise from the mass renormalization of the 
quasiparticles. One should keep in mind that Fermi liq- 
uid theory is devised to account rigorously for any low- 
energy excitations, but cannot be used to calculate, e.g., 
short range properties such as the momentum depen- 
dence of susceptibilities. 

We investigate the Fermi-liquid description of a QCP 
in the case of a Pomeranchuk instability (c.f. Sec. II. C). 
For such a uniform instability (at wavevector Q = 0) the 
momentum dependence of the quasiparticle properties is 
very weak, simplifying the description considerably. 

We note that "hidden" order caused by a Pomer- 
anchuk instability has attracted considerable interest re- 
cently, e.g., to explain the enigmatic ordering transition 
in URu2Si2 (Varma and Zhu, 2006), and in the context of 
the cuprate superconductors. Several calculations for the 
Hubbard model (Grote, Kording, and Wegner, 2002; Hal- 
both and Metzner, 2000; Neumayr and Metzner, 2003), 
and the t— J model (Yamase and Kohno, 2000) indicate 
a Pomeranc;lmk instability in the spin-symmetric rf-wave 
channel. Pomeranchuk phases in isotropic Fermi liquids 
have been discussed by Oganesyan, Kivelson, and Frad- 
kin (2001). The relation to nematic liquid crystals has 
been pointed out (Kivelson, Fradkin, and Emery, 1998). 

In three dimensions in channels with even angular mo- 
mentum, the Pomeranchuk instability is generically of 
first order due to the presence of cubic terms in the 
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Ginzburg-Landau theory. It has been argued that even 
in = 2 strong fluctuations may drive the transition 
first order, thus avoiding critical quantum fluctuations 
(Kee, Kim, and Chung, 2003; Khavkine et al., 2004). 
Full quantum critical behavior is restored in these 2d 
models, if a sufficiently strong repulsive term is added to 
the forward-scattering interaction (Yamase, Oganesyan, 
and Metzner, 2005). 

In the neighborhood of a Pomeranchuk instability 
the electron system shows unusual properties due to a 
"soft" Fermi surface, leading to a strongly enhanced de- 
cay rate for single-particle excitations and non-Fermi- 
liquid behavior (Metzner, Rohe, and Andergassen, 2003). 
The dynamical Fermi-surface fluctuations near a Pomer- 
anchuk instability in d = 2 have been analyzed re- 
cently (Dell'Anna and Metzner, 2006): the electronic self- 
energy scales as w^/^, thus destroying the Fermi liquid at 
all wavevectors. 

Here we sketch the calculation of the critical behavior 
- within Fermi-liquid theory near a spin-symmetric, 
d-wave Pomeranchuk instability, for which the dom- 
inant Fermi liquid interaction component is fuw = 
TTTf-F^'Smd* -d , where d is one of the £ = 2 eigen- 
functions. Following Wolfle and Rosch (2007) the corre- 
sponding susceptibility is of the form 

x£k'(q,^) = d*^^d^^,^2NoSM,^) (119) 
S-\<i,iv) = (W0'+4V-/3(q,c^)-n(q,w) , 



where 



(WO' = -7(1 +^275) 



(120) 



in d = 3 and = ^{1 + F^) in d = 2 with 



^nO(,,o) 



, 7(q,a;) = — Imn2(q,a;) 
g=0 m* 



and 



n°d(q,u;) 



(2^ 



fO 



iO - ek+ + fk- 



where ek± = ek±q/2- Strictly speaking, the formula for 
is only a crude estimate as this high-energy property can- 
not be calculated within Fermi liquid theory (Wolfle and 
Rosch, 2007), in contrast to the low- frequency damping 

It is important to distinguish two types of modes: 
"even" modes (with £ + m even) have 7(g, w) = ^, where 
c is the bare Fermi velocity, and /3{q,u)) =0 . For "odd" 
modes, 7(g,^0 = and P{q,u;) = 2^{^f. 

Consequently, for even modes we have a dynamical crit- 
ical exponent = 3, whereas for odd modes the 'bare' 
dynamical critical exponent (ignoring mass rcnormaliza- 
tions) is z = 2. In the Fermi liquid regime, the modes 
with the highest value of z will dominate. The role of 
the sub-dominant odd modes and their dynamical crit- 
ical exponent z depends on the scaling of the prefactor 



^ in (3 . We have here a situation of multiple critical 
exponents, as discussed in Sec. III.H.l. 

The correlation length ^ diverges in the limit — > —5. 
We will be interested in calculating the quasiparticle ef- 
fective mass and the specific heat as well as the contri- 
bution to the electrical resistivity caused by scattering 
from critical fluctuations. These quantities may be ex- 
tracted from the quasiparticle self-energy 5](k, u), follow- 
ing Dell'Anna and Metzner (2006), with imaginary part 
given by (for a lattice system with only one critical mode) 



ImS(k,w) 



m 



mk I 



duj' 



d'^q 
{2nY 



X — Im Sd(q, io')5{J + io - a+q)- (121) 



In an isotropic system a difii'erent averaging has to be 
used (Wolfle and Rosch, 2007), and for the dominat even 
modes one can replace [(i^kl' 1- After performing 
the integration over frequency w' and the component of 
momentum qr, where q = Q'^ki? -|- qj, and rescaling q^ = 
{u/vf)qr, qt = {i^/ilcY/^qt, one finds 



S(k,a;) = 



2A^o (2^) A ^2 J 
c j dqtqf^n 



d/3 



(122) 



where ( = ^^^(^q/w)^/^. As mentioned above, at the 
critical point {( 0) ImS(fc,a;) oc w''/^. 

The contribution to the effective mass from crit- 
ical fluctuations is obtained from ^ = 1 — 

^£ReE(fc,a;)|,=oas 



m 



(123) 



(The factor of m* /m converts a quasiparticle self energy 
into an electron self energy.) Since the specific-heat co- 
efficient C/T (X m* /m we have as well C/T oc C^"'', in 
agreement with the result obtained within a bosonic de- 
scription of the critical dynamics, see Eq. (67): C/T oc 
(C/Co)^^'', considering that z = 3 in this case. 

The contribution to the resistivity is found from Im E 
by the following qualitative argument: Since the typical 
momentum transfer in electron scattering off a critical 
fluctuation is Ag ~ the weight with which such a 
process contributes to the resistivity is reduced by Ag^ oc 
1 — COS0, where 6 is the scattering angle. As Eq. (122) 
gives ImE oc C^"'', we find for the scattering rate 1/r oc 
and hence 



Ap a e-^iT/To) 



(124) 



Note that a related calculation could be done for a fer- 
romagnetic instability; however, in ferromagnets further 
complications arise, as described below in Sec. III.H.l. 
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H. Breakdown of the Hertz model of a magnetic QCP 

Under what circumstances does the theory of Hertz 

break down? In this section several possible mechanisms 
for such a failure of the Landau-Ginzburg- Wilson (LGW) 
approach in clean systems are highlighted; disorder is 
briefly discussed in Sec. III. J. 

Conceptually, two causes for a failure of the Hertz the- 
ory can be identified: (i) A local analytic expansion of 
the action in terms of the magnetic order parameter does 
not exist, or (ii) additional degrees of freedom other than 
magnetism become critical at the transition. 

The first situation may arise when, in addition to 
the order-parameter fluctuations, other (fermionic) slow 
modes are present in the critical system (as is always 
the case in metallic magnets). Upon integrating out the 
fermions, non-analytic non-local or even singular terms 
may arise, invalidating the approach of Hertz. Then the 
whole concept, namely to consider an effective descrip- 
tion in terms of the critical modes alone, fails - examples 
to be discussed below are the metallic FM in d < 4 (Be- 
litz, Kirkpatrick, and Rollbiihler, 2004) and the metallic 
AFM in d < 2 (Abanov and Chubukov, 2004). A proper 
critical theory should include both order-parameter and 
fermionic modes, but such a coupled RG treatment has 
only been performed in a few cases, see Sec. III.H.l. Note 
that even in cases when a local expansion of the critical 
theory in terms of the order parameter is justified, one 
has make sure that the standard Fermi-liquid form (e.g. 
the \lu\ term from Landau damping) applies. 

The second situation may apply to certain heavy- 
fermion systems and will be discussed in Sec. III. I. 



1. Multiple dynamical exponents: FM QCP 



is the spin density of the electrons, expressed in terms 
of field operators 5'c((r) (Kirkpatrick and Belitz, 1996b; 
Vojta et ai, 1997). If So[^] is the (imaginary-time) ac- 
tion of the Fermi liquid in the absence of the exchange 
interaction Hj, then one can rewrite the partition sum 
formally as a functional integral over the collective field 
(/>(r, r) with the Hubbard-Stratonovich transformation: 



(125) 



with 



( Ul"+1 



X^l </'q,a;„, (126) 



(127) 



where x''"'' are the (connected) n-point spin susceptibil- 
ities of the reference system ^o, i.e., the susceptibilities 
in the absence of the exchange interaction. For a clean 
Fermi liquid in d = 3, Belitz, Kirkpatrick, and Vojta 
(1997) found in perturbation theory in the interactions 
(of the spin-singlet and Cooper channel) 



y(2) 



2 I / 2/1 -'■ r^" 
ci—C2q +C2qin C3 

q q 

U4 — V4£n- 

q 



(128) 
(129) 



The metallic ferromagnet is an example where the 
LGW approach of Hertz fails due to the presence of 
fermionic modes in the system. The idea can be discussed 
in terms of time scales: In a nearly-critical quantum sys- 
tem, the length scale ^ may induce several diverging time 
scales. The order parameter fluctuates on the time scale 

OC with, e.g., zqp = 3 in a clean FM according to 
the theory of Hertz (78). A different time scale is induced 
by the fermions. In a clean system, electrons cross an or- 
dered domain of size ^ ballistically in the much shorter 
time ts oc = 1. In a disordered system, charge 

or spin (if conserved) diffuse over a distance of ^ in the 
time Id oc zd = 2. 

Vojta et al. (1997) [see also Belitz et al. (2000, 2001a); 
Belitz, Kirkpatrick, and Rollbiihler (2004)] have shown 
that for itinerant quantum critical ferromagnets these 
other slow modes are indeed important (for AFM the 
effect is less severe, see Sec. III.H.2.) The problem 
becomes apparent when deriving the LGW functional 
from a microscopic theory in perturbation theory. Con- 
sider a system with an exchange interaction Hj = 
-//S(r)S(r')J(r-r'), where S(r) = M/t,(r)a„a'*a'(r) 



or for arbitrary n, x''"^ ~ \/q^~'^~^. In general, x'"^ 
is a complicated non- analytic function of [n — 1) mo- 
menta and frequencies; we give only the leading sin- 
gularity for certain limits qn,ijJn 0. The presence 
of these non-analytic corrections in the spin suscepti- 
bility of a clean Fermi liquid has later been verified by 
Chubukov and Maslov (2003, 2004); interestingly, these 
non-analyticities are absent in the charge channel. Fur- 
ther we note that related non-analytic corrections also 
show up in the finite-temperature behavior of a Fermi 
liquid, see Sec. II. D. In a diffusive system in d = 3, 
a result similar to (128) holds, with even stronger non- 
analyticities (Kirkpatrick and Belitz, 1996b): 

X^l - ci-C2g2-c^|g|-C3^, (130) 

X^^^ ~ U4 + V4^ (131) 

and x'"^ l/g2n-2-d_ 

The Gaussian action 5*2 close to the QCP is therefore 
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given in leading order by 



^2 - ^E('^ + -^^^'--2'?^^+-3^)|0, 



(132) 



for a clean or diffusive metal, respectively. 

Does the Gaussian action correctly describe the QCP 
in leading order? To decide this question we proceed with 
a scaling analysis. We choose the €3 term in 6*2 and S2 to 
be dimensionless and therefore 



oc f 2 in the clean 

cx in the diffusive metal in d = 3. The 

interaction Sn in d = 3 therefore scales proportional to 



and 



(X 



St OC ^2„-5^1±i„^-(„-l)(3+.) ^^-(f-l)(.-2) (134) 

in a clean or diffusive metal, respectively. 

What is the value of zl From the C3 term in the Gaus- 
sian action (132) one finds that at the QCP the order 
parameter fluctuates very slowly with w oc g^, and there- 
fore zop = 3 both for clean and dirty systems. From 
this argument, it seems that the contributions from the 
interactions (133,134) are irrelevant as they vanish for 
large ^, and therefore it was concluded by Kirkpatrick 
and Belitz (1996b); Vojta et al. (1997) that the critical 
theory is described by the Gaussian model 82- After- 
ward, the authors realized (Belitz et al., 2001a) that the 
simple scaling argument given above is not completely 
correct. The origin of this failure of "naive" scaling is - 
as we discussed in the beginning of this section - that 
other slow time scales with dynamical exponent = 1 
in the clean metal or z/j = 2 in the diffusive system are 
induced by the electrons. If we use these values for z, we 
find that all interactions S*^ and S'fj are marginal, and 
therefore at least the possibility exists that all of them 
have to be considered. 

As mentioned above, one has to conclude that a de- 
scription of the QPT in terms of the critical modes alone 
is not possible for the ferromagnet (Belitz, Kirkpatrick, 
and RoUbiihler, 2004). A generalized critical theory has 
been set up and analyzed for the ferromagnetic QPT in 
a diffusive metal (Belitz et al., 2000, 2001a) - this theory 
involves both the ferromagnetic order-parameter field (p 
and the diffusive modes of the disordered metal. Remark- 
ably, the Gaussian critical theory was found to be correct 
only up to logarithms, as suggested by the scaling with 
z = = 1 or z = ZD = 2, respectively. 

We proceed with a short discussion of the main phys- 
ical consequences of the non-analytic corrections to the 
Hertz theory in quantum critical itinerant ferromagnets. 
In the clean metal, the logarithmic corrections in S2 and 
S% lead probably to an instability of the ferromagnetic 
second-order transition: the 0^ term changes sign at an 



order 




1st order 

FIG. 4 Generic phase diagram of an itinerant ferromagnet, as 
function of temperature T, tuning parameter p, and magnetic 
field h. PM (FM) denote the paramagnetic (ferromagnetic) 
phases, and TCP is a tricritical point. From Belitz, Kirk- 
patrick, and RoUbiihler, 2005. 



exponentially small temperature inducing a weak first- 
order transition (Belitz, Kirkpatrick, and Vojta, 1999). 
Such a first-order transition has, e.g., been reported in 
ZrZn2, see Sec. IV. B. 3. This generic scenario has been 
discussed by Belitz, Kirkpatrick, and RoUbiihler (2005), 
and the resulting phase diagram is shown in Fig. 4. The 
—q^inl/q correction in S2 can also drive the transition to 
an AFM QCP (see Vojta and Sknepnek 2001). At some 
distance from the first-order or AFM transition, the pre- 
dictions of the Hertz approach for the 3d ferrmomagnet 
regarding resistivity and specific heat are likely to hold 
(up to logarithmic corrections), with a resistivity cx T^/^ 
at the QCP and logarithmically diverging specific-heat 
coefficient C/T. The exponents, e.g., for the pressure 
dependence of the Neel temperature (94) (Millis, 1993) 
close to the QCP (but not too close to the first-order 
transition) require more careful considerations ~ outside 
the scope of this review - as they involve directly the 
scaling dimension of the cj)'^ term. Chubukov, Pepin, and 
Rech (2004) have critically studied whether self-energy 
effects and vertex corrections wash out the singularities 
that lead to the first-order transition. Interestingly, they 
found that such corrections can in principle modify crit- 
ical exponents, but in the case of the ferromagnet the 
transition remains of first order. 

In the disordered diffusive ferromagnet, the Gaussian 
fixed point is stable, but the exponents deviate strongly 
from mean-field behavior (Belitz et al., 2000). The quan- 
tum critical behavior shows up, e.g., in the resistivity. 
The Vt cusp in the resistivity of a dirty Fermi liquid is 
modified and a T^/^ temperature dependence (Belitz et 
al, 2000) with logarithmic corrections is expected. Simi- 
larly the tunneling density of states should display a uj^^^ 
anomaly. More details can be found in Belitz et al. (2000, 
2001a). 
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2. Infinitely many marginal operators: AFM QCP in d = 2 

As = 2 is formally the upper critical dimension 
of the spin- fluctuation theory for the antiferromagnet, 
one would expect logarithmic corrections to the mean- 
field behavior as described above. However, Abanov and 
Chubukov (2004) have shown that in d = 2 the deriva- 
tion of the LGW theory breaks down. Somewhat simi- 
lar to the FM case, the low-energy modes of the Fermi 
liquid lead to long-range order-parameter interactions. 
Formally, the coefficients of the high-order interactions 
in the LGW functional diverge, leading to an infinite 
number of marginal operators. An analysis of the re- 
sulting theory is difficult and again requires a treatment 
of a coupled field theory of order-parameter fiuctuations 
and fermions - this has not been done until now. Abanov 
and Chubukov (2004) conclude that the 2d metaUic AFM 
shows a continuous transition with non-trivial exponents, 
but concrete predictions, e.g., for transport are lacking. 

3. Self-energy effects close to QCP 

The scattering from spin fluctuations strongly modifies 
the quasiparticles close to the hot lines et^j = ekH±Q in 
the vicinity of an AFM QCP. In leading order perturba- 
tion theory, the self-energy of those electrons at T = is 
given by 

At 

ImEk(n) ^ giy2 ^ ImXk-k' (w) Im^f^, (w - Q), 
^ Jo 

(135) 

where gs is the vertex of the coupling of electrons to spin 
fiuctuations and g^,{(jj) « l/(^^ — ek + ^O"*") is the Green 
fimction of the (free) fermions. Using 1/xc^±q(uj) ~ q^-|- 
(iw)2/^op and (135) we obtain at the AFM QCP 

ImSk„+,k(0) ~ O^+^/I^Jgl) (136) 

where (5k ~ 5k • Vk^+q is a measure for the distance 
from the hot line, and / is some scaling function with 
f{x ^ 0) ^ const and f{x oo) ^ l/x^~ . 

Obviously a Fermi-liquid description of the electrons, 
which requires ImSk(ek) < Ck for fc — > kp, breaks down 
for d < 3 for quasiparticles with momentum k//. In this 
sense, the critical dimension for a breakdown of FL the- 
ory is d = 3, even within the Hertz approach (which 
formally has d+ = 2). For quantum critical antiferro- 
magnets this breakdown of Fermi liquid for rf < 3 aff'ects, 
however, only a tiny fraction of the quasiparticles. At 
present it is not clear whether this effect will influence 
the spin dynamics. The RG analysis of the theory of 
Hertz (77) suggests that this is not the case, however, this 
question should be addressed in a RG treatment which 
includes both fermionic and bosonic modes. 

Let us contrast this scenario with the Q = situation 
of a metallic ferromagnet or a Pomeranchuk instability. 



At the QCP the susceptibility is of the form l/xq(cj) ~ 
q°' + {iuj)/q^°^~°' with a = 2 and zqp = 3 (supplemented 
by logarithmic corrections in 3d). Then, the self-energy 
is momentum independent, 

ImSk(O) ~ fl^+''~-oT\ (137) 

leading to a marginal behavior in d = 3, and a breakdown 
of the Fermi liquid in d = 2 due to the w^/^ dependence 
of the self-energy, here over the entire Fermi surface due 
to Q = (see also Sec. III.G). Then, it is natural to 
suspect that the precondition for the Hertz theory are 
no longer fulfllled. Chubukov, Pepin, and Rech (2004) 
have investigated in detail the role of self-energy effects 
for ferromagnetic quantum transitions. Indeed, those ef- 
fects arc found to be relevant and would modify critical 
exponents, if the transition were not of first order as con- 
cluded by the authors, see Sec. III.H.l. 

4. Pseudogaps close to QCP 

In the antiferromagnetically ordered phase of a metal, 
gaps open in parts of the Fermi surface with ek « CkiQ ~ 
0, provided that the Fermi volume is sufficiently large. 
(As already mentioned in Sec. III.C.l, this feature is not 
captured by the LGW order-parameter theory of Hertz 
- this theory is not vafid in the ordered phase, e.g., it 
wrongly predicts damped Goldstone modes.) It is there- 
fore important to ask whether precursors of these gaps 
will show up already in the paramagnetic phase close to 
the quantum critical point. In this regime, the behavior 
of the system is dominated by large antiferromagnetic 
domains of size ^, slowly fluctuating on the time scale 

~ ^^op where zqp is the dynamical critical exponent 
of the order parameter. As ^ is diverging when the QCP 
is approached, it is suggestive (Schrieffer, 1995) to as- 
sume that the electrons will adjust their wave functions 
adiabatically to the local antiferromagnetic background 
and will therefore show a similar behavior as in the AFM 
ordered phase. Will precursors of this effect show up and 
induce pseudogaps in the paramagnetic phase for suffi- 
ciently large ^? This would imply a breakdown of the 
LGW approach of Hertz. Below we will try to estimate 
this effect using a simple qualitative scaling analysis, for 
details see Rosch (2001). Pseudogaps play an important 
role in the physics of underdoped cuprates (Anderson, 
1997) and it has been speculated that they are indeed 
precursors of gaps in either superconducting, antiferro- 
magnetic, flux or striped phases (not discussed in this 
review). 

To define the concept of a pseudogap more precisely, 
we first analyze the ordered phase where a proper mean- 
field Hamiltonian of the electrons is of the form 

(138) 
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A is proportional to the staggered order parameter (as- 
sumed to point in z direction) and the k sum extends 
over a magnetic BriUouin zone. Close to the "hot points" 
("hot lines" in three dimensions) with = ek/f±Q = 
a gap opens and the band structure at k = k/f + (5k is 
approximately given by 

4 « ^^^^^-^k ± i ^[(vi - V2)5k]^ + 4A2 (139) 

where vi = Vkj, and V2 = Vk,, +q are the Fermi velocities 
close to the hot points. Interactions will actually induce 
some small weight within these gaps but this does not 
invalidate the mean-field picture. 

Within the LGW theory (77) of the AFM QCP (Hertz, 
1976; Millis, 1993) no precursor of the gap shows up, since 
spin-spin interactions are irrelevant by power count- 
ing. In the discussion of the ferromagnetic QPT in 
Sec. III.H.l, we have seen that these arguments arc not 
reliable in the case of a quantum phase transition in a 
metal, as generally a second dynamical critical exponent 
zf exists which describes that, e.g., ballistic fermions can 
traverse a domain of size ^ in a time tf oc with 
Zf = 1- 

For the following argument, we assume that the sus- 
ceptibility at the QCP has the form l/Xq±Q('^) ~ Q'^ + 
(2Cj)^/^°^. We are mainly interested in the case zqp = 2, 
smaller values for zqp might be relevant if pseudogap for- 
mation takes places, larger values have, e.g., been used to 
fit experiments (Schroder et al, 1998) in CeCu6_j;Auj;. 
For our argument, we assume that the quasiparticles 
move in a (quasi-static) staggered field with the effec- 
tive size A (to be determined later). According to the 
mean-field result (139) a gap of size lu* = A opens in a 
(d — 2) dimensional stripe in momentum space of width 
k* = A/vf- Interactions can change this (see (136)) to 
u* ~ {k*y'' ~ A^^ where Zf = ^ is the mean-field expo- 
nent. Heiscnbcrg's uncertainty relation dictates that the 
electrons have to see a quasi-static AFM background for a 
time T* > l/w* on a length scale of order ^* > 1/fc* per- 
pendicular to the direction of the hot lines to develop the 
pseudo gap. What is the effective size of the quasi-static 
AFM order on these length and time scales? The 

following estimate should at least give an upper limit at 
the QCP: 

{{<t>)f,r') < dco d-^-^ ImXqiQ (140) 

Jo Jq±<k* J-oo 

where the anisotropic integration of q takes into account 
that the momentum of the electrons parallel to the hot 
line can vary on the scale kF- In (141) we assumed zf < 

ZOP- 

If we assume furthermore that A is proportional to 

((/))|?^. as suggested by the mean field analysis (which 
should be valid above the upper critical dimension), we 



obtain the inequality A^ < const- A^'''''^°^^''-'^op+^. This 
implies that, at least in a weak-coupling situation, pseu- 
dogaps can appear only for 

d + zop<4:. (142) 

Note that it is accidental that Eq. (142) coincides with 
the condition for the relevance of the 0^ interaction (79) 
in the Hertz model as is evident from the fact that Zf 
enters Eq. (141). Within the approach of Hertz, zop = 
2 and the critical dimension for pseudogap formation is 
therefore dc = 2, i.e., no pseudogaps are expected in 
d = 3 as long as interactions arc not too strong. 

The derivation of Eq. (142) is based on a number of 
assumptions. The estimate (140) of {(p)^^ and therefore 
(142) is based on the existence of amplitude fluctuations 
of the staggered order parameter which destroy the pseu- 
dogap. Electrons can adjust their wave functions much 
better adiabatically to angular fluctuations of the direc- 
tion of the staggered magnetization than to fluctuations 
of its size. Schrieffer (1995) has argued that pseudogap 
behavior will occur always sufficiently close to an AFM 
QCP. For his argument, he considered models without 
amplitude fluctuations. Within the theory of Hertz, how- 
ever, amplitude fluctuations are present close to the AFM 
QCP md = 3 because the system is above its upper crit- 
ical dimension. Furthermore, strong statistical interac- 
tions of the electrons with the magnetic excitations (see 
Rosch (2001) and references therein) might destroy pseu- 
dogaps even in the absence of amplitude fluctuations. 



5. Itinerant AFM with Q = 2kF 

In our previous discussions of QCP of nearly AFM met- 
als we have assumed a large Fermi volume with Q < 2k f, 
where the spin fluctuations couple directly to quasipar- 
ticles with Ckfj ~ Ckff-i-Q ~ cf along "hot" lines on the 
Fermi surface {d = 3). In the opposite case, Q > 2kF, 
the spin fluctuations decouple in leading order from the 
quasiparticles due to energy and momentum conserva- 
tion and their dynamics follows from the conventional 
term as in magnetic insulators, Eq. (63). 

A special case is Q = 2kF- In d = 3 the hot lines 
shrink to a single point with parallel Fermi velocities 
VknW'^kH+ci for ek„ = einH+Q = ^f- In this situation, 
resonant scattering of the spin fluctuations from the elec- 
trons leads to a complete breakdown (Millis, 1993) of the 
LGW expansion (77) underlying the Hertz theory of a 
QCP. This can be seen from a direct calculation of the 
connected n-point function x" (127) in the limit where 
all momenta are set to Q and all frequencies to 0. In this 
limit, the effective n-paramagnon interaction diverges for 
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even n in the low-T limit 



ri/2 



ZWn — Ck+Q 



l/T"-4 forrf = 3 
l/T"-5/2 ford = 2 



n/2 



(143) 



Our result in d = 3 differs from the formula of Millis 
(1993); we find that the contribution ~ 1/r"-^ van- 
ishes exactly. A simple scaling analysis with k ~ 1/^, 
w ~ T ~ 1/L2, (/)(r,T) - ii-(d+2)/2 siio^s that S'„ in 
(127) diverges with 5„ ~ L"/2-3 in d = 3 and 5„ ~ 
in d = 2. Interactions of arbitrarily high n are therefore 
relevant and the LGW expansion in terms of the order 
parameter (77) breaks down completely. A critical the- 
ory cannot be formulated in terms of order-parameter 
fluctuations alone. 

However, a weak-coupling analysis suggests that Q = 
2kp is not realized in generic three-dimensional systems. 
The reason is that in (126) the polarizability of non- 
interacting fermions 



xf(-=o)=i: 



/(ek) - /(£k+q) 

Ck — Ek-l-q 



where /(w) is the Fermi function, does not peak at 
q = 2fci?. This is not only true for a quadratic disper- 
sion, but also for any band structure in the absence of 
perfect nesting. Whether interaction effects can stabilize 
an AFM QCP with Q = 2kF in the absence of perfect 
nesting is not known. 

In d = 2, however, a spin-density wave transition with 
Q = 2kF is very likely, as the polarizability of non- 
interacting electrons Xq is typically peaked at 2kp. Alt- 
shuler, loffe, and Millis (1995) analyzed such a situation 
and concluded that the strong interactions will probably 
induce a first-order transition. In this sense, the QCP is 
destroyed. 



6. Superconductivity 

Generically, the spin fluctuations induce an attrac- 
tive interaction between the quasiparticles (Abanov, 
Chubukov, and Finkelstein, 2001; Monthoux and Lon- 
zarich, 1999). Accordingly, one can expect a supercon- 
ducting phase close to a magnetic QCP of a metal as it 
is observed in sufliciently clean samples, see Sec. IV. C. 
Superconductivity is outside the scope of this review; 
here we just note that the order parameter is proba- 
bly unconventional and that the superconducting phase 
will change the dynamical critical exponent zq^ and 
therefore the critical behavior of the antiferromagnctic 
QCP due to a suppression of the spin-wave damping in 
the presence of gaps. In this sense the LGW theory 
(77) of Hertz breaks down due to superconductivity. In 
d = 3, the superconducting phase appears typically at 



very low temperature (Sec. IV. C), and the Hertz theory 
remains valid at temperatures above the superconducting 
Tc- The situation may be different in (quasi-)2d systems, 
where Cooper pairs form at much higher temperatures 
(Abanov, Chubukov, and Finkelstein, 2001; Monthoux 
and Lonzarich, 1999), effectively reducing the damping 
of spin fluctuations. 



I. Breakdown of the Kondo effect in heavy-fermion metals 

For heavy-fermion systems (HFS) it is generally ac- 
cepted that the magnetic transition is driven by a com- 
petition of the lattice Kondo effect, which favors a param- 
agnetic ground state, and a magnetic RKKY or superex- 
change interaction between the local / moments (Do- 
niach, 1977). In the heavy Fermi-liquid state the local 
moments contribute to the Fermi volume leading to a 
"large" Fermi volume, see Sec. II. F. 2 - the / electrons 
are usually termed "delocalized" in this situation. 

Before discussing scenarios for the magnetic transition 
we have to think about the nature of the ordered phase in 
HFS (the state on the l.h.s. of Fig. 1). Two distinct types 
of magnetically ordered metals appear possible, (i) The 
magnetism can arise from a spin-density wave instal)ility 
of the parent heavy FL state. Here, Kondo screening is 
essentially intact, with a weak polarization of the local 
moments which are still "delocalized". We will refer to 
such a state as the SDW metal, (ii) A different kind 
of magnetic metal is possible where the localized mo- 
ments order due to RKKY exchange interactions, and do 
not participate in the Fermi volume, i.e., Kondo screen- 
ing is absent. We will denote this second state, which 
can be expected to be realized deep in the ordered phase 
(Yamamoto and Si, 2007), as a "local- moment magnetic 
(LMM) metal" . The distinction between these two kinds 
of states can be drawn sharply, if the Fermi surfaces have 
different topologies (albeit the same volume modulo that 
of the Brillouin zone of the ordered state), such that they 
cannot be smoothly connected to one another. 

Returning to the transition from the paramagnet to the 
antiferromagnet, one possibility is that the Fermi liquid 
undergoes a transition to a SDW metal - this QCP is well 
described by the LGW approach of Hertz, Eq. (77). In 
this situation the local moments remain screened across 
the phase transition, i.e., a suitably defined lattice Kondo 
temperature stays finite at the QCP. The anomalous 
behavior close to AFM QCP in heavy-fermion systems 
like CeCu6-a:Aua; and YbRh2Si2 (discussed in detail in 
Sec. IV. A) - inconsistent with the Hertz scenario - has 
stimulated discussions about a different transition (Cole- 
man, 1999; Schroder et al, 1998, 2000; Senthfl, Sachdev, 
and Vojta, 2003; Senthil, Vojta, and Sachdev, 2004; Si, 
Smith, and Ingersent, 1999; Si et at, 2001, 2003): If the 
ordered state is a LMM metal, the transition to be con- 
sidered now involves the breakdown of Kondo screening 
(due to competing magnetic fiuctuations) , accompanied 
by an abrupt change of the Fermi surface. This is an 
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exciting scenario, as the complete collapse of the Fermi 
surface is in a sense the most drastic violation of the as- 
sumptions of the Hertz theory. No local order parameter 
can be defined, and the LGW approach fails. Instead, 
the criticality is carried by emergent degrees of freedom 
associated with the Kondo effect. 

Our present theoretical understanding of such tran- 
sitions is limited, and wc will describe a few theoret- 
ical approaches below. An obvious question then is: 
Can there be a continuous transition where the Kondo 
screening disappears concomitantly with the appearance 
of magnetic long-range order? This will be discussed in 
Sec. III. 1. 4. Parenthetically, we note that some materials 
show a first-order volume collapse transition at finite T 
(McMahan et ai, 1998) in contrast, we are interested 
here in a continuous transition at T = 0. 

A zero-temperature transition involving the break- 
down of Kondo screening ("Kondo transition") implies 
a collapse of the Fermi surface - in fact, we use this as a 
defining criterion for a Kondo transition. Experimentally, 
the collapse of the Fermi surface may be detected via pho- 
toemission or de Haas-van Alphen measurements, and 
transport properties like the Hall conductivity will show 
a jump upon crossing the transition at lowest T (Cole- 
man et al., 2001; Coleman, Marston, and Schofield, 2005; 
Si, Smith, and Ingerscnt, 1999). (At a SDW transition 
the Fermi surface evolves continuously, and the Hall coef- 
ficient displays a kink, but no jump; only at a magnetic- 
field driven transition docs the derivative of the Hall cur- 
rent with respect to the magnetic field jump.) As the 
Kondo transition is not associated with a single critical 
(fermionic) wavevector, one may expect critical fluctua- 
tions in an extended range of the reciprocal space. These 
qualitative theoretical considerations fit remarkably well 
some recent experiments: In CeCue-xAux (Sec. IV. A. 1) 
the susceptibility at the AFM QCP was found to obey 
1/X(q,w) ~ /(q) + {-ioj + aTY" (Eq. 149 below) with 
an anomalous exponent a ~ 0.8, obtained from fits to 
susceptibility measurements and inelastic neutron scat- 
tering data at various positions in momentum space in- 
cluding q = (Schroder et al, 1998, 2000). Momen- 
tum and frequency dependence "separate", this favors 
an interpretation in terms of a Kondo transition. For 
YbRh2Si2 (Sec. IV. A. 5) no neutron scattering data are 
available to date. However, a recent magneto-transport 
measurement may indicate a jump in the Hall coefficient 
at the magnetic QCP (Paschen et al., 2004); in addition, 
magncto-striction data of Gegcnwart et al. (2007) show 
the vanishing of several energy scales at the same QCP 
of YbRhaSia. 

Let us emphasize that the breakdown of Kondo screen- 
ing does not imply that the local moments are free to 
fluctuate at the QCP: The critical behavior will be man- 
ifested in anomalous power laws in the spin correlations, 
as shown explicitly, e.g., within the scenario of "local crit- 
icality" described below. Similarly, there will be no £n2 
entropy per spin at the QCP or in the quantum critical 
region. Thus the characteristic temperature T1/2, where 



the magnetic entropy equals 0.5/'n2, is not expected to 
go to zero at the QCP - this is in fact consistent with 
specific heat data on both CeCue-xAu^ (Fig. 8 below) 
and YbRh2Si2. 

In a simple scenario, where the low-temperature state 
of the conduction electrons with small Fermi volume is 
adiabatically connected to high temperatures, a char- 
acteristic signature of a breakdown of Kondo screening 
may be a shift of the maximum temperature in the resis- 
tivity, Tm, to lower temperatures upon approaching the 
QCP (as this may signal the crossover from "small" to 
"large" Fermi volume). However, concrete calculations 
for the transport crossover are lacking, and moreover this 
picture seems not to be supported experimentally: In 
CeCue-ajAu^;, Tm (for p measured along the a direction) 
decreases smoothly across Xc and vanishes at a; « 0.16 
(v. Lohneysen et al., 2002). 

We note that arguments have been put forward (Mae- 
bashi, Miyake, and Varma, 2005) for a generic breakdown 
of the single-impurity Kondo effect at an antiferromag- 
netic QCP; the consequences for lattice models have not 
been studied in detail. 



1. "Local" QCP within extended DIVIFT 

A first approach designed to capture the breakdown 
of the lattice Kondo effect due to magnetic bulk fluctu- 
ations employs an extension of the dynamical mean- field 
theory (DMFT) and has been worked out by Si, Smith, 
and Ingcrsent (1999): Si et al. (2001, 2003); Smith and 
Si (2000). It led to the proposal of a "local" QCP (to be 
made precise below), based on the idea that the break- 
down of Kondo screening is a spatially local phenomenon, 
i.e., it affects every spin of the underlying Kondo lattice 
independently. 

The starting point is the Kondo lattice model, where 
localized spins Sj couple to the spin density of conduction 
electrons at lattice site i, Sj = cl^aai3Ci0/2, with 

H = Y1 + J^^ S, • s, + Y,IijSi ■ Sj (144) 

ka i i,j 

where a direct spin-spin exchange term (7) has been 
added to the usual Kondo lattice model, see Sec. II. F. 2. 
While the usual DMFT maps the lattice problem to a 
single impurity in a fermionic bath, the extended DMFT 
(EDMFT) uses a mapping to a so-called Bose-Fermi 
Kondo model with both a fermionic and bosonic baths 
(represented by operators Cfc and hk)' 

H\oc = ^ Ekc\^Cka + JS-SO + 

k(7 

+70 ^ S • {hk + bLfe) + Wfcbibfe,(145) 

k k 

see also Eq. (56) in Sec. II.F.4. Within EDMFT, the 
Green's functions and susceptibilities of the lattice model 
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are approximated by l/gk('-^) ~ ^ ^ ^ ^('-^) s-iid 
l/XqC^^) ~ -fq + M{co), where and M{ui) arc the 

electron and b self-energies of the local problem. The 
free parameters E)., Uk and g are determined from the 
self-consistency condition that the local Green's function 
and susceptibility in the global and local model, (144) 
and (145), match. Formally, EDMFT can be justified 
within a certain d ^ oo limit (Smith and Si, 2000) but it 
may be used as an approximation to a finitc-dimcnsional 
system as long as the physics is not dominated by long- 
range spatial fluctuations. 

Within EDMFT, it is possible to describe situations 
where collective magnetic fluctuations destroy the Kondo 
effect. The Bose-Fermi Kondo model (145) is known 
to have a continuous QPT, due to the competition of 
the two baths, between a phase with Kondo screening 
and one with universal local-moment fluctuations - see 
Sec. 11.F.4. The QCP of the lattice model (144) is thus 
mapped - via EDMFT - onto the impurity QCP of 
Eq. (145), where the magnetic instability of the lattice 
drives the Kondo effect critical. At this "local" QCP all 
self-energies are momentum-independent, and the non- 
local dynamics of the magnetic fluctuations is Gaussian. 

A qualitative analysis of the EDMFT equations ind = 
2 (Si et al, 2003) shows that the logarithmic divergence of 
the local susceptibility at the QCP can cause a power-law 
behavior of M{co) at T = 0, M(w) = -/q -|- (-iw/A)", 
where Q is the ordering vector, A a cutoff, and a a 
non-universal exponent. A numerical solution of a sim- 
plified EDMFT (without fcrmionic self-consistency and 
with Ising magnetic symmetry) has confirmed this result 
(Grempel and Si, 2003; Zhu, Grempel, and Si, 2003). 
The impurity critical point has been shown to feature 
uj/T scaling in x- These results are in remarkable agree- 
ment with what has been found in the experiments of 
Schroder et al. (2000), see Eq. (149) below - in particular 
the anomalous exponent of the susceptibility is obtained 
as a « 0.72 (Grempel and Si, 2003) while the value from 
fitting the experimental data is a « 0.74. Let us point 
out that, in this theory, the occurrence of lo/T scaling, 
despite the non-local magnetic dynamics being Gaussian, 
is caused by the non-LGW character of the critical point, 
where the leading singularities are driven by local physics 
controlled by an interacting impurity QCP. 

As explicitly discussed by Si, Smith, and Ingcrsent 
(1999); Si et al. (2003), the collapse of the Kondo scale 
necessarily leads to a jump in the Fermi volume upon 
crossing the transition at T = 0. The destruction of the 
Fermi surface is expected to cause non-Fermi liquid be- 
havior in the resistivity at the QCP, however, concrete 
theoretical predictions for transport at finite tempera- 
tures are lacking to date. 

One issue in the DMFT description of bulk critical- 
ity is related to the zero-point entropy: Impurity crit- 
ical points generically display a finite residual entropy 
(Vojta, 2006a); this implies an extensive entropy for the 
bulk system (which would render the corresponding fixed 
point extremely unstable). In the case of the Bose-Fermi 



Kondo model (145) a reliable calculation of the impurity 
entropy is not available to date, however, it is conceivable 
that it vanishes in the limit d — > 2+, circumventing this 
problem. 

2. Fractionalized Fermi liquid and deconfined criticality 

A different approach to the breakdown of Kondo 
screening, without any assumptions about spatial lo- 
cality, starts by identifying the zero-temperature phase 
which arises when Kondo screening breaks down without 
the simultaneous onset of magnetic order (or other types 
of symmetry breaking) . As has been detailed by Senthil, 
Sachdev, and Vojta (2003), the resulting state is a para- 
magnet where the conduction electrons form well-defined 
quasiparticles on their own and the local moments are in 
a fractionalized spin-liquid state - this phase represents a 
metallic spin-liquid state and has been termed "fraction- 
alized Fermi liquid" (FL*). The spin-liquid component 
may be gapped or gapless, and may feature a secondary 
instability to an ordered state, see Sec. III. 1. 4. 

Technically, the transition from FL* to FL, at T = 
as a function of some control parameter like pressure, 
can be analyzed in slave-boson mean-field theory plus 
Gaussian fluctuations around the saddle point (Senthil, 
Vojta, and Sachdev, 2004). The mean- field Hamiltonian 
for the Kondo-lattice model (144) reads 

Hmf = ^<^kcl^Cka - XO X] (frafr'a + h.C.) 
k {rr') 

+ frJra - 6o ^ (4a/fea + ^-C.) (146) 

r k 

where = \ fra'^aa'fra' IS thc auxiliary-fcrmion repre- 
sentation of the local moments, and xo^ l^fi ^-nd are 
mcan-fielcl parameters. 

The QPT from FL* to FL involves a change from a 
small Fermi volume, containing only conduction elec- 
trons, to a large Fermi volume including all Kondo spins 

- on the mean-field level this is signaled by the conden- 
sation of the slave boson 60 measuring the hybridization 
between the c and / bands. (Beyond mean field, a com- 
pact gauge field needs to be introduced to implement the 
local constraint of the / fermions; this, e.g., suppresses 
a finite-temperature phase transition.) It is illuminat- 
ing to discuss the Fermi surface properties: Close to the 
transition the Fermi surface consists of two sheets with 
primarily c and / character, respectively. Approaching 
the QPT from the FL side, the quasiparticle weight on 
an entire sheet of the Fermi surface vanishes continuously 

- this illustrates how a discontinuous jump in the Fermi 
volume can happen at a continuous transition. Clearly, 
the transition is not associated with a specific critical 
(fermionic) wavevector, but it is also not spatially local, 
as all self-energies retain their momentum dependence. 

The critical theory of the FL-FL* transition can 
be derived starting from the slave-particle formulation. 
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FIG. 5 Phase diagram near the conjectured FL-LMM quan- 
tum transition of a heavy- fermion metal (see also Fig. 1). The 
primary phase transition, characterized by the breakdown of 
Kondo screening, is between the heavy Fermi liquid (FL) and 
fractionalized Fermi liquid (FL*). FL* is unstable at low en- 
ergies towards local-moment magnetism. Two distinct energy 
scales are present on the l.h.s., manifested by the differing ex- 
ponents by which Tn and Tcoh (or T*^^) approach the quan- 
tum critical point. Prom Senthil, Sachdev, and Vojta, 2005. 



Eq. (146), supplemented by a gauge field. Provided that 
the two Fermi surfaces do not overlap, the fermions can 
be integrated out, and one ends up with a theory for 
dilute bosons h coupled to a compact U(l) gauge field. 
The transition is tuned by the chemical potential of the 
bosons; it occurs at the (bosonic) wavevector (5 = 
and has dynamical exponent 2; = 2, it is thus above 
its upper critical dimension. The FL coherence tem- 
perature vanishes as the transition is approached from 
the FL side as Tcoh oc |r|. The specific heat acquires a 
singular contribution from gauge-field fluctuations with 
C/T ~ £n{l/T) in = 3, resembling the experimental 
result on CeCue-a^Auj,. A preliminary transport calcula- 
tion (Senthil, Vojta, and Sachdev, 2004), taking into ac- 
count the scattering of the critical bosons off gauge-field 
fluctuations, led to a resistivity p ~ l/in{l/T), incon- 
sistent with experiments. Interestingly, the decay of the 
bosons into particle hole pairs becomes possible above an 
energy E* which can be small if the distance between the 
two Fermi surfaces is small; above this energy the the- 
ory obeys z — 3, and an additional in{l/T) contribution 
in C/T appears (Paul et al, 2007). In this regime, the 
resistivity has been estimated as p ~ T£n{T). Coleman, 
Marston, and Schofield (2005) have calculated the T = 
Hall coefficient using the model (146), and found a jump 
when passing through the QPT. Clearly, more detailed 
transport studies, also taking into account impurity scat- 
tering, are required. 

So far, magnetism is not involved in this scenario. 
Clearly, the spin liquid in the FL* phase is potentially 
unstable towards magnetic order at low T the resiilt- 
ing state will be a LMM metal. A particularly appeal- 
ing scenario is that this instability arises as a secondary 
one. driven by an operator that is irrelevant at the QCP 
(Fig. 5). This naturally leads to a Landau-forbidden 



transition within the concept of deconfined criticality, 
discussed below in Sec. III.I.4. We caution, however, that 
explicit calculations using slave-particle theories, taking 
into account magnetism, result in two separate critical 
points for magnetism and the Kondo effect, at variance 
with Fig. 5 (Senthil, Vojta, and Sachdev, 2004). 



3. Spin-charge separation at the QCP 

A related scenario for the breakdown of Kondo screen- 
ing has been proposed by Pepin (2005). It is based on 
the idea that the heavy quasiparticle fractionalizes into 
a spinon and a spinless fermion % at the QCP (Coleman 
et al, 2001). 

Formally, the Kondo interaction of the Kondo lattice 
model is decoupled here with a fermionic field (in contrast 
to the slave boson in standard mean-field theory for the 
Kondo effect), and the dynamics of this fermion is key for 
the critical behavior. In contrast to the ideas sketched 
above, in the approach of Pepin (2005) the Fermi vol- 
ume does not jump, but evolves continuously through 
the QCP. Building on a number of phenomenological as- 
sumptions on the dynamics and dispersion of the x mode, 
it is possible to describe various properties of YbRh2Si2, 
like a T~^/'^ upturn of the specific-heat coefficient at low 
temperatures which is not reflected in transport measure- 
ments. However, the theory of Pepin (2005) is not able 
to describe the heavy Fermi-liquid state itself, as the x 
fermion cannot condense. 



4. One vs. two transitions 

For all the above approaches, two possibilities arise. 
Generically, the breakdown of Kondo screening will not 
occur at the same point as the magnetic QCP which is 
associated with the formation of long-range order. There- 
fore the Hertz theory (77) for an antiferromagnetic QPT 
seems to remain valid. The situation is different when 
the two QCP coincide (or are in close proximity). 

Within the EDMFT model, this happens in the case 
of two-dimensional magnetic fluctuations: the local sus- 
ceptibility at the magnetic QCP diverges which drives 
the Kondo effect critical, leading to a power-law behav- 
ior of M{lu). (In d = 3 the magnetic transition precedes 
a possible breakdown of Kondo screening.) Interestingly, 
the momentum dependence of the critical fluctuations 
in CeCue-xAux appears to be two-dimensional (Stock- 
ert et al., 1998), which also implies that the self-energy 
of the electrons is weakly momentum-dependent (Roscli 
et al., 1997). Thus, the EDMFT model reproduces a 
number of important aspects of the experimental data on 
CcCufi-^Aua;, but it remains to understand why the mag- 
netic fluctuations are two-dimensional (no obvious struc- 
tural reason for this behavior is evident), and whether 
this fact is generic. 

A different scenario, based on the idea of a Landau- 
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forbidden transition, has been proposed within the frac- 
tionahzcd Fermi-hquid concept (Scnthil, Sachdev, and 
Vojta, 2005). Provided that the primary transition is 
the Kondo breakdown (leading to locahzed / moments 
on the r.h.s. of the phase diagram in Fig. 5), magnetism 
can arise as a secondary instabihty of the FL* state. The 
RG flow is similar to that proposed in the scenario of de- 
confined quantum criticality (Senthil et at., 2004a, b), i.e., 
an operator which destabilizes the dcconfincd phase (in 
our case towards magnetism) is irrelevant at the critical 
point. A consequence would be the presence of two differ- 
ent energy scales on the magnetic side of the QPT: fluc- 
tuations associated to the Kondo effect (loosely speaking 
Fermi surface fluctuations) exist on a much higher energy 
scale than magnetic fluctuations; this is accompanied by 
rather weak magnetism (i.e. an anomalously small or- 
dered moment) close to the QCP. Then, the Ncel temper- 
ature, Tn, will vanish faster upon approaching the QPT 
than the temperature scale, Tcoh, at which well-defined 
quasiparticlcs appear at the large Fermi surface on the 
FL side (see Fig 5). A microscopic calculation verifying 
this proposal is not available to date; a "naive" slave- 
particle calculation (which is blind to the mechanism of 
deconfined criticality) yields two separate critical points 
(Senthfl, Vojta, and Sachdev, 2004). 



J. Disorder effects close to quantum phase transitions 

When dealing with real materials, the influence of 
static or quenched disorder on the properties of a quan- 
tum phase transition is an important aspect. Remark- 
ably, the effect of disorder is not completely understood 
even for classical phase transitions. 

In a theoretical description, quenched disorder can oc- 
cur in different ways: on a microscopic level, e.g., ran- 
dom site energies or bond couplings, or randomly dis- 
tributed scattering centers are possible. In an order- 
parameter field theory, disorder usually translates into 
a random mass tcirm for the order-parameter fiuctua- 
tions. Importantly, the quantum statistical description 
of a quantum problem with quenched disorder leads to a 
(d+2;)-dimcnsional field tli(X)ry with strongly anisotropic 
correlated disorder because the disorder is frozen in the 
time direction. In some cases, lattice eff'ects not captured 
by the field theory can be important, this applies, e.g., to 
all types of percolation problems. Moreover, disordering 
a quantum model can lead to random Berry phase terms 
which have no classical analogue, an example are diluted 
Heisenberg magnets. 

If disorder is added to a system which displays a con- 
tinuous (classical or quantum) phase transition, obvi- 
ous questions arise: (i) Will the phase transition remain 
sharp or become smeared? (ii) Will the critical behavior 
change? (iii) What happens in the vicinity of the "dirty 
transition" ? 



1. Harris criterion and fixed points 

To answer the first two of the above questions one has 
to investigate the stability of a critical fixed point with 
respect to a small amount of disorder. The Harris crite- 
rion (Chayes et al., 1986; Harris, 1974) states that dis- 
order will induce qualitative changes (i.e. is a relevant 
perturbation) ifvd <2, where u is the correlation-length 
exponent. Note that the Harris criterion is identical for 
classical and quantum phase transitions (i.e., d is not re- 
placed by d + z) because the disorder is frozen in the 
time direction. Within the LGW theory (77) of Hertz, 

= 1/2 for d+z > 4, therefore disorder is always relevant 
sufficiently close to the QCP in a system with quenched 
disorder for d = 2 and 3 (Harris, 1974; Kirkpatrick and 
Belitz, 1996b). 

Combining results of neutron scattering and thermo- 
dynamic measurements, one can use the Harris argument 
to obtain an order-of-magnitude estimate of the temper- 
ature and doping regime where disorder will affect the 
system. For our argument, we consider a material where 
a magnetic QCP is reached by doping, e.g. a ternary 
compound ABi^xCx with a QCP at a; = Xc, Xc < 1/2. 
On the non-magnetic side of the phase diagram, one 
can think of the system as consisting of many fluctu- 
ating domains of size ^ with a volume £f. The 
number Nq of impurities in such a domain is approxi- 
mately given by Nc = xV^/Vuc ~ cx{x — Xc)~'"^, where 
Vuc is the volume of the unit cell and the constant c « 
(xo — XcY'^V^{xo) /V\jc can be estimated from inelastic 
neutron scattering experiments at a doping xq ^ x,,. Ob- 
viously, fluctuates statistically with variance 
Therefore the typical fluctuations fa o f the doping x are 
of the order &x « x^/ ^/Wl « yf^Tfcix - Xry^l"^ . The 
Harris criterion is equivalent to the statement that dis- 
order is relevant if 5a; > |a; — a;c| or 

A/a~7c > |a;- Xcl^""^ (147) 

\i vd < 2, disorder changes the critical behavior for 
|a; — Xcl < bx* « [xc/c) '^-"'^ . This doping scale 8x* can 
be translated into a temperature scale T* below which 
disorder changes the thermodynamics at the QCP, e.g., 
by using the x dependence of the ordering temperature 
or of other relevant crossover scales. (However, an order- 
of-magnitude estimate along these lines for CeCue-xAu^ 
with a QCP at = 0.1 turns out to be inconclusive, but 
suggests that disorder could be important in the experi- 
mentally relevant regime.) 

What will happen if quenched disorder is relevant suf- 
ficiently close to the QCP? Recent work has shown that 
three possibilities exist: (a) Disorder leads to a new 
conventional (finite-disorder) critical point, with power- 
law behavior and exponents fulfllling the Harris crite- 
rion, vd > 2. An example is the rung-diluted bilayer 
Heisenberg magnet in d = 2 (Sknepnek, Vojta, and Vo- 
jta, 2004). (b) Disorder leads to a so-called infinite- 
disorder fixed point. Here, the dynamics is extremely 
slow, £n^T ~ S,'^ (replacing the conventional ~ S,^), 
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and the statistical distributions of obscrvables become 
very broad. Such a behavior has been estabhshed for the 
random quantum Ising model in d = 1 (see Sec. III. J. 2). 
(c) Disorder can destroy the sharp transition, replacing 
it by a smooth crossover. This interesting scenario is rel- 
evant for certain metallic magnets and will be discussed 
in Sec. III.J.3. 

Independent of the fate of the phase transition point 
itself, the third of the above questions is still to be an- 
swered: What happens in the vicinity of the transition? 
Interestingly, even at some distance from the QCP, the 
system can show power-law behavior (e.g. as a hmction 
of temperature) with non-universal exponents. These so- 
called Griffiths effects are discussed in the next section. 



2. Rare regions and quantum Griffiths singularities 

Disorder in a magnet usually suppresses the ordering 
tendency and thus changes the location of the phase tran- 
sition. In a parameter regime where the clean system 
would order but the disordered does not, one will find 
(arbitrarily large) regions that are accidentally devoid 
of impurities, and hence show local order, with a small 
but non-zero probability that usually decreases exponen- 
tially with the size of the region. These static disorder 
fluctuations are known as "rare regions" , and the order- 
parameter fluctuations induced by them as "local mo- 
ments" or "instantons" . Since they arc weakly coupled, 
and flipping them requires to change the order parameter 
in a whole region, the local moments have very slow dy- 
namics. Griffiths (1969) was the first to show that rare 
regions lead to a non-analytic free energy in the whole 
region between the transition points of the clean and 
disordered system, known as the Griffiths (or GrifRths- 
McCoy) region. A review has been recently given by 
Vojta (2006b). 

In generic classical systems Griffiths effects are weak, 
since the singularity in the free energy is only an essen- 
tial one. Near quantum phase transitions Griffiths sin- 
gularities are enhanced compared to the classical case as 
disorder is frozen in the time direction. Interestingly, the 
three cases listed above regarding the fate of the tran- 
sition yield different quantum Griffiths behavior as well 
(Vojta, 2006b): (a) In the vicinity of a finite-disorder 
fixed point the Griffiths effects lead to weak exponential 
corrections, (b) For infinite-disorder fixed points Grif- 
fiths effects are strong, and obscrvables display power- 
law singularities with continuously varying exponents, 
(c) If the rare regions become static, the transition is 
smeared, and conventional Griffiths behavior does not 
exist. Griffiths singularities occur in principle also on 
the ordered side of a QCP, but their signatures are much 
weaker. Notably, the cases (a/b/c) corresponds to sit- 
uations where the rare regions are below/at /above the 
lower-critical dimension of their ordering transition (Vo- 
jta and Schmalian, 2005) . In the following we sketch the 
physics of situation (b) which has been thoroughly inves- 



tigated for spin models. Situation (c) which is relevant 
for the damped order-parameter dynamics of metals is 
discussed in Sec. III. J. 3.) 

The random Heisenberg and transverse-field Ising 
models have been studied in detail md= 1, but more re- 
cently also in higher dimensions (Fisher 1995, Motrunich 
et al. 2000, Pich et al. 1998, Senthil and Sachdev 1996). 
A transparent physical picture emerges from a real-space 
RG analysis (Fisher, 1995; Ma, Dasgupta, and Hu, 1979; 
Motrunich et al., 2000). For a strongly disordered AFM 
Heisenberg chain, H = Y^JijSiSj, the decimation RG 
scheme proceeds as follows: In each step, the strongest 
bond is eliminated (i.e. frozen as a singlet), which induces 
a new coupling between the adjacent spins via second- 
order perturbation theory. The RG procedure follows 
the flow of the distribution of couplings P{J) upon suc- 
cessive elimination of spins and bonds. The renormaliza- 
tion scheme is based on a strong-coupling expansion and 
perturbation theory, it is valid if the strongest coupling is 
typically much larger than neighboring couplings, i.e., if 
the distribution P{J) is very broad. For the Heisenberg 
chain a typical initial distribution of disorder gets broader 
and broader: it flows towards an infinite-randomness 
flxed point, corresponding to a so-called random-singlet 
phase, and the method described above is asymptotically 
exact. For the transverse-field Ising chain, where a sim- 
ilar scheme can be applied, the two stable phases are 
conventional, but the flow to infinite randomness occurs 
at the zero-temperature phase transition point. The crit- 
ical dynamics at such an infinite-randomness fixed point 
turns out to be extremely slow, with IuS^t oc (so-called 
activated scaling), and the distributions of macroscopic 
observables become infinitely broad. 

If the QPT is controlled by such an infinite-randomness 
fixed point [situation (b)] the following picture emerges 
(Fig. 6): On the paramagnetic side of the phase dia- 
gram, a distribution of magnetic domains induces a dis- 
tribution of local susceptibilities Xi O'' typical energies 
Ai ~ l/xi with probabilities P(Ai < A) - A''/^', 
where is the dimension. The smallest possible en- 
ergy A,nin(i) in a region of size L is therefore given by 
P(Ai < Amin) ~ 1/L'' (as there are L'^ sites within this 
region) and Amin ~ . In this sense, z' is a dynamical 
critical exponent. However, one should keep in mind that 
minimal, typical, and average energies can be very differ- 
ent in this Griffiths regime. Note that the characteristic 
low-energy scale of a domain consisting of N spins, e.g. 
the tunnel splitting between different magnetic configu- 
rations, is exponentially small in N . The relevant domain 
sizes A'' ^ InA are therefore rather small, and the size of 
domains enters most physical properties only logarithmi- 
cally. From the distribution of energies one finds for the 
specific-heat coefficient and the average susceptibility 

X ~ cv/T ~ T'^/^'-i . (148) 

The divergence of the average non-linear susceptibility is 

even stronger, x^ni ~ T'^/^ x^ni '^^^ diverge even if x 
is regular. As the typical size of domain is of order tnN 
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FIG. 6 Schematic phase diagrams near a QCP in the presence 
of quenched disorder. Left; infinite-randomness fixed point 
with Griffiths region [situation (b)]. In the Griffiths region, 
thermodynamic quantities display power-law singularities in 
a finite region around the QCP. The dashed line indicates the 
divergence of the dynamical exponent z' upon approaching 
the QCP. Right: smeared phase transition [situation (c)], with 
exponentially small transition temperature on the disordered 
side of the clean system (r > 0). 



and therefore small, the leading temperature dependence 
of the order-parameter susceptibility and the static sus- 
ceptibility is typic;ally the; same. The exponents d/ z' in 
the Grifhths region are non-universal as they depend on 
microscopic details and the distance r from the QCP. For 
the infinite-randomness fixed points one finds z' ^ r""*, 
with universal exponents \E' and v. (For numerical values 
see Motrunich et al, 2000; Pich et al, 1998.) This char- 
acteristic dependence of the Griffiths exponents on the 
distance to the QCP is in our opinion the most impor- 
tant signature of the quantum Griffiths scenario. Exper- 
iments on UCus-ajPdj: may possibly provide an example 
of this dependence (Vollmer et al., 2000). 

Numerical simulations (Motrunich et al., 2000; Pich et 
al., 1998) suggest that infinite-randomness fixed points 
may not be restricted to systems in d = 1, raising the 
possibility that exotic critical behavior dominated by rare 
regions may be common to certain quenched- disorder 
quantum systems, in particular those with Ising symme- 
try. Recent investigations of 2d diluted antiferromagnets 
with Heisenberg symmetry indicate conventional critical 
behavior [situation (a)], and in addition an interesting 
interplay of quantum and geometric criticality at the per- 
colation threshold (Sandvik, 2002; Sknepnek, Vojta, and 
Vojta, 2004). 

Castro Neto, Castilla, and Jones (1998); Castro 
Neto and Jones (2000) have proposed that Griffiths 
singularities (and the related Kondo disorder sce- 
nario) can explain the anomalous behavior in cer- 
tain strongly disordered heavy-fermion systems. In- 
deed, de Andrade et al. (1998) fitted C/T and x 
of a number of systems over a certain temperature 
range with Eq. (148), e.g., Thi.^U^PdaAlg, Yi.^U^Pds 
or UCu5_2;Pdj:. More recently, the low- field ac sus- 
ceptibility of Ce(Rui_a;Rhj,)2Si2 was found to exhibit 
B/T scaling compatible with a quantum Griffiths sce- 
nario (Tabata et al., 2004), while "^^Si NMR data on 
CePtSii_xGex near a magnetic instability appear to dis- 



agree with this scenario (Young et al, 2004). 

As we will discuss in the next section, in metallic sys- 
tems with Ising symmetry the quantum tunneling of the 
magnetic domains - which is at the heart of the quan- 
tum Griffiths effect - is prohibited at lowest temper- 
atures by the coupling to the feniiioiis (Millis, Morr, 
and Schmalian, 2002). Whether the characteristic power 
laws are nevertheless observable in a certain temperature 
regime may depend on non-universal details and is not 
completely clear (Castro Neto, Castilla, and Jones, 1998; 
Millis, Morr, and Schmalian, 2002). 

Miranda and Dobrosavljevic (2001) succeeded in cal- 
culating Griffiths singularities close to a metal-insulator 
transition in a strongly correlated electron system. The 
authors attacked the challenging problem to describe 
both the formation of localized magnetic moments and 
the metal insulator transition using the so-called "sta- 
tistical dynamical mean-field" (SDMFT) approximation. 
These developments have been recently reviewed by Mi- 
randa and Dobrosavljevic (2005). 



3. Effects of rare regions on metallic QPT 

How do disorder and rare regions influence the QCP in 
a metal? In general, the situation is different from insu- 
lators due to the overdamped order-parameter dynamics. 
Early approaches, e.g., by Narayanan et al. (1999) who 
used a RG formulation in terms of the order parame- 
ter (Kirkpatrick and Belitz, 1996a) and found run-away 
flow to strong disorder, remained inconclusive regarding 
the nature of the transition. By now, only a few results 
are available, and the answer may depend on the order- 
parameter symmetry. 

Remarkably, for a metallic antiferromagnet with Ising 
symmetry the effects of disorder are so strong that a di- 
rect transition from the antiferromagnetic to a paramag- 
netic (Griffiths) phase is prohibited [situation (c)] (Vo- 
jta, 2003b). Generically, a spin-glass phase emerges at 
the lowest temperatures - for the specific model consid- 
ered by Vojta (2003b) the antiferromagnetic transition is 
smeared. This effect can be understood starting from the 
properties of a single rare region, a small antiferromag- 
netic domain, in the (nominally) paramagnetic phase. As 
disorder is frozen in time direction, a rare region in space 
translates into a rod-like object in space-time, which has 
an effective 1/t^ interaction arising from Landau damp- 
ing, i.e., the |w| term in Eq. (78). The effective one- 
dimensional Ising model for each rare region displays an 
ordered phase due to this long-range interaction, in other 
words, isolated rare regions will develop static order. In 
this situation, quantum Griffiths behavior (as described 
above) does not exist as there is no quantum tunneling. 
Isolated rare regions are coupled weakly by a RKKY in- 
teraction with random sign which leads to the formation 
of a classical spin glass (or rather "domain" or "cluster" 
glass), with the transition temperature being exponen- 
tially suppressed by the distance from the quantum crit- 
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ical point, sec Fig. 6. It is not known how the spin glass 
and antiferromagnetic phases merge; a Ukely scenario is 
a first-order phase transition. 

For continuous symmetries the situation is less clear. 
The rare-region physics of Vojta (2003b), applied to an 
itinerant magnet with Heisenberg (instead of Ising) sym- 
metry, has been shown (Vojta and Schmalian, 2005) to 
give rise to power-law Griffiths singularities near criti- 
cality, and the QCP is then likely of infinite-randomness 
type. However, Dobrosavljevic and Miranda (2005) have 
argued that the coupling between the rare regions can 
also load to a "freezing" of magnetic domains at very 
low T, destroying the Griffiths physics. Again this is 
expected to lead to a glassy state for T ^ 0. 



4. Metallic quantum glasses 

Strong disorder can lead to the formation of itinerant 
glass phases, with glassiness in the spin sector (metallic 
spin glasses, SG) or the charge sector (electron glasses). 
We will not discuss the physics of metallic glasses in detail 
and refer the reader to the theory reviews of Sachdev 
and Read (1996) and Miranda and Dobrosavljevic (2005) 
(for an experimental introduction see Mydosh 1993), and 
restrict the discussions here to a few short remarks. 

Unfortunately, not much is known theoretically about 
metallic spin glasses (SG) in 3d systems. While the de- 
scription of the ordered phase of metallic AFM or FM 
within Fermi-liquid theory is straightforward, a theory 
of a metallic SG has been developed mainly for infinite- 
range models (Sachdev, Read, and Oppermann, 1995; 
Sengupta and Georges, 1995). Close to the QCP of a 
metallic SG, one can expect both non-Fermi liquid be- 
havior due to critical fluctuations and Griffiths singular- 
ities (see section above and Sachdev 1998). Surprisingly, 
the QCP of a metallic SG is very similar to the d = 3, 
z = 2 theory of Hertz (77) (Sachdev, Read, and Opper- 
mann, 1995; Sengupta and Georges, 1995): the authors 
predict at the QCP a Vt correction to the specific-heat 
coefficient and a T^/^ law in the temperature dependence 
of the resistivity. We note that these mean-field theories 
do not satisfy hyperscaling. The actual situation for re- 
alistic spin glasses remains an open problem; an analysis 
of fluctuations around the mean-field solution leads to a 
runaway flow to strong coupling. 

Glassy behavior in the charge sector has long been in- 
vestigated in disordered insulators, where studies have 
concentrated on the physics of the Coulomb glass. Re- 
cent work has discussed the zero-temperature melting of 
such a glassy phase due to quantum fluctuations (Pas- 
tor and Dobrosavljevic, 1999). A mean- field theory for 
the transition from a Fermi liquid to an electron glass 
(Dalidovich and Dobrosavljevic, 2002) predicts a T^/^ 
behavior in the resistivity; corrections to the mean-field 
behavior have not been determined. 



IV. FERMI-LIQUID INSTABILITIES AT QUANTUM 
PHASE TRANSITIONS: EXPERIMENT 

As already mentioned in the Introduction, there are 
several material classes where non-Fermi-liquid behav- 
ior is observed. In the spirit of this review, we will fo- 
cus on quantum phase transitions and will not discuss 
single-ion NFL scenarios (see Sec. I.D). Suffice it to 
mention that, historically, the first 4/-electron NFL sys- 
tem was UxYi_a,Pdx with x « 0.2, where the observed 
specific- heat behavior C/T r-^ £n{To/T) was interpreted 
by Seaman et al. (1991) as arising from the two-channel 
quadrupolar Kondo effect, while essentially similar data 
were suggested by Andraka and Tsvelik (1991) to sig- 
nal a quantum phase transition. Further work by the 
Maple group underlined the metallurgical and disorder- 
related issues with this material by a comparative study 
of Ux^i-ajPda;, with A = Y, La and Sc (Gajewski, Chau, 
and Maple, 2000). A recent report indicated some scal- 
ing of the dynamical susceptibility x(w, T) (Wilson et al., 
2005) although the scatter of the data was substantial. 

From the beginning, there has been a lot of discussion 
about Kondo disorder leading to NFL behavior. Here 
UCug-xPda; was in focus (Bernal et al, 1995), which was 
the first material where a new type of anomalous scaling 
of the (wavcvector-integrated) susceptibility x(a;,T) was 
observed (Aronson et al, 1995). UCus-^Pda, has also 
been a candidate for the Griffiths scenario (Sec. HI. J. 2) 
(de Andradc et al, 1998) as corroborated by a depen- 
dence of the temperature exponents of x and C/T, 
Eq. (148), on the distance from the quantum critical 
points at X = 1 and 1.5 (VoUmer et al, 2000). How- 
ever, one has to bear in mind that the detailed behavior 
of this system is strongly affected by the Cu/Pd site in- 
terchange, the amount of which depends on annealing 
(Booth et al, 2002; MacLaughlin et al, 2006). For in- 
stance. Booth et al (2002) reported instead of a power 
law - a logarithmic temperature dependence of C/T in 
annealed UCu4Pd. 

Leaving these matters aside, we now briefiy describe 
the main focus of this section, i.e., experiments on mag- 
netic QPT in intermetallic rare-earth and transition- 
metal alloys and compounds. (Superconductivity will be 
discussed in so far as it is relevant to quantum criticality.) 
Rather than trying to be exhaustive in treating every sys- 
tem for which NFL behavior near a magnetic QPT has 
been claimed, we will discuss mostly systems that have 
been characterized thoroughly. These comprise hcavy- 
fermion systems, e.g., CeCue-xAu^, YbRh2Si2, and the 
Ce-115 compounds, as well as transition metal com- 
pounds such as MnSi and ZrZn2. A comprehensive re- 
view of NFL behavior in a large variety of heavy-fermion 
systems has been given by Stewart (2001, 2006). The is- 
sue of (one or several) QPT in cuprate high-temperature 
superconductors and other transition-metal oxides is be- 
yond the scope of this review. Of course, the linear resis- 
tivity p{T) in cuprates at optimal doping constitutes an 
early manifestation of NFL behavior (cf. the marginal 
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Fermi-liquid phenomenology of Varma et al. 1989). 

In discussing QPT, two generic cases have to be dis- 
tinguished. The prototype is a continuous, e.g., second- 
order transition between two ground states tuned by a 
non-thermal parameter. If the transition is discontinu- 
ous (first order), the evolution of the length and time 
scales fluctuations will be limited. An interesting case 
arises when a line of first-order transitions terminates in 
a critical cndpoint and if this critical point can be driven 
to zero temperature by another tuning parameter, lead- 
ing to a quantum critical end point. This might be the 
case for Sr3Ru2 07 where tuning parameters are magnetic 
field and the field orientation (Grigera et al., 2001; Millis 
et al, 2002). 



A. Quantum critical beliavior in heavy-fermion systems 

In general, Fermi-liquid theory has been very success- 
ful in describing the low-temperature behavior of met- 
als with strong electronic correlations, including many 
heavy-fermion systems (HFS). These materials often 
comprise a regular sublattice of 4/ or 5/ atoms, notably 
Ce, Yb, or U. Under certain conditions, a crossover of 
the magnetic behavior occurs with decreasing tempera- 
ture T, from that of a collection of "free" localized 4/ or 
5/ magnetic moments (subjected to the crystalline elec- 
tric field and to spin-orbit interaction) coupled weakly to 
the conduction electrons, to low-T local singlets where 
the localized moment is screened completely by the con- 
duction electrons by virtue of the Kondo effect. Sec. II. F. 

The energy gain of singlet formation Tk ~ 
Dexp{—l/NoJ) sets the temperature scale where this 
crossover occurs, i.e., the Kondo temperature Tk- Here 
Nq is the (unrenormalized) conduction-electron (c) den- 
sity of states at the Fermi level and J is the exchange 
constant between c and / electrons. The onset of a coher- 
ent state leading to a Fermi liquid at usually still lower 
temperatures is mostly signaled by a maximum in the 
electrical resistivity p{T). Although there exists a lot of 
discussion on the onset of coherence at a coherence tem- 
perature Tcoh thus defined, its relation to Tk remains 
under debate, see Sec. II. F. 2. Here we use for simplicity 
Tk as a measure of the energy gain due to singlet forma- 
tion, keeping in mind that Tk in HFS might be modified 
by interactions with respect to the single-ion Kondo tem- 
perature in dilute magnetic alloys. Experimentally, this 
Tk is often determined through fits of the intermediate- 
temperature data for C(T) or x(^) to single-impurity 
results. (Frequently the symbol T* is used for the char- 
acteristic temperature of a Kondo lattice, instead of Tk.) 

At sufficiently low T <C Tk, Fermi-liquid (FL) proper- 
ties are observed in many HFS with a very large effec- 
tive mass m* derived from the huge linear specific-heat 
coefficient ^ = C/T and a correspondingly large Pauli 
susceptibility, both being only weakly dependent on T. 
The electrical resistivity of a FL exhibits a contribution 
Ap = AT^ arising from electron-electron collisions, see 



Sec. II. E. 3. The phenomenological correlations 7 ^ x 
(Fisk, Ott, and Acppli, 1987) and A ^ 7^ (Kadowaki and 
Woods 1986, see Sec. II. E. 3) observed approximately for 
different HFS, do suggest the validity of the FL descrip- 
tion. The Wilson ratio R = (x/7)(7r^A:^//io/ioff) deviates 
from the free-electron value i? = 1 (Fisk, Ott, and Aep- 
pli, 1987). The observed values of i? ^ 2 to 5 can be 
accounted for in the frame of FL theory by a negative 
Landau parameter Fq of the order of —0.5 to —0.8, see 
Sees. II.B and II.F.2, Eq. (44). The problem of prov- 
ing or disproving FL behavior in HFS lies in the low 
energy scale set by the Kondo temperature Tk ^ 10 to 
100 K compared to conventional metals where Tp ^ 10^ 
to 10^ K. The first heavy-fermion system discovered was 
CeAls, with a 7 coefficient of ~ 1.5J/molK^ (Andres, 
Graebner, and Ott, 1975). The interest in these mate- 
rials increased tremendously when Stcglich et al. (1979) 
reported superconductivity in CeCu2Si2, with the heavy 
quasiparticles being responsible for the superconductiv- 
ity. The early work on heavy-fermion systems has been 
reviewed by Grewe and Steglich (1991). 

The competition between on-site Kondo interaction, 
quenching the localized magnetic moments, and inter- 
site RKKY interaction between these moments allows for 
non-magnetic or magnetically ordered ground states in 
HFS. In the Doniach picture (Doniach, 1977), this com- 
petition is governed by a single parameter, namely the 
effective exchange constant J between conduction elec- 
trons and local moments, which enters the characteristic 
energy scales Tk and Trkky ~ J^Nq) for Kondo and 
RKKY interactions, respectively (see Sec. II.F.2). The 
strength of the exchange interaction is usually tuned by 
composition or pressure. In addition, a large magnetic 
field can suppress Kondo screening. Owing to the ex- 
tremely strong dependence of the Kondo energy scale on 
the interatomic distance d, which arises from the expo- 
nential dependence of Tk on J, volume changes are of- 
ten the dominant effect in producing the magnetic-non- 
magnetic transition if isoelectronic constituents are sub- 
stituted against each other. 

Many HFS show some kind of static magnetic order 
with often very tiny ordered moments (of the order of or 
smaller than IO^^/ib)- A particularly intriguing exam- 
ple is CeAla, long considered the archetypal HFS with- 
out magnetic order. However, CeAla has been shown 
to order antiferromagnetically when produced in single- 
crystalline form (Lapertot et al., 1993). This is probably 
due to subtle structural differences, e.g., strains, between 
polycrystals and single crystals. An inhomogeneous dis- 
tribution of Kondo temperatures was inferred from NMR 
measurements in this material (Gavilano, Hunziker, Ott, 
1995). 

This example shows that disorder can have a deci- 
sive influence on the low-temperature properties even 
in nominally stoichiometric HFS samples. For non- 
stoichiometric substitutional alloys one can distinguish 

between substitution of the 4/ or 5/ site and on the lig- 
and site. In both pronounced effect may be an- 
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ticipated. In the former case, replacement of Ce or U by 
a non-magnetic atom in a otherwise stoichiometric HFS 
might lead to the formation of a Kondo hole. In the latter 
case, replacement of the ligand atom around a given Ce 
or U site might change the hybridization and hence of the 
local Tk- In general, predictions are difficult as to which 
effect will be stronger in a given system. Although much 
work on QPT that were tuned by composition of substi- 
tutional alloys has been done, stoichiometric compounds 
avoiding disorder are preferable, as many of the complica- 
tions, theoretically anticipated for samples with sizeable 
disorder (see Sec. III. J), will be absent. Alternatively, 
different tuning parameters should be employed to check 
the role of disorder, as has been done for CeCue-xAux- 



1. CeCue-xAux and CeCue-xAgx 

CeCue has been established as a HFS showing no long- 
range magnetic order down to the range of ~ 20 mK (Am- 
ato et aL, 1987; Onuki and Komatsubara, 1987). CcCug 
crystallizes in the orthorhombic Pnma structure and un- 
dergoes an orthorhombic-monoclinic distortion around 
200 K (Gratz et aL, 1987). The change of the orthorhom- 
bic angle is only small 1.5°). In order to avoid con- 
fusion, we always use the orthorhombic notation for the 
direction of the lattice vectors. CeCue exhibits a pro- 
nounced magnetic anisotropy with the magnetization ra- 
tios along the three axes Mc : Ma : Mt w 10 : 2 : 1 
at low T (Amato et aL, 1987). Schuberth et aL (1995) 
have extended the measurements of the specific heat C 
down to 10 mK and of the magnetic susceptibility % to 
even below 1 mK. Their analysis of x very low T (af- 
ter subtraction of an impurity contribution attributed to 
Gd) suggests magnetic order around 5 mK. This is backed 
by NQR measurements which likewise hint at (possibly 
nuclear) magnetic order (Pollack et aL, 1995). Direct ev- 
idence for magnetic order below 2 mK was found in the 
ac magnetic susceptibility and thermal expansion (Tsu- 
jii et aL, 2000). /iSR measurements have put an upper 
limit for a static moment of 10~^ to 10""^ /xs/Ce-atom 
(depending on the assumption of long-range magnetic vs. 
spin-glass order) above 40 mK (Amato et aL, 1993). 

Although CeCue does not order magnetically above 
5mK, the expectation C/T w const for a FL is not met 
very well (see Fig. 7). The single- ion Kondo model with 
Tk = 6.2 K does not fit the data below ^ 0.4 K (Schlagcr 
et aL, 1993). Instead C/T increases slightly towards low 
T which might be a precursor of the 5-mK order. On the 
other hand, the T^ dependence of the electrical resistiv- 
ity is rather well obeyed between 40 and 200 mK (Amato 
et aL, 1987). Antiferromagnetic fluctuations were ob- 
served in inelastic neutron scattering by peaks in the dy- 
namic structure factor S'(Q,w) for energy transfer Hlu = 
0.3 meV at Q = (1 0) and (0 1±0.15 0) (Rossat-Mignod 
et aL, 1988). The rather large widths of these peaks cor- 
respond to correlation lengths extending roughly only to 
the nearest Ce neighbors. These correlations vanish at 
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FIG. 7 Specific heat C of CeCue-xAua; in the vicinity of the 
QPT plotted as C/T vs. temperature T (logarithmic scale). 
Application of hydrostatic pressure at the respective critical 
value pc shifts C/T of the antiferromaguetie samples x — 0.2 
and 0.3 towards NFL behavior for x = 0.1 at ambient pres- 
sure. From v. Lohneysen et aL, 1996a, 1998a. The inset 
shows the Neel temperature Tn of CeCue-xAux vs. Au con- 
centration X as determined from specific heat (triangles) and 
magnetic susceptibility (circles). From Pietrus et aL, 1995. 



a field of ~ 2T applied along the easy c direction, also 
associated with a shallow maximum at 2 T in the differ- 
ential magnetic susceptibility dM/dB at very low T (v. 
Lohneysen, Schlager, and Schroder, 1993). This maxi- 
mum has been identified with the "mctamagnctic transi- 
tion" in loose analogy to the metamagnetic transition in 
strongly anisotropic antiferromagnets. 

Upon alloying with Au the CeCug lattice expands 
while retaining - in fact: stabilizing - the orthorhombic 
(at room temperature) Pnma structure (Pietrus et aL, 
1995). Thus the hybridization between Ce 4/ electrons 
and conduction electrons, and hence J, decrease, leading 
to a stabilization of localized magnetic moments which 
can now interact via the RKKY interaction. The result is 
incommensurate antiferromagnetic order in CeCug-xAux 
beyond a threshold concentration Xc ~ 0.1. This was first 
inferred from sharp maxima in the specific heat C(T) 
and magnetization M{T) (Germann et aL, 1988) and 
confirmed by neutron scattering (Chattopadhyay et aL, 
1990; Schroder et aL, 1994; v. Lohneysen et aL, 1998b). 

For 0.1 < a; < 1 where Au exclusively occupies the 
Cu(2) position in the CeCue structure, the Neel temper- 
ature Tn varies linearly with x (Fig. 7). For the stoichio- 
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FIG. 8 Temperature T1/2 and staggered moment /u for 
CeCu6-a;Aua:. Left: Concentration dependence of the tem- 
perature Txji where the magnetic entropy reaches 0.5R£n2, 
obtained from a temperature integration of the specific-hcat 
data. The decrease in T1/2 with x for x < 0.5 might reflect the 
decrease in Tk inferred an analysis of the high-field specific 
heat in the frame of a single-ion Kondo model (v. Lohneysen 
et al., 1996a), while the larger T1/2 at a; = 1 likely results from 
the strong AF order. Right: x dependence of the staggered 
moment per Co ion extracted from elastic neutron scattering. 
From V. Lohneysen et al, 1998b. Lines are to guide the eye. 



metric compound CcCusAu a complex magnetic phase 
diagram has been mapped out (Paschke et ai, 1994b). 
Beyond x = 1, In decreases again. This is due to a sub- 
tle change within the orthorhombic structure: for x < 1 
the lattice parameters a and c increase while b decreases 
with growing Au content, whereas for a; > 1 all a, b, 
and c increase (Pietrus et ai, 1995). The orthorhombic- 
monoclinic transition is quickly suppressed with increas- 
ing X and vanishes around a; = 0.14 (Grube et al., 1999). 
The order parameter, i.e., the staggered magnetic mo- 
ment per Ce atom as extracted from elastic neutron scat- 
tering data (v. Lohneysen et ai, 1998b), is shown in 
Fig. 8. The pronounced rise between x = 0.3 and 0.5 
may be related to the change of the magnetic ordering 
wavevector occurring in the same x range. 

Fig. 7 shows specific-heat data for concentrations in 
the vicinity of the critical concentration Xc « 0.1 plotted 
as C/T. The Neel temperature Tn manifests itself as a 
sharp kink in C/T which becomes less pronounced as Tn 
decreases. It is, however, still clearly visible for x = 0.15 
where Tn = 0.080 K as confirmed by a maximum in the 
susceptibility. For x = 0.1 we observe the non- Fermi- 
liquid behavior C/T = a£n(To/T), with a = 0.578 J/mol 
and Tq = 6.2 K, between 0.06 K and ~ 2.5 K, i.e., 
over almost two decades in temperature, rendering it one 
of the best examples of a logarithmic divergence of C/T 
for T ^ 0. (The positive deviations above 2.5 K can be 
attributed to phonon and crystal-field contributions to 
C.) Concerning the NFL behavior at Xc it is important 
to verify that it does not arise from some inhomogene- 
ity of the alloys, i.e., a distribution of magnetic ordering 
temperatures. A muon spin relaxation (/zSR) study has 
shown that there is no ordered magnetic moment in a x 
= 0.1 sample, with the detection limit /z < 1Q~^ ^s/Ce- 



atom (Amato et al, 1995). Likewise, a wide distribution 
of Kondo temperatures appears to be unlikely (Bernal et 
al., 1996). 

The specific-heat data can be temperature-integrated 
to obtain the entropy, where the dominant low- 
temperature contribution arises from the Ce moments. 
As a free moments (with effective spin 1/2) would give 
an entropy of EIn2, it is convenient to define a temper- 
ature, Ti/2, where the entropy reaches 0.5R£n2 - this 
defines a scale which is a measure of the (Kondo) tem- 
perature, below which the moments are screened. The x 
dependence of Tx/2 for CeCug-xArij; is shown in Fig. 8. 
Importantly, T1/2 does not vanish at Xg. Its decrease 
with X can be essentially understood from the weakening 
of the Kondo effect because of the lattice expansion upon 
alloying Au; the increase of T1/2 towards x = 1 is likely 
due to the quenching of the moments by the ordered an- 
tiferromagnetism. 

As expected from the correlation between molar vol- 
ume and Tn discussed above, Tn of CeCug-xAu^ de- 
creases under hydrostatic pressure p (Bogenberger and 
V. Lohneysen, 1995; Germann and v. Lohneysen, 1989; 
Sieck et at, 1997). The Neel temperature (again as deter- 
mined from the inflection point of C(T) above the max- 
imum) decreases linearly with increasing p for x = 0.3. 
For X = 0.2 a linear T^{p) decrease is also compatible 
with the data. Tn ~ is reached for 7 - 8 kbar and 3.2 
- 4 kbar for x = 0.3 and 0.2, respectively. At these pres- 
sures both alloys exhibit NFL behavior with, surprisingly, 
the same coefficients a and Tq for both, coinciding with 
the NFL alloy x = 0.1 and p = (Fig. 7). At 6.9 kbar 
for X = 0.2, the clear suppression of the low-T increase 
of C/T towards the data for CeCue indicates restora- 
tion of the FL (not shown), i.e., one can pressure-tune 
CeCu5.8Auo.2 all the way from antiferromagnetic order 
through the QCP to FL (v. Lohneysen, 1996b). Like- 
wise, application oi p = 6.0 kbar for x = 0.1 drives this 
alloy even further towards FL behavior (Fig. 7), as C/T 
at low T now falls even below the p = data of CeCue- 

We point out one peculiar feature of the C/T data: 
Both for the concentration- and pressure-driven transi- 
tions, C/T at e.g. 100 mK continues to increase - as a 
function of x or p - beyond the QCP when moving from 
the disordered to the ordered phase. Fig. 7. This behav- 
ior parallels that of the entropy. Fig. 8, discussed above. 

The magnetization M(T) measured in a field B = 

0. 1 T for X ~ 0.1 exhibits a cusp for T ^ which 
can be modeled as x ~ M/B — xo — ots/T between 
80 mK and 3 K. Roughly the same T dependence of M/ B 
is found upon reduction of the field to 0.01 T with a 
slightly stronger upturn towards low T below 0.3 K (v. 
Lohneysen et ai, 1998a). The 0.1 T data can be very 
well described also by a different functional dependence, 

1. e., x(T)-i -x(0)-i oc T« with a = 0.8 (Schroder et al., 
1998), see below. 

Fig. 9 shows p{T) for different CeCug-^Aua; alloys 
for current parallel to the orthorhombic o,-axis. For 
X < Xc f« 0.1, p{T) increases at low T as p{T) = po + AT"^ 
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FIG. 9 Electrical resistivity p of CcCufi-j. Au^, vs. tempera- 
ture T, with current applied to the a direction. Arrows indi- 
cate the Neel temperature. Inset shows data for a; = 0.1 along 
the b direction. For all directions, p = po + A'T is observed. 
From V. Lohneysen et al, 1998a. 



which is expected for a FL with dominant quasiparticlc- 
quasiparticle scattering for T ^ 0. This has been ob- 
served before for CeCue (Amato et al., 1987). For a; = 0.1 
a linear T dependence of p is observed between 20 mK 
and 0.6 K (sec Fig. 9), signahng NFL behavior. The 
anisotropic p{T) dependence of the magnetically ordered 
alloys can be qualitatively interpreted in terms of the ob- 
served magnetic order: p{T) for all alloys except x — 1 
increases below Tn for current directions with a non-zero 
projection of the magnetic ordering vector Q determined 
from clastic neutron scattering (v. Lohneysen et at, 
1998b). An increase of p{T) below Tn has been observed 
before in other HFS, for example, in CeRu2_xRha;Si2 as 
will be discussed below (Miyako et al, 1997). 

The abundance of low-energy magnetic excitations as 
Tn — > has been suggested to cause the NFL behav- 
ior at the magnetic instability (v. Lohneysen et al., 
1994). This is supported by the recovery of FL behav- 
ior in high magnetic fields B (Finsterbusch et al., 1996; 
V. Lohneysen et al., 1994). A negative deviation from 
the C/T r-^ £n{To/T) divergence is seen for all fields 
B > 0.2 T, with a crossover temperature roughly obeying 
B. A similar systematic recovery of FL behavior 
of a quantum critical system upon application of a mag- 
netic field has been observed in many other systems. We 
add that the high-field specific heat of all CeCu6-a:Aua; 



alloys including a; = 0.1 can be reasonably well described 
(Schlager et al., 1993; v. Lohneysen et al., 1996a) within 
a single-ion Kondo model. 

The £n{To/T) dependence of C/T and the linear T 
dependence of p in CeCue-ajAua; at the magnetic insta- 
bility have constituted a major puzzle ever since they 
were first reported. The LGW theories for 3d itinerant 
fermion systems predict C/T = jQ-fSVT and Ap ~ T^/^ 
for antiferromagnets (z ~ 2), while C/T = £n{To/T) and 
Ap ~ T^/^ are expected for ferromagnets {z = 3), see 
Sees. IILC-IILF. In addition, Tn should depend on the 
control parameter = x — Xc or rp = p^Pc as Tn ^ |r|''^ 
with V = z/{d + z-2) = z/iz + 1), Eq. (91), for d = 3, 
while for CeCu6-a;Auj; '0 = 1 for both Vx and Vp. Rosch 
et al. (1997) showed in an analysis similar in spirit to 
that of Millis (1993) that 2d critical fluctuations coupled 
to quasiparticles with 3d dynamics lead to the observed 
behavior C/T ~ £n{To/T), Ap ~ T and Tn ~ |r|. 

Let us discuss the question of 2d vs. 3d magnetism 
in CeCug-ajAuj;. CeCug-jjAuj; does exhibit 3d antifcr- 
romagnetic ordering, and the anisotropy of the electri- 
cal resistivity along different crystallographic directions 
does not exceed a factor of 2. Therefore CeCuc-xAu-j. 
looks like a 3d antiferromagnetic metal. The magnetic 
structure of CcCug-ajAux (0.15 < x < 1) has been in- 
vestigated with elastic neutron scattering (Okumura et 
al., 1998; v. Lohneysen et al, 1998b). An example of 
resolution-limited magnetic Bragg reflections is shown in 
Fig. 10. The magnetic ordering vector is Q = (0.625 
0.253) for X = 0.2 and remains almost constant up to x 
— 0.4. For larger x it jumps onto the a* axis, Q = (0.56 
0) for X = 0.5 and (0.59 0) for a; = 1. 

A detailed investigation of the critical fluctuations at 
Xc = 0.1 using inelastic neutron scattering (Stockert 
et al., 1998) showed that the critical fluctuations are 
strongly anisotropic and extend into the a*c* plane. This 
is inferred from a large number of / scans in the a*c* 
plane, some of which are shown in Fig. 11. Hence the 
dynamical structure factor S{q,fuj = 0.15meV) has the 
form of rods (see Fig. 10). Since a quasi-ld feature in re- 
ciprocal space corresponds to quasi-2d fluctuations in real 
space, the 2d LGW scenario (Rosch et at, 1997) appears 
to be applicable. The width of ^(q, iv) perpendicular to 
the rods is roughly a factor of five smaller than along 
the rods. It is an issue of current debate whether this 
anisotropy of the correlation length is enough to qual- 
ify the fluctuations as being 2d. The 3d ordering peaks 
for X = 0.2 and 0.3 fall on the rods for a; = 0.1 which 
therefore can be viewed as a precursor to 3d ordering 
(Fig. 10). Fig. 11 demonstrates the essentially similar, 
albeit broader S{q^,hco = const) dependence for samples 
away from the critical concentration, i.e., for a; = and 
0.2 (v. Lohneysen et al., 2002). 

The dynamic structure factor S{c[ = const, fkv) of 
CeCu6-a:Aux was investigated around Q = (0.8 0), i.e., 
on the rods (Fig. 10), by Schroder et al. (1998). They 
found a scaling of the dynamical susceptibility of the form 

X-i(q, E, T) = c-i [/(q) + {-lE + aT)«] (149) 
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FIG. 10 Neutron scattering results for CeCue-xAua,. Up- 
per panel: Resolution-limited magnetic Bragg reflections for 
CeCuB.8Auo.2 corresponding to an incommensurate magnetic 
ordering wavevector Q = (0.625 0.275). Prom v. Lohneysen 
et ai, 1998b. Lower panel: Position of the dynamic cor- 
relations {x — 0.1, hu = 0.1 mcV, T < 100 mK) and mag- 
netic Bragg peaks (0.2 < j: < 1.0) in the a*c* plane in 
CeCue-xAux. Closed symbols for x = 0.2 represents short- 
range order peaks. The vertical and horizontal bars indicate 
the Lorentzian linewidths for x = 0.1. The four shaded rods 
are related by the orthorhombic symmetry (we ignore the 
small monoclinic distortion). The inset shows a schematic 
projection of the CeCue-xAux structure onto the ac plane 
where only the Ce atoms axe shown. The rods in reciprocal 
space correspond to planes in real space spanned by b and 
lines in the inset. Prom Stockert et al., 1998. 



with an anomalous scaling exponent a = 0.74 (a 
Lorentzian fluctuation spectrum would be described by 
a = 1). This translates to 

X"{E,T) = T-'^giE/ksT), (150) 
g{y) = csin[atan-i(y)]/(y2 + l)«/2. (151) 

Fig. 12 shows the scaling obtained by Schroder et al. 
(1998), which is confirmed by more recent data (Schroder 
et al., 2000). The exponent a ~ 0.8 implies z = 2/a = 
2.5. It is interesting to note that Eq. (149) implies for 
the static uniform susceptibility (E = 0,q = 0) 



[x'{T)-'-x'{0)-']=c-\aTr 



(152) 



Schroder et al. (1998) showed that indeed the static imi- 
form susceptibility can be described by Eq. (152) with 



a = 0.8 between 0.08 and 8 K in agreement with a = 0.74 
from neutron scattering. The fact that this quantum 
critical scaling persists to temperatures well above Tk 
appears surprising. 

The two experiments by Stockert et al. (1998) and 
Schroder et al. (1998) are complementary in that the for- 
mer focuses on the q dependence and the latter on the 
E dependence of x(q, E) of the critical fluctuations in 
CeCue-xAux- The data agree qualitatively in the range 
where they overlap. An open question is the origin of the 
fluctuation spectrum being effectively two-dimensional, 
and whether this fact may be generic for certain Ce-based 
materials. It is clear that the 2d planes are spanned by 
the b axis and the connection between next-nearest neigh- 
bor Ce atoms (sec Fig. 10). Only more detailed investi- 
gations can establish if, perhaps, the low-dimensionality 
arises from a strong spatial anisotropy or frustration of 
the RKKY interaction or the c / hybridization. 

The unusual type of scaling in the dynamical suscep- 
tibility of CeCug-ajAua. at the QPT is incompatible with 
the LGW model - there, E/T scaling as in Eq. (149) is 
only expected below the upper critical dimension (which 
is d = 2 for the metallic antiferromagnet). The ex- 
periments therefore prompted new theoretical concepts, 
based on the idea that at the QCP the (lattice) Kondo 
effect breaks down (Sec. 111.1). Specific proposals in- 
clude a so-called "local" quantum critical point (Si et 
al., 2001, 2003), spin-charge separation (Coleman et al., 
2001; Pepin, 2005), or fractionalization of the Fermi liq- 
uid (Senthil, Sachdev, and Vojta, 2003; Senthil, Vojta, 
and Sachdev, 2004). The local QCP scenario of Si et 
al. (2001) requires the existence of two-dimensional crit- 
ical spin fluctuations; as mentioned above, this is not 
inconsistent with the experiments. Among the currently 
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FIG. II Inelastic neutron scattering scans in the reciprocal 
a*c* plane of GeCue-xAux for (a) a: = with an energy 
transfer hw = 0.15 meV, (b) x = 0.1, hw = 0.10 meV, and 
(c) x = 0.2, hu; = 0.15 meV. Data were talcen at the triple 
axis spectrometers IN12 (x = 0) and IN14 {x = O.I; 0.2) at 
the ILL Grenoble. The individual scans are shifted by 100 
counts {x = 0; 0.2) for 150 counts (x = O.I) with respect to 
each other. Scans along (2.8 — ^0/) for x — 0.2 are symmetry 
equivalent to (1.2 ^ /) scans (a; = 0; 0.1). From v. Lohneysen 
et al, 2002. 
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FIG. 12 Scaling plot of inelastic neutron scattering data for 
CeCuB.8Auo.2 at g = (0.8 0) vs. E/k^T. Solid line corre- 
sponds to a fit of the scaling function Eq. (150) with a = 0.74. 
Inset shows how the quality of the scaling collapse varies with 
a. From Schroder et al, 1998. 



available scenarios, only the local QCP picture contains 
E/T scaling in the dynamical susceptibility (by virtue 
of the model at all wavevectors). Moreover, the frac- 
tional exponent obtained from numerical calculations of 
the EDMFT equations (Grempcl and Si, 2003), in the 
experimentally relevant situation of Ising symmetry, is 
a Ri 0.72, i.e., close to the experimental value of ~ 0.74. 

As discussed in Sec. III. I, the breakdown of the Kondo 
effect should be connected with an additional energy scale 
(associated with Fermi surface fluctuations) vanishing at 
the QCP, apart from that of magnetic ordering. Experi- 
mentally, the situation is ambiguous. In CeCug-xAuj, es- 
sentially all experimental results indicate that some scale 
of the order of a few Kelvin stays finite at the QCP for 
a doping of x = 0.1. This can be seen from the ex- 
istence of a maximum in the resistivity at Tm sa 4K 
(for current fiow along the a direction) (v. Lohneysen 
et al., 1998a, 2002), or the temperature T1/2, where the 
entropy reaches 0.5R£n2 (Fig. 8). (For YbRh2Si2 simi- 
lar conclusions on T1/2 can be drawn from the specific- 
heat data. However, for this material recent magneto- 
striction measurements have indicated the existence of 
thermodynamic low-energy scales which vanish at the 
QCP (Gegenwart et al, 2007).) In CeCue-xAu^, the 
specific-heat coefficient is large even in the ordered phase: 
As is apparent from Fig. 7 and already noted above, C /T 
below 100 mK is larger for the samples just above the 
critical concentration than for Xc, before decreasing for 
concentrations approaching x — 1 (see also Mock et al. 
1994). 

Despite these open questions it should be stressed that 
CeCue-xAua; is one of the best characterized heavy- 
fermion systems exhibiting NFL behavior. It is reward- 
ing that the unusual behavior of the thermodynamic and 
transport quantities at the QCP can be traced back in a 
consistent fashion to an unusual low-dimensional fluctu- 
ation spectrum. 

Recent measurements of the ^^Cu nuclear spin-lattice 
relaxation time Ti on one of the five inequivalent Cu sites 
in CeCus.gAuo.i down to O.IK, employing the 6.25 MHz 



Cu NQR line, revealed a non-exponential nuclear mag- 
netization recovery which was attributed to different Cu- 
neighbor configurations in the alloy (Walstedt et al, 
2003). More importantly, Tj~^ was found to vary as T° ''^ 
at low T. Since Tj~^ is proportional to the weighted 
squared average (A(q)^) of transverse hyperfine cou- 
plings over the Brillouin zone, the "agreement" with the 
exponent a found in the T and E dependence of x{<l: ^) 
(149) is entirely accidental. 

Although CeCug alloyed with Au is the most thor- 
oughly studied system, other alloys derived from CeCue 
have been widely studied as well. Magnetic ordering in 
the system CeCu6-a;Aga, was discovered simultaneously 
with CeCue-xAuj, (Praunberger et al., 1989; Gangopad- 
hyay et al., 1988; Germann et al., 1988). CeCue-iAgj, 
displays critical behavior around Xc = 0.2, with thermo- 
dynamic properties being very similar to CeCus.gAuo.i, 
e.g., a logarithmic divergence oiC/T. Interestingly, ther- 
mal expansion measurements have uncovered that the 
critical behavior is incompatible with both the 3d and 
2d LGW model: the Griineisen parameter P, Eq. (70), 
is expected to diverge in the quantum critical regime for 
T ^ as F ~ T-i/(^'') with zu = I. However, the ex- 
periment shows a much weaker (logarithmic) divergence 
of F (Kiichler et al, 2004). 

For field-tuned transitions the behavior appears 
to be different: NFL behavior in a poly crystalline 
CeCu4.8Agi.2 alloy, subjected to a magnetic field, was 
reported, i.e., approximately C/T ^ iniTo/T) between 
0.35 and 2.5 K (Heuser et al, 1998a). However, in the 
light of the strong magnetic anisotropy of this system, 
with B affecting crystallites of different orientations quite 
differently, this result should be viewed with caution. 
Subsequently, Scheidt et al. (1998) reported specific-heat 
data down to 0.07 K on a CeCu5.2Ago.8 single crystal 
with Tn = 0.7K. For a critical magnetic field Bc = 2.3T 
applied along the easy direction, C'/T varies logarithmi- 
cally from ~ 1.5 K down to 0.2 K and then levels off. 
Furthermore, the resistivity p{T) of CeCu5.2Ago.8 can be 
described by a T^-^ dependence at Be- The authors in- 
terpret the data within the LGW scenario of Moriya and 
Takimoto (1995), Sec. III.E, with = 3 and = 2. 

A detailed comparison of pressure-tuned and field- 
tuned QCP on the same system CeCu5.8Auo.2 (v. 
Lohneysen et al, 2001) demonstrated that field, as op- 
posed to pressure, drives the system to a 3d LGW 
quantum critical point. The data p{T) and C/T for 
CeCu5.2Ago.8 for Be look quite similar to those for 
CeCu5.8Auo.2 for B = 0.3 or 0.5 T, i.e., in the region 
of Be ~ 0.4 T as determined from neutron scattering. 
This difference to concentration and - very likely also 
to pressure tuning was recently directly confirmed by 
inelastic neutron scattering where the observed co/T^'^ 
scaling corresponds to the LGW scenario (Stockert et al, 
2006b). It is open whether the apparent similarity to the 
spin-fluctuation theory represents cross-over phenomena 
or the approach to a new QCP. 

A few experiments have been performed on fur- 
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ther CeCug-dcrivcd alloys. For CcCuG-xPda; and 
CeCug_2:Pta; the NFL behavior was found in the spe- 
cific heat at the critical concentration for the onset of 
magnetic order for Pd and «0.1 for Pt), while 

for CeCue-ajNia; a non-magnetic ground state is stabi- 
lized for smaU x (Sieck et al., 1996). Notably, the C/T 
data for the Pd and Pt alloys at their respective Xc agree 
well with those for CeCue-a^Aua: at the quantum critical 
point, suggesting universality for all QPT in CeCue al- 
loys in zero magnetic field. Note, however, that recent 
data on a CeCug-xAgj, polycrystal with x = 0.2, be- 
lieved to be close to the quantum-critical concentration, 
exhibit a somewhat steeper slope of C/T vs. £n{TQ/T), 
with, however, a gentle leveling-off towards low T, lead- 
ing in effect to a C/T larger by 20% at 0.1 K compared 
to CeCus.gAuo.i (Kiichler et al, 2004). 
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FIG. 13 Magnetic contribution to the specific heat divided 
by T, Cm/T, vs. temperature T (logarithmic scale) for 
Cei-a;Laa;Ru2Si2 for x = 0,0.05 and 0.075. Solid lines are 
fits to the LGW model. Prom Kambe et al., 1996. 



2. Cei-a:Laa;Ru2Si2 and Ce(Rhi_a;Rua:)2Si2 

CeRu2Si2 is, like CeCug, a canonical heavy-fermion 
system without apparent magnetic order. It crystallizes 
in the tetragonal I/4mmm structure. The gromid state 
doublet of Ce'^+ in CeRu2Si2 exhibits Ising character 
(Regnault et al., 1988). Static magnetism with very small 
ordered moment was detected by muon-spin-rotation ex- 
periments (Amato et al., 1994). The magnetic order de- 
velops below about 2 K and the ordered moment reaches 
a value of ^ 10"'^ /is/Ce-atom at low temperatures (T < 
0.1 K). The Kondo coupling between Ce 4/ electrons and 
conduction electrons is stronger than in CeCue, this is 
also suggested by the rather low (for a heavy-fermion sys- 
tem) 7 value of 0.36 J/molK^ (Fisher et al., 1991). From 
specific-heat data Tk is found to be 20 K. Magnetic in- 
tersite correlations are also much stronger in CeRu2Si2 
than in CeCue- This has been shown by inelastic neutron 
scattering where g-dependent correlations are peaked at 
ki = (0.3 0) and k2 = (0.3 0.3 0) (Rossat-Mignod et 
al., 1988). In a recent inelastic neutron scattering study, 
the spin fluctuations were interpreted within the LGW 
model (Kadowaki, Sato, and Kawarazaki, 2004). 

A metamagnetic transition occurs in a field of Bm ~ 
8T (Hacn et at, 1987; Holtmeier et al, 1995), which 
manifests itself as a steep jump in the magnetization 
curve M{B). The corresponding maximum in dM/dB at 
Bm = 8T is much more pronounced and much sharper 
than in CeCue at Bm = 2 T. It is interesting to note in 
the present context that all the way through the meta- 
magnetic transition, CeRu2Si2 exhibits well-defined FL 
behavior as evidenced by the prevalence of de Haas-van 
Alphen (dHvA) oscillations very close to the transition 
(Tautz et al, 1995). Some of the branches could be de- 
tected to within 0.005 T of Bm = 7.8 T. However, the 
branch of the heavy quasiparticles is strongly affected by 
the metamagnetic transition; its effective mass changes 
from m* = 12 mo at 6 T to 4 mo at 12.5 T. 

Alloying CeRu2Si2 on the Ce site, i.e., diluting with 
non-magnetic La, leads to long-range antiferromagnetic 



order in Cei_j,Laj,Ru2Si2. This might at first sight be 

even more siirprising than the observation of antiferro- 
magnetism induced by Au-doping of CeCug . It may, how- 
ever, be simply attributed to a lattice expansion caused 
by replacing Ce by La (the famous lanthanide contraction 
backwards), leading to a strong reduction of Tk which 
outweighs the reduction of the number of Ce 4/ mo- 
ments for small x. This is nicely supported by measure- 
ments on Cei_2;(Lao.63Yo.37)a;Ru2Si2 which is isochoric 
to CcRu2Si2 and does not show magnetic order (Sakak- 
ibara et al, 1993). 

The long-range antiferromagnetic order for x > 0.1 
manifests itself by maxima in specific heat and magnetic 
susceptibility (Fisher et al., 1991) and in an increase of 
the electrical resistivity below the ordering temperature 
(Djerbi et al., 1988). Applying hydrostatic pressure de- 
creases Tn as shown for x = 0.2 where Tn for pc ~ 
7kbar. At the same time, the ordered moment as derived 
from elastic neutron scattering vanishes (Regnault et al, 
1990). However, the pressure dependence Tn ~ \p—pc\'^ 
is clearly sub-linear. This matches the sub- linear T^{x) 
dependence where Tn vanishes for a critical concentra- 
tion Xc « 0.075 (Quezel et al, 1988). 

Upon approaching Xc ~ 0.075 where Tn 0, the 
specific-heat coefficient increases towards low T (Flou- 
quet et al., 1995). Kambe et al. (1996) performed ad- 
ditional experiments for x = 0.075 which clearly show 
the strong enhancement of C /T . They describe the data 
consistently in terms of the LGW model as applied to 
heavy-fermion systems by Moriya and Takimoto (1995). 
In contrast to CeCue-xAua;, the data can be quantita- 
tively interpreted within the latter model for a 3d an- 
tiferromagnet, i.e., the specific-heat coefficient 7 crosses 
over from a logarithmic T dependence in an intermedi- 
ate T range to C/T = 70 - /3\/T for T ^ (Fig. 13), 
which is the asymptotic behavior within the LGW model. 
Sec. III.D. Likewise, the electrical resistivity at the criti- 
cal concentration shows the LGW dependence for a disor- 
dered antiferromagnet, Ap ~ T^/^ (Kambe et al, 1996). 
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Inelastic neutron-scattering experiments on 
Ceo.925Lao.o75Ru2Si2 showed three-dimensional cor- 
relations for this system (Knafo et ai, 2004; Raymond 
et al., 1997). At the wavevector Qi = (0.69 1) which 
corresponds to antiferromagnetic order, an approximate 
scaling of the form 

X"{E, T) = T-'g{E/{kBTf^) (153) 

with g{y) = y/{l + y)"^ was observed, however, only for 
T > 5K. As discussed by Knafo et al. (2004), this ap- 
parent scaling is probably related to a crossover regime, 
see Sec. I.C, as the Kondo temperature in this system is 
about 17 K. Furthermore, the linewidth of the dynamical 
fluctuations does not vanish for T ^ as it should for a 
quantum critical point. This may indicate that x — 0.075 
is a little off the QCP. Nevertheless, the linewidth is re- 
duced by a factor of 5 compared to the pure compound 
while the correlation length for T ^ is only 1.5 times 
larger than for pure CeRu2Si2. The dynamic suscepti- 
bility x(q, ^^) is interpreted in terms of the RPA, which 
yields an overall satisfactory agreement between the pa- 
rameters of the LGW model as derived from inelastic 
neutron scattering and those derived from the specific 
heat (Raymond et al., 1997). As a final point we men- 
tion that in addition to the inelastic signal, Raymond et 
al. (1997) observed below 1.8 K a very small elastic signal 
at Q = (0.69 1 0) for x = 0.075. This is surprising be- 
cause this sample is believed to be just below the critical 
concentration Xc = 0.08 for the onset of magnetic order. 
The linewidth is only slightly larger than the instrumen- 
tal resolution, and the ordered moment is ~ 0.02 /xb at 
the lowest temperature which may be due to concentra- 
tion fluctuations. The evolution of spin dynamics was 
also studied for an ordered sample with x = 0.13 around 
the critical pressure pc = 2.6kbar (Raymond et al., 2001). 

To summarize, Cci_2;Laa.Ru2Si2 appears to be a sys- 
tem with a QCP exhibiting 3d critical fluctuations. 
Hence the LGW model describes the thermodynamic 
transport and neutron-scattering data rather well. 

CeRu2Ge2 exhibits magnetic order, due to its larger 
Ce-Ce separation as compared to CeRu2Si2. Unlike 
many other Ce compounds, it orders ferromagnetically 
below Tc = 8K (Bohm et al, 1988; Pontes et at, 1996). 
Upon applying pressure, the ferromagnetism is quickly 
suppressed, and a rather complex magnetic phase dia- 
gram with several phases evolves. By ac susceptibility 
measurements under pressure, and in analogy with dop- 
ing studies of CeRu2(Gei_a;Sia;)2 where neutron scatter- 
ing has been done, these phases are identified as anti- 
ferromagnetic. Magnetic order in CeRu2Ge2 disappears 
at Pc ~ 6.5 GPa; at Pc NFL behavior is observed in the 
electrical resistivity, p = po + A'T (Siillow et al., 1999; 
Wilhelm et al., 1999). The system provides a nice exam- 
ple of the equivalence between hydrostatic pressure and 
chemical pressure exerted by smaller Si atoms. 

CeRh2Si2, which is isostructural to CeRu2Si2, is a 
local-moment antiferromagnet (Tn = 36 K) with a rel- 
atively large ordered moment oi n = LS/Us/Ce-atom 



(Quezel et al., 1984). Tn can be driven to zero by hy- 
drostatic pressure oi Pc = 9kbar (Thompson, Parks, and 
Borges, 1986). T^{p) varies slowly first and then precipi- 
tously drops to zero upon approaching Pc (Movshovich et 
al., 1996). This may be suggestive of a first-order tran- 
sition near T = 0, perhaps even stronger than the one 
observed in MnSi (Pfleiderer et al., 1997), Sec. IV.B.2. 
Clearly, a first-order transition would cut off quantum 
fiuctuations. Consequently, anomalies in the specific 
heat associated with NFL behavior were not observed 
in CeRh2Si2 around Pc (Graf et al., 1997). In addition, 
the resistivity above pc exhibits a PL-like behavior 
(Movshovich et al., 1996). It is interesting to note that in 
the vicinity of pc and beyond, CeRh2Si2 becomes super- 
conducting with Tc between 0.2 and 0.4 K (Movshovich 
et al, 1996), see also Sec. IV. C. 

The magnetic phase diagram of Ce(Rui_a;Rha;)2Si2 has 
been studied by several groups (Calemczuk et al., 1990; 
Kawarazaki et at, 1995; Lloret et al., 1987) with a num- 
ber of techniques, including neutron scattering. The data 
reveal a complex T^{x) phase diagram with two differ- 
ent antiferromagnetic phases, one between x « 0.05 and 
0.27, with a maximum Tn of ~ 5K, and the other phase 
emanating from CeRh2Si2 {x = 1), with Tn decreasing 
steeply to ~ 13 K for x = 0.95 and then exhibiting a 
plateau at Tn = 12K imtil x = 0.6, when it drops again 
and reaches Tn = around Xc ~ 0.5. The magnetic 
structure in both concentration ranges has been deter- 
mined by elastic neutron scattering. For the CeRh2Si2- 
derived phase, a wavevector Q = (5 5 0) corresponding 
to an antiferromagnetic coupling within (0 1) planes 
(Quezel et al, 1984) was reported to exist down x = 0.6 
(Lloret et al., 1987), with ordered moments 1.5^B/Ce- 
atom for x = 1 and 0.65/iB/Ce-atom for x ~ 0.65. This 
structure was essentially confirmed by Kawarazaki et al. 
(1995), although these authors observed for x = 1 addi- 
tional reflections corresponding to Q' = (^ ^ i) which 
already had been seen by Grier et al. (1984). It is not 
clear whether this behavior arises from a homogeneous 
single phase or a domain structure. 

A large effort has been devoted to analyze the low- a; 
magnetic phase centered around a; = 0.15. A sinusoidal 
spin-density wave with the moments ordered along the 
propagation axis has been inferred (Kawarazaki et al, 
1995; Miyako et al., 1996) which shifts from incommen- 
surate to commensurate: Q = (0 0.42), (0 0.45) 
and (0 0.5) for x = 0.15, 0.2, and 0.25, respectively. 
This magnetic structure is quite different from that of 
La-diluted CeRu2Si2 discussed above. Recent inelastic 
neutron scattering experiments for x = 0.03 revealed a 
y3/2 dependence of the linewidth at the antiferromag- 
netic wavevector Q = (0 0.35), compatible with the 
LGW scenario (Kadowaki et al., 2006). 

It should be mentioned that p{T) increases strongly be- 
low Tn which is interpreted as arising from a spin-density 
wave (Miyako et al., 1997; Murayama et al., 1997). A 
detailed scenario for Fermi-surface nesting has been de- 
veloped for X = 0.15 (Miyako et al., 1996). Deviations 
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from FL behavior have been observed in the region where 
Tn « for X 0.5 (Graf et ai, 1997) and x = 0.4 
(Taniguchi et ai, 1998) in both specific heat and mag- 
netic susceptibihty. While initial data for x = 0.4 show 
C/T ~ iniTo/T) between 10 K and 0.15 K, data taken 
by the same group with higher precision indicate a slight 
leveling off below 1 K (Gu et al, 2002). C/T for x = 0.5 
levels off more strongly below ~ 1 K. At least part of the 
data can be described by a Kondo disorder model with a 
distribution of Tk (Graf et al, 1997). 

Let us point out that the T^{x) functional dependence 
of Ce(Rui_3,Rh2,)2Si2 for x > 0.5 does not at all resemble 
the Tn(p) dependence of CeRh2Si2. This might serve as a 
warning in taking "chemical pressure" effects too literally 
and reducing changes brought about by alloying to sim- 
ple volume changes. Even for CeCu6-rrA\i,r where the 
functional dependence of Tt^{p) and Tj^{x) arc the same 
(both with an exponent tp = 1) as discussed above, and 
where isoelectronic substituents are used (Au vs. Cu) as 
opposed to the present case (Ru vs. Rh), the calculated 
change of Tn with volume from Tn(p) via the compress- 
ibility does not correspond to the change calculated from 
Tjq[x) with the x-ray determined volume V{x) (Germann 
and V. Lohneysen, 1989; Pietrus et al., 1995). The com- 
parison with the CeRu2(Sii_a;Gea;)2 system above might 
suggest that the chemical pressure exerted by the poly- 
valent metal is more akin to hydrostatic pressure than 
that of transition metals with partly covalent bonding. 



3. CeCu2Si2 and CeNi2Ge2 

GeCu2Si2 was the first heavy-fermion superconduc- 
tor, discovered nearly 30 years ago by Steglich et al. 
(1979). Like the CcRu2Si2 family, it exhibits the tetrag- 
onal I/4mmm structure. This is not the place to re- 
view the tremendous amount of work that has been done 
ever since on this compound, revealing a very complex 
ternary phase diagram in the vicinity of the 1:2:2 stoi- 
chiometric composition (StegHch et at, 1996, 1997). In 
a small region around this point, three different types 
of crystals exist with different ground states: A (mag- 
netically ordered ground state with superconducting mi- 
nority phase), S (superconducting ground state), and AS 
(both states almost degenerate). Which type of crystal 
is prepared depends on the stoichiometry, on the prepa- 
ration method (polycrystals versus single crystals), and 
heat treatment (Steglich et al., 1996). Unfortunately, 
the different types could up to now not be traced to 
different crystallographic and/or defect- induced proper- 
ties. Recently, the type of magnetic order in A-type 
crystals was identified by neutron diffraction (Stockert 
et al., 2004). The observed incommensurate magnetic 
order is compatible with a spin-density wave, with the 
ordering wavevector determined by the band structure 
of heavy fermions as suggested by renormalized band- 
structure calculations. These findings support the early 
suggestion that the A/S transition can be viewed as a 



QCP (Steglich et at, 1996, 1997), where a variety of NFL 
features compatible with the 3D LGW scenario are ob- 
served. 

The A phase in a polycrystallinc sample can be sup- 
pressed by a moderate pressure ~ 2kbar (Gegenwart 
et al, 1998) where the resistivity exhibits Ap ~ ^3/2. 
At 6.7kbar the specific heat can be consistently inter- 
preted within the LGW theory for an antiferromagnet. 
(Specific-heat data at 2kbar were not reported.) How- 
ever, the zero-pressure starting point does not seem to 
correspond to a FL: while Ap ~ above the temper- 
ature Ta = 0.75 K of the A-phase transition, 7 = C/T 
is strongly T-dcpendent. Another puzzle results when 
an S-typc single crystal lacking the A-phase signature is 
investigated (under zero pressure). In order to suppress 
superconductivity, a moderate field has to be applied. At 
T > 0.2 K and for 0.2 T < B < 6 T, Ap varies as T^/"^ 
and 7 almost as 70 — P\/T suggesting the proximity of 
the QCP, and corroborating the high-pressure data of the 
A-phasc crystal discussed above. However, when cooling 
below 0.2 K, 7 shows a very steep upturn towards low 
T that cannot be ascribed to nuclear Zeeman contribu- 
tions. In the same T range, Ap ~ T^/"^ is still observed. 
The latter is in keeping with spin-fluctuation theory for a 
rather disordered system (po ~ 30 /iilcm) as pointed out 
by Rosch (1999, 2000), cf. Sec. III.F. For a disordered 
system at the QCP, a recovery of FL behavior Ap ~ 
towards T — > is not expected (Rosch, 1999). Gegen- 
wart et al. (1998) suggest that the "disparate" behavior 
of transport and specific heat near the QCP in CeCu2Si2 
is due to a "break up" of the heavy quasiparticles. 

The compound CeCu2Gc2, isostructural to CeCu2Si2, 
orders antiferromagnetically with Tn = 4.15 K (de Boer 
et al., 1987) which can be attributed to the larger vol- 
ume and weaker c-/ hybridization, in keeping with the 
Doniach model. Upon application of a rather large pres- 
sure, superconductivity occurs for p « 70 - 80kbar with 
Tc « 0.6 K (Jaccard, Behnia, and Sierro, 1992). Tc re- 
mains independent of pressure up to 120kbar, reaches a 
maximum Tc = 2K at 160kbar and drops again (Var- 
goz and Jaccard, 1998). The magnetic transition tem- 
perature as determined from a kink of p{T) decreases 
only slightly with increasing pressure and appears to sat- 
urate around 2K at 94kbar; beyond this pressure a sig- 
nature in p{T) can no longer be observed (Jaccard et al., 
1999). The appearance of superconductivity near the 
pressure where magnetic order is suppressed is very in- 
teresting. As in CeCu2Si2 (Bruls et al., 1994; Stockert 
et al., 2006a), superconductivity and magnetic order (A 
phase) appear to be mutually exclusive. Because of the 
possible precipitous drop of Tn(p), this transition might 
be of first order, thus suppressing quantum fiuctuations 
and the concomitant NFL behavior expected only for a 
continuous magnetic non-magnetic transition. It was al- 
ready pointed out by Jaccard et al. (1999) that quantum 
critical behavior is not observed: the resistivity exponent 
a in Ap ^ T" has almost exactly the FL value 2 at the 
critical pressure. 
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CeNi2Ge2 is a stoichiometric HFS where pronounced 
NFL features are observed in thermodynamic properties 
(Gegcnwart et al., 1999; Grosche et al, 2000; Kiichler et 
ai, 2003): The resistivity exponent a, Ap ~ T", varies 
between 1.2 and 1.5 below 4K dependent on the residual 
resistivity po- This variation of a(po) may be understood 
in terms of the competition between magnetic scatter- 
ing and impurity scattering (Rosch, 2000), see Sec. III.F. 
While early specific-heat data between 0.35 K and 5 K 
were found to follow C/T oc —inT (Gegenwart et al., 

1999) , later measurements by Kiichler et al. (2003) could 
be described by C/T = 70 — PVt as expected for the 
LGW scenario, when corrected for an (unexplained) steep 
low-T upturn C ~ T~^. (For a discussion of the sam- 
ple dependence of the specific heat of this material, see 
Kiichler et al. 2003.) The divergence of the Griineisen 
parameter F ^ for T ^ is also in line with this 
model (Kiichler et al., 2003). It should be mentioned 
that in very pure CeNi2Ge2 samples traces of supercon- 
ductivity (Gegenwart et al, 1999), or even a complete 
resistive transition to superconductivity (Grosche et al., 

2000) were found. Thus, CeNi2Ge2 appears to be one 
of the few systems following the predictions of the LGW 
model for 3d antiferromagnets. 



4. CeTlns (T = Co, Rh, Ir) 

A few years ago, members of the series CcTIns (T 
= Ir, Rh, Co) were discovered as HFS with interesting 
ground states such as superconductivity at ambient pres- 
sure in CeColng with = 2.3 K (Petrovic et al, 2001b) 
as well as antiferromagnetic order in CeRhlns with Tn = 
3.8 K that under pressure gives way to superconductiv- 
ity (Hegger et al, 2000). The crystal structure of these 
compounds is derived from cubic Celns intercalated with 
nn2 layers along the (0 1) direction. This gives rise to a 
rather pronounced two-dimensionality in various physical 
properties. Alloys across this isoelectronic series display 
a rich variety of interplay between superconductivity and 
magnetic order, as can be seen in Fig. 28 below (Pagliuso 
et al., 2002). We will just mention a few findings about 
quantum criticality in these systems, focusing on the stoi- 
chiometric parent compounds. Space limitations force us 
to neglect the many interesting studies on alloy series, 
e.g., by Jeffries et al. (2005). 

Although CeCoIns does not order magnetically, in- 
tersite magnetic correlations seem to be rather strong: 
transport and thermodynamic data on La-diluted 
GeCoIns were interpreted in terms of a single-ion Tk = 
1.7K andTcoh = 45K > Tk (Nakatsuji et al, 2002). Pro- 
nounced NFL features are observed in CeColna above the 
upper critical field Bc2 = 4.95 T when superconductivity 
is suppressed (Bianchi et ai, 2003b). Fig. 14 shows C/T 
vs. £nT for 5 T < S < 9 T. The data can be reasonably 
well described by the LGW model where the parameter 
Do indicates the distance from the QCP. {yo is equivalent 
to r of Sec. III.) 
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FIG. 14 Magnetic specific- heat data of CeCoIns. (a) 7(T) = 
Cm{T)/T in magnetic fields B = ^oH \\ [001]. Dashed lines 
for 8 and 9T emphasize the FL behavior with constant 7. 
Left(right) arrows indicate the crossover to FL behavior for 
8T (9T). Solid hnes are fits to the LGW spin-fluctuation 
model for each field, with the corresponding values of j/o in- 
dicating the distance from the critical field Be- (b) Scaling 
analysis of the data in (a) for a = 0.71 and (3 = 2.5. Inset: 
Plots of Q {(3) which minimize for a given /3, and for 
these a and (5. From Bianchi et al., 2003b. 



At the same time, the data are shown to exhibit scaling 
in the form of 7(B,T) - 7(Bc,T) w (AS)"/(AS/T^) 
where AB — B — Be and Be ~ Bc2 is the (quantum-) 
critical field. Best scaling collapse is obtained for a = 
0.71 and /3 = 2.5. Note, however, that a value of a ^ is 
not consistent with the expected recovery of FL behavior 
for T — > and B > Be, which seems to be a more likely 
interpretation of the experiment. Scaling with a = has 
been previously obtained for U0.2Y0.8Pd3 with (3 = 1.3 
(Andraka and Tsvelik, 1991), for YbRh2Si2 with (3 = 1.05 
(Trovarelli et ai, 2000), for CeCua;Aga; with exponents 
/3 = 0.85, 1.35, 1.6, 1.7 for x = 0.09, 0.48, 0.8, 1.2 (Heuser 
et al., 1998a, b). At least in YbRh2Si2 and CeCu5.8Ago.2 
the scaling does not hold at the lowest T. In view of the 
different values of a and /?, the physical significance of 
the scaling is not clear. 

The electrical resistivity p{T) was shown by Bianchi 
et al. (2003b) to follow the LGW model as discussed by 
Moriya and Takimoto (1995), Sec. IILE, with a set of yo 
parameters consistent with the specific heat. The small 
value yo ^ 0.01 suggests that CeColns is indeed close to 
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a QCP. The existence of a QCP in CcCoIns very close 
to Bc2 was also inferred independently from resistivity 
and magneto-resistivity measurements by Paglione et al. 
(2003). Surprisingly, when 3^2 is suppressed upon alloy- 
ing CeCoIn5_a;Snj;, the field range where NFL behavior 
is observed tracks Bc2 (Bauer et al, 2005). However, 
recent resistivity measurements on CeCoIns under pres- 
sure have shown that Bc2 decreases much faster than the 
quantum critical field, indicating that the two phenom- 
ena are not related (Ronning et al, 2006). 

While application of hydrostatic pressure drives 
CeCoIns away from quantum criticality (Nicklas et 
al., 2001), the antiferromagnetic compound CeRhlns is 
driven towards an instability and incipient superconduc- 
tivity upon application of hydrostatic pressure of pc ~ 
14.5 kbar (Fisher et al, 2002; Hegger et al, 2000). The 
compounds CeTIus consist of nearly two-dimensional 
Celns layers, separated by intercalated TIn2 layers. How- 
ever, the magnetic fluctuation of CeRhIn5 as determined 
from neutron scattering are three-dimensional albeit with 
some anisotropy (Bao et al., 2002). The pressure-induced 
transition to superconductivity appears to be of first or- 
der, thus "avoiding" a QCP. On the other hand, pro- 
nounced deviations from the Fermi-liquid dependence 
of p{T) were reported by Muramatsu et al. (2001) to oc- 
cur over a wide pressure range. Recent studies of the an- 
tiferromagnetic/superconducting phase boundary under 
pressure by Park et al. (2006) and Knebel et al. (2006) 
revealed a quantum critical line between a phase of co- 
existence and a purely superconducting phase. Applying 
additionally a magnetic field, this line is suggested to end 
in a quantum tetracritical point. 

Turning to Celrlns, this material exhibits bulk super- 
conductivity below 0.4 K, with p{T) vanishing already 
at 1.2 K (Petrovic et al., 2001a). The Ce-derived spe- 
cific heat Cm, with lattice and nuclear contributions sub- 
tracted, exhibits when plotted as Cm/T vs. £nT a broad 
hump for magnetic fields applied both along and perpen- 
dicular to the c axis. For fields || c larger than 12 T, 
a divergent Cm/T for T ^ is observed which is inter- 
preted as NFL behavior. Below 0.5 K with field applied 
to suppress superconductivity, the resistivity follows a 
T^-^ dependence, also signaling NFL behavior. These 
findings are interpreted as being due to a field-induced 
metamagnetic transition (Capan et al., 2004). 



5. YbRh2(Sil-a,Gea,)2 

The compound YbRh2Si2 was the first Yb compound 
to show pronounced non-Fermi-liquid effects near a mag- 
netic ordering temperature at very low temperature 
(Trovarelli et al., 2000). It crystallizes in the tetragxj- 
nal I4/mmm structure (same structure as CcCu2Si2). At 
high temperatures T > 200 K the susceptibility follows a 
Curie- Weiss law with /Xeff = 4.5/UB/Yb-atom but differ- 
ent Weiss temperatures 0±^ = — 9K and 0[| = — 180K, 
due to magnetocrystalline anisotropy, where || and _L re- 
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FIG. 15 Yb contribution to the specific-heat coefficient, 
AC/T, of YbRh2Si2 vs. T (logarithmic scale) for differ- 
ent magnetic fields B applied along the a axis. The dot- 
ted straight lino represents AC ~ IniTo/T), the horizon- 
tal line the B = 6T data after subtraction of the hyperflne 
contribution. The inset demonstrates scaling of the data as 
C{B)/T-C{Q)/T = f{B/T'^) with /? = 1.05. Prom TrovareUi 
et al., 2000. 



fer to directions relative to the c axis. At O.IK, the large 
ratio x±/X\\ ~ 100 classifies YbRh2Si2 as an easy-plane 
system, as opposed to the easy-axis system CcCue-sAua;. 

Maxima in the ac susceptibility xTC^) (Trovarelli et 
al, 2000) and specific heat C(T) (Gegenwart et at, 2002) 
as well as a kink in the resistivity (Gegenwart et al., 
2002) around 70 mK signal the onset of antiferromagnetic 
ordering, although to date no neutron scattering data 
are available to corroborate this assignment. (Recently 
Gegenwart et al. (2005) reported evidence for ferromag- 
netic quantum critical fluctuations, see below). Well 
above the magnetic ordering temperature, i.e., between 
0.3 and 10 K, the speciflc heat varies as C/T = aln{To/T) 
with a = 0.17J/molK2 and Tq = 24K (Fig. 15). The 
electrical resistivity was found to vary as p = po + A'T 
between 20 mK and IK, see Fig. 16 (Trovarelli et al., 
2000). (The above-mentioned kink at the ordering tran- 
sition was only found later, likely due to increased sample 
quality.) These NFL features look at flrst sight very sim- 
ilar to those in CeCue-^Aua; and prompted Trovarelli 
et al. (2000) to suggest quasi-2d critical fluctuations by 
analogy with CeCug-ajAua;. Indeed, there is a surpris- 
ing scaling when C/T is plotted as function of £n{To/T): 
the data for YbRh2(Sio.95Ge)o.o5 and CeCu5.8Ago.2 fall 
on top of each other in the range where C/T oc in{To/T) 
(Kuchler et al, 2004). The upturn of C/T below 0.4 K for 
YbRh2Si2 (Fig. 15) can be modelled as C/T ~ T"" with 
a ~ 0.3 (Custers et al., 2003). This upturn suggests the 
presence of an additional low-energy scale below which 
the £nT behavior of C/T is cut off. It is an open question 
whether such crossover exists in CeCue-xAuj, at corre- 
spondingly low reduced temperatures T/Tq where Tq w 
6K is a factor of 4 smaller. Interestingly, the low-T up- 
turn is suppressed below 0.35 K when YbRh2Si2 is diluted 
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FIG. 16 Low-temperature electrical resistivity p of YbRh2Si2 

at p = measured along the a axis as a function of temper- 
ature T, obeying p{T) — po + hT" , with a ~ 1. (a) Temper- 
ature dependence of the effective exponent e = d(:n/S.p/d£nT 
with Ap = p-po. (b) p(T) plotted as p vs. T^ for B < 14T 
applied along the c axis. Prom Trovarelli et al., 2000. 



with La (Ferstl et al, 2005): For Ybo.95Lao.o5Rh2Si2, 
C/T - £n{To/T) is observed down to 0.35 K, the lowest 

T measured. 

Measurements on YbRh2(Sii_a;Gea;)2 with a nominal 
Ge concentration x = 0.05 showed that the magnetic 
ordering is suppressed down to 20 mK where a rather 
broad maximum in the specific heat is observed (Custcrs 
et al., 2003). Above 100 mK, the C/T curves for x = 
and X = 0.05 are nearly identical, indicating that the 
T^" "* upturn is not associated with classical short-range- 
order effects in the proximity to magnetic ordering. 

A further interesting discovery for YbRh2Si2, following 
the prediction by Zhu et al. (2003) discussed in Sec. III.D, 
was the observation of a divergent Griinciscn parameter 
r near the QCP by Kiichlcr et al. (2003). As shown 
in Fig. 17, the P data for YbRh2(Sio.95 Geo. 05)2 follow a 
T~° '' dependence at lowest temperatures, i.e., between 
50 mK and 0.6 K, with a slightly steeper slope approach- 
ing an exponent -1 at higher T (Fig. 17). The T~°-^ 
dependence of F at low T is certainly weaker that 
predicted for the 3d LGW scenario (and actually ob- 
served for CeNi2Ge2 by Kiichler et al. 2003), but also 
weaker than T''^ (eninT) / (inT) expected for a 2d LGW 
model. This observation may support the idea that de- 
grees of freedom other than magnetism become critical 
at the QPT, see Sec. III. I. However, the data would im- 
ply critical exponents with l/izv) = 0.7 (at least if the 
transition is below its upper critical dimension), which 
has not been verified by other measurements. 

Rather detailed experiments were carried out to study 
the behavior of YbRh2(Sii_xGea;)2 in magnetic fields, ex- 
ploring the field timing of a quantum critical point, as 
previously done for the CeCue-aAga; and CeCue-xAux 
series. The magnetic B — T phase diagram of YbRh2Si2 
shown in Fig. 18 is universal with respect to the dif- 
ferent field orientations when the field axes are scaled 



appropriately, similar to the behavior shown before for 
the anisotropic easy-axis system CeCue-iAuj; (Schlager 
et al., 1993). While for high magnetic fields constant 
Kadowaki- Woods ratios A/^q and A/x^ typical for a 
Fermi liquid are observed, the ratio A/^q ~ {B—Bc)~^'^^ 
appears to diverge close to Be, as has been shown for 
YbRh2(Sio.95Geo.o5)2 by Custers et al. (2003). In this 
alloy, B^ = 0.027 T only and 70 diverges as 70 ~ 
[B - Bc)""-^-\ Yet, the coefficient A of the resis- 
tivity diverges more strongly. For a 2d LGW scenario, 
A/^Q diverges much faster, i.e., A/^q ~ 1/{B — B^) with 
logarithmic corrections. 

As discussed in detail in Sec. III.I, one of the cen- 
tral questions for quantum critical points in HFS is 

whether the Kondo effect breaks down at the critical 
point and whether the Fermi volume changes abruptly 
at the second-order transition (Coleman et al., 2001; 
Senthil, Sachdev, and Vojta, 2005; Si, Smith, and In- 
gersent, 1999; Si et al., 2003). Indeed, Paschen et al. 
(2004) found indications of such a behavior in measure- 
ments of the Hall effect across the field-driven transition 
in YbRh2Si2. In Fig. 19 the Hall constant Rh is shown 
for a geometry (see inset) where the field B2 driving the 
QPT is parallel to the current /, and an additional small 
probing field Bx perpendicular to I induces the Hall volt- 
age. As can be seen in Fig. 19, Rh varies quite strongly, 
consistent with a change of the Fermi volume by one elec- 
tron per unit cell (assuming a free-electron single-band 
picture). For lower T the crossover gets sharper, consis- 
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FIG. 17 Electronic specific heat as Cei/T (left axis) and vol- 
ume thermal expansion as —P/T (right axis) vs T (on a log- 
arithmic scale) for YbRh2(Sio.95Geo.o5)2 at B = 0. The solid 
lines indicate £n{To/T) dependencies with To = 30 K and 13 K 
for Cei/T and —P/T, respectively. The dotted line repre- 
sents -p/T = ao + ai/T with ao = 3.4 x 10"^K"^ and oi = 
1.34 X IQ-'^K"^. The inset displays the log-log plot of TcriT) 
with Fcr = Vm/KT ' Pcv/Ccr usiug Kt = 5.3 X 10~^^Pa~^, 
/3„ = p{T) + aoT, and Ccr = Cei(r). The sohd and dotted 
lines represent Per oc 1/T^ with x = 0.7 and x = 1, respec- 
tively. From Kiichler et al., 2003. 
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FIG. 18 Temperature field phase diagram of YbRh2Si2. Full 
and dotted black curves represent the field dependence of the 
Neel temperature Tn and the crossover temperature T* below 
which the resistivity exhibits a Fermi-liquid dependence. 
Symbols present the crossover field where the Hall resistivity 
changes the slope. Light symbols: B = Bi parallel to the c 
axis, dark symbols: crossed-field experiment with the tuning 
field B = B-z perpendicular to the c axis. The relation B = 
Bi = IIB2 reflects the magnetic anisotropy of YbRh2Si2. 
Inset: Width of the transition at Be vs. T, indicating that it 
vanishes as T — » 0. From Paschen et al., 2004. 



tent with a scaling of the half-width with Vt, see dark 
symbols in the inset of Fig. 18. Paschen et al. (2004) 
have also measured the Hall effect in a different geome- 
try, where the field driving the QPT is identical to the 
field which induces the Hall voltage. This geometry is, 
however, more difficult to interpret, as the resulting dif- 
ferential Hall resistivity is expected to jump even for a 
standard SDW quantum critical point. (Note that such 
a jump is expected to trace Tn, opposite to what is ob- 
served in YbRh2Si2.) 

By using fits to the theory of the anomalous Hall ef- 
fect arising from skew-scattering (Fert and Levy, 1987), 




FIG. 19 Hall effect and magnetizations of YbRh2Si2. Left: 
Hall effect as a function of a magnetic field parallel to the 
current. From Paschen et al, 2004. The crossover at the 
critical field Bq = 60 mT gets sharper towards lower T. Right: 
Magnetization as a function of B, arrows indicate the critical 
field. From Gegenwart et al., 2002. 



Paschen et al. (2004) concluded that the former is suf- 
ficiently small to be either subtracted or neglected (de- 
pending on the geometry). In this case the extrapolated 
jump of the Hall coefficient at T = would indeed prove 
rather unambiguously a jump of the Fermi volume at a 
second-order phase transition, certainly a very spectacu- 
lar result. To establish this firmly, it would be important 
to track the sharpening of the crossover in Fig. 19 to- 
wards lower temperature, and check for similar effects in 
other systems as well. When interpreting the crossover 
at finite T one may also have to take into account that 
YbRh2Si2 is almost ferromagnetic (sec below). Even the 
tiny critical field of 60 mT already induces a sizable mag- 
netization of almost O.I/Ub per Yb which is strongly tem- 
perature dependent, see right panel of Fig. 19. Such a 
magnetization could at least in principle lead to large 
changes of the Fermi surface and therefore of the Hall 
effect. The change in Pxy{B) has also been interpreted in 
terms of quantum critical valence fluctuations (Norman, 
2005). 

We now turn to the issue of critical fluctuations in 
YbRh2Si2 near the QPT. Unfortunately, neither elastic 
nor inelastic neutron scattering data exist to date which 
could identify the type of magnetic order and the nature 
of the critical fluctuations. Hence for a microscopic view 
one has to rely on local probes such as NMR and ESR, 
and on macroscopic magnetization data. Ishida et al. 
(2002) reported ^^Si NMR data on aligned single crystals 
of YbRh2Si2 in fields parallel to the easy axis. While the 
nuclear spin-lattice relaxation rate divided by T. \/TiT, 
was found to vary as ^ T^^/^ at small magnetic field of 
B_\_ = 0.15T down to the lowest temperature (~ 50 mK), 
the Knight shift K first also increases towards T ^ 
but turns over to a constant value at T « 200 mK in the 
same field. While K probes the uniform static suscepti- 
bility x'(ct = 0), l/TiT gives information over the imagi- 
nary part of the q-averaged dynamical spin susceptibility 
x"(q, w) at the (quasi-static) frequency corresponding to 
the nuclear Zeeman splitting. The difference between K 
and 1/TiT has been attributed to the presence of finite- 
q critical fluctuations. For higher fields, K and 1/TiT 
tend to a T-independcnt value and track each other, i.e., 
the Korringa relation 1/TiT ~ SK^ is approximately 
fulfilled between 0.25 and 2.4 T. The strong deviation of 
S from the free-electron value So = ■7r?i7^fcs/MB' 
5 w 0.1 as measured at 100 mK, suggests the presence 
of dominant q = 0, i.e., ferromagnetic fluctuations. 

Similarly the approximately constant ratio A/xq dis- 
cussed above and observed for B > 23^ is not easy to 
understand for dominant AFM fluctuations. In addition, 
the strong increase of the Wilson ratio xo /lo by a factor 
of up to 30 compared to the free-electron value "high- 
lights the importance of FM fluctuations in the approach 
to the QCP" (Gegenwart et al., 2005). However, by com- 
paring tendencies to ferro- and antiferromagnetism in a 
Ge-doped sample with an undoped one, Gegenwart et 
al. (2005) also suggested that ferromagnetic tendencies 
are "not directly correlated to the AFM QCP". Recent 
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magncto-striction data, showing the vanishing of several 
crossover energies at the QCP, have also been discussed in 
relation to uniform magnetism (Gegenwart et al., 2007). 

There is a further intriguing observation in YbRh2Si2, 
possibly related to ferromagnetic fluctuations, i.e., the 
observation of an ESR signal by Sichelschmidt et al. 
(2003) that can be ascribed to a Yb^+ resonance. This 
observation is very surprising because the large Kondo 
temperature Tk ~ 25 K is associated with a very high 
spin-fluctuation rate that would according to current 
wisdom render any ESR signal below Tk unobservable. 
However, Sichelschmidt et al. (2003) observe an ESR 
Yb^"'"4/^'^ signal well below Tk that suggests that at least 
60% of the Yh^+ ions in YbRh2Si2 contribute. 

It is tantalizing that the T""-^ dependence of F and the 
T-O-3 dependence of C/T in YbRh2(Sii_^Ge^)2 at very 
low temperatures, but above the ordering temperature 
correspond to a 2d ferromagnetic (z = 3) scenario of a 
conventional LGW QPT. A speculative view could be 
that ferromagnetic planes are the source of the strong FM 
component of the fluctuations that have been discussed 
above, until the magnetic transition intervenes, perhaps 
induced by a weak antiferromagnetic interplane coupling. 



B. Quantum critical beliavior of itinerant transition-metal 
magnets 

Intcrmetallic transition-metal compounds have been 
traditionally viewed as systems close to a magnetic in- 
stability, and the LGW model has been formulated in 
the sclf-consistently rcnormalizcd (SCR) theory of spin 
fluctuations (Sec. IILE) to describe these systems (Lon- 
zarich and Taillefer, 1985; Moriya, 1985). While MnSi 
and ZrZn2 have been outstanding examples under in- 
tense investigation recently - to be discussed below - 
there have been a number of other systems that should 
be mentioned. 

The nearly ferromagnetic metal Ni can be driven to 
static ferromagnetic order by introducing Pd, with a Pd 
critical concentration of x = 0.026. Here, the resistivity 
p follows a T^/^ dependence (Nicklas et al, 1999). This 
T dependence, as well as the Tc ^ (x — Xc)"^/^ depen- 
dence are in accord with the FM LGW model. Likewise, 
the specific heat AC (after subtraction of the phonon 
contribution) shows AC/T - £n{Tn/T) at x^. The fcrro- 
magnct NiaAl has been investigated with resistivity mea- 
surements under pressure up to 10 GPa (Niklowitz et al, 
2005). The critical pressure Pc where Tg ^ is estimated 
to be around 8 GPa. The T-dependent part Ap of the re- 
sistivity does not reveal pronounced critical behavior, the 
exponent a depends on T and varies for 3.2 GPa between 
1.9 and 1.5 at IK and 5K, respectively, and smoothly 
shifts to lower values with increasing p. Below 1 K, the 
data can be approximated by Ap = AT"^, and A exhibits 
a rather sharp maximum at Pc- This behavior has been 
tentatively associated with a first-order transition near 
Pc (Niklowitz et al., 2005). YMn2 is an itinerant antifer- 



romagnct with a first-order transition around Tn ~100K 
(Frcltoft, 1988). On the other hand, Mossbauer measure- 
ments under pressure indicated a continuous reduction 
of the ordered magnetic moment (Block, Abd-Elmeguid, 
and Micklitz, 1994) suggesting a second-order QCP. 

We finally mention experiments on an, albeit complex, 
elemental metal, /3-Mn (Stewart et al., 2002), where the 
resistivity varies as Ap ~ T^/^ between 5 K (the lowest 
temperature measured) and 25 K, suggesting the prox- 
imity to a AFM QCP. The spin fluctuations, however, 
defy a simple interpretation in terms of a unique scaling. 
As already mentioned at the beginning of this chapter, 
space limitations inhibit the discussion of QCP in transi- 
tion metal oxides, such as high-Tc cuprates, manganites, 
or Sr3Ru207. 



1. Cri_xVj, 

Before diving into ferromagnetic systems, let us discuss 
the material Cri_a;Va;, which derives from the itinerant 
aniiferromagnet Cr, and is one of the very few "non- 
Kondo" systems which can be driven through an antifer- 
romagnetic QPT. 

In particular, the Hall effect near the QCP in Cri^j^V^; 
was investigated with concentration tuning (xc = 0.035) 
(Yeh et al., 2002), and by pressure tuning of a sample 
with X = 0.032 close to Xc (Lee et al, 2004). For this 
concentration Tq = 52 K compared to Tc = 311 K for 
pure Cr. While the transition upon concentration varia- 
tion appears very sharp (i.e. almost of first order), using 
a finely spaced pressure tuning reveals it to be contin- 
uous. The decrease of the inverse Hall constant, i?^^, 
upon opening of the spin-density wave gap, tracks the 
increase of p above the paramagnetic background resis- 
tivity. The Hall effect evolution was interpreted in terms 
of an almost nested Fermi surface (Norman et al, 2003). 
The available data on Cri_a,Va; appear to be consistent 
with a 3d LGW scenario of the transition. 



2. MnSi 

The itinerant-electron magnet MnSi orders below Tc = 
29.5 K in a spiral magnetic structure with a long pitch of 
about 175 A. Early on it was realized (Bak and Jensen, 
1980; Nakanishi et ai, 1980) that the helical order arises 
as the cubic B20 structure (space group P2i3) lacks in- 
version symmetry. Therefore the weak spin-orbit interac- 
tions assume a Dzyaloshinsky-Moriya form / S-(VxS)dr 
which twists the ferromagnetic alignment into a helix. 
The ordered moment of about 0.4 /Lts/Mn-atom is sizable 
for an itinerant magnet. With the possible exception of 
high-frequency optical-conductivity experiments (Mena 
et al., 2003), the ground state below Tc seems to be a 
standard three-dimensional weakly spin-polarized Fermi 
liquid (Fawcett et at, 1970; Ishikawa et at, 1985; Taille- 
fer, Lonzarich, and Strange, 1986). MnSi, being almost 
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FIG. 20 Temperature pressure phase diagram of MnSi. 
Above p* , the transition is first order, see Fig. 21 and van- 
ishes at pc. In the shaded region, the resistivity shows an 
anomalous temperature dependence, Ap{T) ~ y^/s 'pjjg 
sets qualitatively show main features of the neutron scattering 
intensity, see text. Prom Pfleiderer et al., 2004. 



a ferromagnet, played an important role in the devel- 
opment of spin-fluctuation theory for complex transition 
metal compounds (Lonzarich and Taillefer, 1985; Moriya, 
1985), see Sec. III.E. 

When pressure is applied (see Fig. 20) the magnetic or- 
der is suppressed (Thompson, Fisk, and Lonzarich, 1989) 
and vanishes above pc = 14.6 kbar. Above this pressure, 
the resistivity shows an anomalous power-law, Ap ~ T^/^ 
(Doiron-Leyraud et al, 2003; Pfleiderer et al, 2001b). It 
has to be emphasized that this power-law holds over al- 
most three decades in temperature from about 6 K down 
to a few mK and over a large pressure range. 

Therefore the question arises whether the anomalous 
resistivity can be associated with critical fluctuations. 
For example, due to the long pitch of the spiral, one 
might expect a quantum critical behavior similar to that 
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PIG. 21 Magnetic susceptibility of MnSi vs. T at different 
pressures (ambient, 1.80, 3.80, 6.90, 8.60, 10.15, 11.25, 12.15, 
13.45, 13.90, 14.45, 15.20, 15.70, and 16.10 kbar going down, 
starting from the top curve at 30 K). Close to the critical 
pressure, x{T — > 0) shows a pronounced jump indicative of a 
first-order phase transition. Prom Pfleiderer et al., 1997. 
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PIG. 22 Scans of the elastic neutron scattering intensity in 
MnSi at 16 kbar perpendicular (left) and parallel (right) to 
the surface of the sphere shown in Pig. 20 [in (110) direction]. 
Prom Pfieiderer et al., 2004. 

of a ferromagnet with Ap T^/^ (Sec.II.E), not too far 
from the observed Ap ~ T^/'^. There are two main argu- 
ments against such an interpretation. First, the anoma- 
lous transport is observed far away from the putative 
quantum critical point. While a pressure of 14.6 kbar is 
sufficient to suppress Tc from almost 30 K down to zero, 
an up to threefold pressure increase is not sufflcient to 
recover Fermi-liquid behavior at temperatures of the or- 
der of 50 mK. (Recent measurements by Pedrazzini et al. 
(2006) indicate, however, that for p > 3pc the exponent 
c;hanges and Fernii-Uquid behavior might be recovered at 
the lowest temperatures.) More relevant is a second ar- 
gument: Susceptibility measurements, see Fig. 21, show 
that close to pc the phase transition is of first order and 
critical fluctuations should be absent. Indeed, Pfleiderer 
et al. (2001b) and Doiron-Leyraud et al. (2003) have tried 
to estimate the consequence of such a first-order transi- 
tion quantitatively and concluded that within a standard 
spin-fluctuation scenario a resistivity should be ob- 
servable below 1 or 2 K, in contradiction to experiment. 
The situation is, however, complicated by the fact that 
not all experimental probes show indications of a strong 
first-order transition. For example, the hysteretic part 
of the signal in elastic neutron scattering (Pfleiderer et 
al., 2004) is rather small close to pc, and on the ordered 
side of the phase diagram the A coeflScient of the T ^ 
resistivity p = po + AT"^ seems to diverge upon approach- 
ing Pc indicative of a second-order transition (possibly of 
percolative type). 

A consistent explanation of the precise nature of the 
quantum phase transition is presently lacking. Theoret- 
ically, it has been argued by Schmalian and Turlakov 
(2004) that strong fluctuations of the direction of the or- 
dering vector generically drive the QPT towards first or- 
der. If quantum critical fluctuations can be ruled out as 
the origin of the NFL behavior, the alternative scenario 
is the existence of an extended new NFL phase which is 
stable against changes of pressure and temperature in a 
sizable region. This would be remarkable for a three- 
dimensional cubic system where disorder effects should 
be negligible due to high-quality single crystals with a 
mean free path of several thousand A. 

Recent elastic neutron scattering experiments by Pflei- 
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derer et al. (2004) have surprisingly shown that even 
above p,, the magnetic order has not vanished com- 
pletely. While for p < pc one observes Bragg reflec- 
tions in the (111) directions with an ordering vector of 
IQI « 0.037 A-i ( increasing under pressure), one finds 
for p ^ Pc and below a crossover temperature To (shown 
in Fig. 20) a strange type of "partial order" : the neutron 
scattering intensity is concentrated on the surface of a 
sphere in reciprocal space with a radius Q » 0.0437 A~^. 
Resolution- limited scans along the radial direction of this 
sphere, see Fig. 22, reveal that the helical order survives 
above pc on length scales of at least ^ 2000 A. The broad 
distribution of the Q vector in tangential direction on 
the surface of the sphere (Fig. 22 and inset of Fig. 20) 
implies the absence of long-range order: the helices have 
lost their directional order. Interestingly, the total inten- 
sity of the elastic neutron signal is comparable to that at 
ambient pressure implying that a sizable fraction of the 
local order survives above Pc- 

Likewise, zero- field ^^Si NMR experiments by Yu et 
al. (2004) show this local order below Tq. The broad- 
ening of the distribution of NMR frequencies shown in 
Fig. 23 is consistent with a picture of intrinsic phase in- 
homogeneities above p* (note that the powdered samples 
used in this experiment may have extra pinning centers). 
Several theoretical attempts have been made to iden- 
tify possible unconventional order-parameter structures 
and lattices of defects (Binz, Vishwanath, and Aji, 2006; 
Fischer, Shah, and Rosch, 2007; Rofiler, Bogdanov, and 
Pflciderer, 2006) which may serve as a starting point to 
understand the observed partial order. The role of dis- 
order and phase inhomogeneities in this extremely clean 
system remains, however, presently unclear. 

The question arises whether the anomalous partial or- 
der or a new QCP, associated with the vanishing of the 
crossover scale Tq (see Fig. 20), can be the origin of the 
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FIG. 23 Distribution of ^^Si NMR frequencies in MnSi for 
pressures p = 0.58, 1.15, 1.23, 1.40, 1.49, and 1.62 GPa (right 
to left). Inset: The corresponding total intensities. From Yu 
et al, 2004. 



anomalous transport properties in MnSi. For example, 
anomalously soft Goldstonc modes of a hclimagnet give 
rise to singular behavior of thermodynamic and trans- 
port quantities (Belitz, Kirkpatrick, and Rosch, 2006a). 
However, for clean systems Ap ~ j'5/2 i.a,ther than the 
observed T^/^ has been predicted (Belitz, Kirkpatrick, 
and Rosch, 2006b) [neglecting, however, the role of spin- 
orbit coupling (Fischer and Rosch, 2004)]. A challenge 
is to understand the role of the scale To above which the 
partial order seems to vanish. Surprisingly, no signatures 
of Tq can be seen in transport, Ap ~ T^/^ both above and 
below Tq. This is unexpected as at Tq, when long-range 
order sets in, the resistivity shows a pronounced kink 
and a strong reduction. It has therefore been speculated 
by Pfleiderer et al. (2004) that the partial order survives 
on intermediate length and time scales and gives rise to 
an anomalous scattering of electrons (e.g. from fluctuat- 
ing topological defects) even above Tq everywhere in the 
shaded area in Fig. 20. Within this interpretation Tq is 
just a freezing temperature below which the partial order 
gets static and is observable by elastic neutron scattering 
or NMR. As both inhomogeneities and fluctuating order 
on intermediate time and length scales arc expected to be 
of importance in other quantum critical systems as well, 
a further understanding of the NFL behavior in MnSi is 
highly desirable. 



3. ZrZna 

ZrZn2 is an itinerant electron ferromagnet. It has a cu- 
bic structure (C15, Fd3m). The small ordered magnetic 
moment /xq = 0.17/XB/Zr-atom contrasts with the large 
Curie- Weiss moment /ieff = l.Q/Us/Zr-atom. This differ- 
ence and the low Curie temperature Tc = 28.5 K classify 
the system as weak itinerant magnet. It has long been 
known that Tc and fio can be suppressed by hydrostatic 
pressure while jUeflf stays practically constant (Huber et 
al., 1975). However, the detailed Tc(p) dependence de- 
pends strongly on the purity of the sample. ZrZn2 has 
long been considered a candidate for p-wave supercon- 
ductivity (Fay and Appel, 1980). Recent experiments re- 
ported superconductivity in very pure ZrZn2 (Pflciderer 
et al, 2001a) which was, however, subsequently shown to 
arise from a surface layer treated by spark erosion (Yel- 
land et at, 2005). Magnetization experiments on very 
pure crystals (RRR w 100) suggest that po and Tc de- 
crease linearly and drop discontinuously at a pressure Pc 
= 16.5 kbar as shown in Fig. 24 (Uhlarz, Pfleiderer, and 
Hayden, 2004). This is in contrast with the continuous 
transition observed as a function of temperature at p = 0. 
De Haas-van Alphen experiments under pressure show 
that the exchange splitting between two nearly spherical 
Fermi-surface sheets that can be identified with minor- 
ity and majority spin sheets, respectively, decreases with 
increasing pressure (Kimura et al., 2004). 

The magnetization data as function of p and B sug- 
gest the existence of two distinct ferromagnetic phases 
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FIG. 24 Pressure dependence (a) of the spontaneously or- 
dered moment M in units oi hb per formula unit and (b) of 
the Curie temperature Tc of ZrZn2. Dashed line in (b) in- 
dicates the lowest temperature measured, (c) (B, P) phase 
diagram at low T with two ferromagnetic phases FMl and 
FM2 and paramagnetic phase PM, separated by phase lines 
MMTl and MMT2. Inset: Continuous decrease of M as a 
function of T. Prom Uhlarz, Pfleiderer, and Hayden, 2004. 



FMl and FM2. This rather comphcated phase diagram 
may be a consequence of a double-peak structure^ in the 
electronic density of states close to Ep (Sandcman et 
al., 2003), which also possibly causes the first-order tran- 
sition near pc. On more general grounds, the transition 
close to the QPT may be generically of first order due to a 
coupling of long-wavelength magnetic modes to particle- 
hole excitations, see Sec. III.H.l. In this scenario (Be- 
litz, Kirkpatrick, and Rollbiihler, 2005) the second-order 
transition at low T and zero field becomes first order 
at a tricritical point, while second-order phase transition 
lines seam a surface of first-order transitions in the three- 
dimensional {p,B,T) space, see Fig. 4 in Sec. III.H.l. 

A systematic study of Zri_2:Nba;Zn2 (Fig. 25) revealed 
a divergence of x for T ^ at a critical concentration 
Xc = 0.083 where x ^ T-^/^ (Sokolov et at, 2006). The 
spontaneous moment vanishes linearly when x Xc- 
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FIG. 25 Temperature dependence of the inverse initial sus- 
ceptibility of Zri_a;Nba;Zn2 indicating the critical behav- 
ior of X near X « 0.08. From Sokolov et al., 2006. 



C. Superconductivity near the magnetic-non-magnetic 
quantum phase transition 

For many years, CeCu2Si2 stood out as the single 
example of superconductivity in HFS, followed by the 
discovery of superconductivity in UPts, UBei3, and 
URu2Si2. CeCu2Ge2 becomes superconducting above a 
pressure of 8GPa (Jaccard, Behnia, and Sierro, 1992; 
Jaccard et ai, 1999). A systematic search by several 
groups, notably Lonzarich and coworkers, for supercon- 
duc;tivity in Ce compounds at the brink of magnetic order 
led to the discovery of several new heavy-fermion super- 
conductors, CePd2Si2 (Grosche et al., 1997; Julian et al., 
1996, 1998; Mathur et al, 1998), CeRh2Si2 (Movshovich 
et al., 1996), Celna (Julian et al., 1996, 1998; Mathur et 
al., 1998), and possibly CeCu2 (Vargoz, Link, and Jac- 
card, 1997). In order to observe superconductivity in 
these systems at the magnetic instability where Tn ^ 0, 
pure samples are often essential. This, together with 
the very fact that superconductivity appears at a point 
where low-lying magnetic fluctuations abound, suggest 
that the superconductivity may be magnetically medi- 
ated (Mathur et al, 1998). 

Fig. 26 shows the phase diagram for CePd2Si2 where 
Tn drops from 10.5 K at ambient pressure to below 1.6 K 
around 25 kbar. Tn extrapolates to zero by continuing 
the linear T^ip) dependence at Pc = 27 kbar. Around this 
critical pressure the superconducting transition temper- 
ature Tc has its maximum of « 0.6 K and extends almost 
symmetrically to ± 5 kbar around Pc- However, it is not 
clear whether antiferromagnetism gives way to supercon- 
ductivity just below Pc or both types of order coexist. 

At Pc the electrical resistivity p varies with T^-^ over 
almost two orders of magnitude up to T > 30 K (see in- 
set of Fig. 26). This quasi-linear T dependence of p, 
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together with the hnear relationship Tn ^ \p ^ Pc\^i 
ip = 1, have led to the suggestion of 2d fluctuations 
(Mathur et al., 1998), as in CeCue-xAux- On the other 
hand, the anomalous T^-^ dependence may arise from 
an effective crossover of p{T) in the 3d spin-fluctuation 
scenario (Rosch, 1999, 2000) as discussed in Sec. III.F. 
CeNi2Ge2 exhibits a similar anomalous power law of p{T) 
and a C/T = 7o — PVt dependence at ambient pressure 
(Grosche et al, 1997; Julian et al, 1998; Kiichler et al, 
2003). Some very pure samples even exhibit traces of su- 
perconductivity at p = (Gegenwart et al, 1999). These 
findings suggest that CcNi2Gc2 at p = is right at the 
magnetic instability, although other samples become su- 
perconducting at p > 15kbar only (Grosche et al, 1997). 

It should bo mentioned that besides systems exhibit- 
ing a narrow dome of superconductivity near a magnetic 
QCP, exemplified by CePd2Si2, wide pressure ranges of 
superconductivity are observed, e.g. in CeCu2Ge2. In- 
deed, it has been suggested that superconductivity in 
CeCu2Ge2 is mediated by valence fluctuations rather 
than spin fluctuations which may cause a rather wide 
range of superconductivity. This idea is corroborated by 
recent experiments on CeCu2Si2 doped with Ge where 
indeed two disconnected superconducting regions are ob- 
served as a function of pressure as shown in Fig. 27 (Yuan 
et al, 2003, 2006). The relation to quantum criticality is 
manifest through the pronounced pressure dependence of 
the residual resistivity po and of the resistivity exponent 
a. 

Celns has been suggested to be an antiferromagnet 
with 3d critical fluctuations, giving way to superconduc- 
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FIG. 26 Temperature-pressure phase diagram of a high- 
purity CePd2Si2 single crystal. Superconductivity appears 
below Tc in a narrow window where the Neel temperature Tn 

tends to absolute zero. For clarity, the values of Tc have boon 
scaled by a factor of three, and the origin of the inset has been 
set below absolute zero. The inset shows that the resistivity 
p exhibits a T^'^ dependence on temperature over a wide T 
range. From Mathur et al, 1998. 



tivity again below 1 K under pressure, on account of the 
T^-^ dependence of A/7(T) close to ~ 25 kbar (Julian 
et al, 1998; Mathur et al, 1998). In hue with this in- 
terpretation, the dependence of Tn on |p — Pc| appears 
to be sub-linear; an exponent ^/^ = 2/3 is predicted for a 
3d antiferromagnet, Eq. (94). However, '^^^In NQR mea- 
surements around the critical pressure observed no trace 
of NFL spin fluctuations, perhaps indicating a first-order 
magnetic transition (Kawasaki et al, 2004). 

The Celns-derived CeTIn5 samples exhibit supercon- 
ductivity over a wide range in alloying among each other 
(T = Co, Ir, Rh), as shown in Fig. 28. Particularly in- 
teresting is CeCoIns (see also Sec. IV. A. 4), where a mag- 
netic field induces a second superconducting phase below 
Bc2 (Bianchi et al, 2003a). This phase is a prime can- 
didate for a modulated superconducting state, originally 
proposed by Fulde and FerrcU (1964) and Larkin and 
Ovchinnikov (1965), dubbed FFLO state. Last not least, 
we wish to mention CeRh2Si2 which becomes supercon- 
ducting above Pc = 6 kbar (Movshovich et al, 1996), how- 
ever, here antiferromagnetism disappears through a first- 
order transition, see Sec. IV.A.2. 

A few general remarks are in order. Although the var- 
ious Ce compounds have vastly different magnetic order- 
ing temperatures, the superconducting Tc are quite simi- 
lar of the order of 0.5 K with the exception of the "high" 
Tc = 2.3 K for CeCoIn5. Further, the diverse role of im- 
purities in HFS superconductors is not understood: while 
some of them arc extremely sensitive to impurities, e.g. 
CePd2Si2, others are not. In most of the HFS, the ques- 
tion of whether magnetism and superconductivity coex- 
ist cooperatively has not been investigated in detail. A 
prominent counterexample is CeCu2Si2 where supercon- 
ductivity competes with A-phase magnetism. We finally 
note that it has been speculated by Steglich (2005) that 




2 
p-pjv.) (GPa) 

FIG. 27 Combined {p,x,T) diagram of CeCu2(Sii_a;Ge:c)2. 
Pci(x) denotes the prcssuro-tunod QCP where the Nccl tem- 
perature Tn 0; Pci(O) = 0, pci(2.5) = 2.4 GPa. Open sym- 
bols denote Tn, closed symbols the superconducting Tc. In 
the wide shaded area, non-Fermi-liquid behavior is observed 
in the electrical resistivity. From Yuan et al, 2006. 
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FIG. 28 Phase diagram of Ce(Rh,Ir,Co)In5 displaying the in- 
terplay of antiferromagnetic order (AFM) and superconduc- 
tivity (SC). Prom Pagliuso et al, 2002. 

conventional LGW critical points favor unconventional 
superconductivity, while non-LGW QCP disfavor super- 
conductivity in their vicinity. 

While all the above systems exhibit superconductiv- 
ity in the vicinity of antiferromagnetic instabilities, there 
has been a long-standing quest for superconductivity of 
an incipient ferromagnet (Fay and Appel, 1980). Ini- 
tially such type of superconductivity, presumably of spin- 
parallel pairing, was sought among weak itinerant ferro- 
magnets, e.g., ZrZn2. However, the first unambiguous 
observation of a ferromagnetic superconductor was found 
in UGe2 under pressure close to the FM instability, with 
a maximum Tc ~ 0.4 K (Saxena et al., 2000), while the 
related system URhGe becomes superconducting even 
at ambient pressure (Tc ~ 0.3 K) (Aoki et al, 2001). 
Fig. 29 shows the pressure-temperature phase diagram 
of UGe2 as investigated by Pfleiderer and Huxley (2002), 
revealing two different ferromagnetic phases, overall sim- 
ilar to ZrZn2 (Fig. 24). The different phases are clearly 
identified by their different ordered moments. As in 
ZrZn2, the pressure-driven transitions between the FM 
phases and from ferromagnet to paramagnet are clearly 
first order. Consequently, the electrical resistivity shows 
FL-like behavior across the transition, with Ap = AT^ 
and A peaking at the critical pressure Pc (Saxena et al., 
2000). The superconductive phase in UGe2 resides com- 
pletely inside the FM domain, giving strong evidence for 
spin-parallel pairing in coexistence with fcrromagnctism. 
These observations have instigated numerous theoretical 
studies (among others Chubukov et al. 2003; Kirkpatrick 
et al. 2001; Kirkpatrick and Belitz 2003; Roussev and 
Millis 2001), since the early weak-coupling calculation 
(Fay and Appel, 1980) suggested two superconductive 
"domes" just below and above Pc- 

Of particular interest is the observation of two distinct 
superconductive phases in URhGe as a function of mag- 
netic field (Levy et al., 2005). The low-field phase sub- 
sides at a critical field of Bc2 = 2 T applied along the 
b axis of the orthorhombic crystal structure. In a field 
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PIG. 29 Phase diagrams of UGe2. (a) Pressure-temperature 
phase diagram indicating the Curie temperature Tc and the 

transition temperature Tx between two ferromagnetic phases 
PMl and FM2. Shaded area is the superconductive transi- 
tion temperature Ts. (b) Pressure dependence of the sponta- 
neously ordered moment M in units of /is per formula unit at 
2.3 K. (c) Pressure-field phase diagram at T = 2.3 K indicat- 
ing the metamagnetic transition and the transition between 
FMl and FM2. From Pfleiderer and Huxley, 2002. 

of 11.7 T the spin direction changes abruptly. A sec- 
ond superconductive phase appears in the vicinity of this 
transition, extending between 8 and 13 T. Hence, view- 
ing the spin reorientation as a quantum phase transition, 
the siiperconductivity is clearly linked to this QPT. 

It has been suggested that the absence of supercon- 
ductivity in MnSi, even in very pure samples, is due to 
the lack of inversion symmetry in this material (Sax- 
ena et al., 2000). However, superconductivity below 
Tc = 0.75 K has been observed recently in a HFS lacking 
inversion symmetry, CePtsSi, albeit in coexistence with 
anMerromagnetism, with Tc « 2.2 K (Bauer et al., 2004). 

V. COIMCLUSIOIMS 

In general, metallic systems of interacting fermions at 
low temperature are well described by the Landau Fermi- 
liquid theory. Non-Fermi-liquid behavior over an ex- 
tended range of temperatures down to absolute zero may 
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occur near the borderline between two qualitatively dif- 
ferent ground states. The order-parameter fluctuations 
in the neighborhood of such a quantum critical point can 
induce singular scattering between the fermions, leading 
to a breakdown of the usual phase-space arguments on 
which Fermi-liquid theory is based. 

In this review we aimed to give a systematic, balanced 
and critical account of the present knowledge in the area 
of Fermi-liquid instabilities near quantum critical points. 
We focused on magnetic QCP, which are perhaps the best 
studied cases at present. While a large body of experi- 
mental and theoretical work on such systems has been 
accumulated over the past 20 years, several principal 
questions appear to be far from answered. Theoretically 
well established examples of non-Fermi-liquid behavior 
induced by quantum critical fluctuations are provided 
by single-site quantum impurity models, a prominent ex- 
ample being the multi-channel Kondo model. Unfortu- 
nately, on the experimental side a definite realization of 
this model in a bulk material is yet to be confirmed. 

For lattice systems, the best understood situations are 
those with only a single variable relevant for the crit- 
ical behavior, the order parameter. Here, a Landau- 
Ginzburg- Wilson field-theoretical description allows one 
to calculate critical exponents and scaling functions, and 
a self-consistent RPA-type theory provides in addition 
approximations for the full dependence of observables on 
parameters like pressure, temperature, or field, as shown 
in the pioneering works of Hertz and Moriya, extended 
by Millis. In general, different universality classes are ex- 
pected for ferromagnets and antiferromagnets. Further, 
one should distinguish nominally clean systems, where 
the critical thermodynamics is not influenced by disorder, 
from disordered ones. (The latter case is more compli- 
cated, as disorder may modify critical exponents or even 
destroy the zero-temperature phase transition.) Experi- 
mentally well characterized examples of systems follow- 
ing the LGW predictions for clean systems appear to be 
CeNi2Ge2 and Cei_2;Laa;Ru2Si2, which are close to an 
antiferromagnetic instability. 

Somewhat unexpectedly, a growing number of systems 
has been discovered showing properties inconsistent with 
those of LGW theory. A common theme is that further 
soft variables exist such that the LGW approach is no 
longer applicable. These soft variables may be either the 
fermionic particle-hole excitations which strongly cou- 
ple to the order parameter (as is the case for ferromag- 
nets), or additional degrees of freedom like those asso- 
ciated with the Kondo effect in heavy-fermion systems. 
A full understanding of these more complex situations 
of coupled slow modes is not available at present, al- 
though several interesting proposals exist and have been 
reviewed above. 

In the following we try to summarize what we think 
are pressing questions in the field, starting with the 
theory side. (Tl) Can one formulate a theory for the 
critical behavior in metallic magnets for the approach 
from the ordered side, along the works of Hertz and 



Moriya? (T2) Can one analyze a coupled theory of order- 
parameter fluctuations and fermions for clean metallic 
magnets using RG techniques? (T3) Can one prove that 
the LGW theory for 3d antiferromagnets is stable and 
self-consistent? Is there room for non-LGW criticality 
in 3d antiferromagnets? (T4) What are the character- 
istics of the metallic antiferromagnetic QPT in d = 2? 
(T5) Can one develop one of the scenarios of the break- 
down of the Kondo effect (Sec. III. I) to consistently de- 
scribe the phenomenology of materials like CeCug-j^Aua, 
or YbRh2Si2? (T6) Can one understand in more detail 
the physics of QPT smeared by disorder? (T7) Are there 
scenarios for stable NFL phases in 3d systems, which may 
explain, e.g., the physics of MnSi? 

Important directions for experimental work are: (El) 
A detailed analysis of energy- and momentum-resolved 
magnetic fluctuations of materials with non-LGW criti- 
cality (in addition to CeCu6_a;Aux) would be instructive. 
Is there w/T scaling? Arc quasi-2d fluctuations generic? 
(E2) Given a nominally clean material with a well- 
characterized QCP, it would be important to systemati- 
cally study the influence of disorder on the critical prop- 
erties. (E3) Can one find a "dirty ferromagnet" which fol- 
lows the predictions of the theory in Sec. III.H.l (which is 
the only well-developed theory of order-parameter fluc- 
tuations coupled to fermionic modes)? (E4) Can one 
single out heavy-fermion systems showing clear-cut ev- 
idence for a jump in the Fermi volume as T — > 0? 
(E5) Can one identify two or more distinct diverging 
time/length scales near certain magnetic heavy-fermion 
critical points? (E6) Can one establish a link between the 
presence or absence of superconductivity near a QCP and 
the universality class of the QCP? 
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